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Hydromechanische Grundlagen fiir die Behandlung von stationiren 
und instationaéren Grundwasserstrémungen. 


Von G. Heinrich und K. Desoyer. 


1. Einleitung. In den Lehrbiichern der Hydraulik geht man bei der Behandlung der 
Grundwasserbewegung im allgemeinen vom Darcy-Gesetz aus, ohne auf eine genaue Analyse 
der auf das Wasser und den Festkérper wirkenden Krafte einzugehen. Dies fiihrt dazu, daB 
bei den verschiedenen Autoren ganz verschiedene Ansichten tiber den ,,Strémungsdruck“ 
herrschen, obwohl die Hydromechanik klare und eindeutig definierte Begriffe zur Erfassung 
des hier vorliegenden Kraftespieles zur Verfiigung stellt. Auch werden die bei der Ableitung 
der Gleichungen fiir die Grundwasserstrémung zugrunde gelegten vereinfachenden Annahmen 
meist nicht geniigend genau prazisiert. Inshesondere werden bei der Untersuchung insta- 
tionarer Grundwasserstrémungen gewisse Vereinfachungen eingefiihrt, die in vielen Fallen 
praktisch nicht zutreffen. Man kommt jedoch ohne diese Vereinfachungen zu einer bedeutend 
ubersichtlicheren und theoretisch weitaus mehr befriedigenden Lésung der Probleme der 
instationéren Grundwasserstrémung. 

Es soll daher im folgenden nach einer Analyse des bei der Grundwasserstrémung auftreten- 
den Kraftespiels eine verallgemeinerte Theorie der Grundwasserstrémung, insbesondere der 
Spiegelbewegung bei instationéren Strémungen entwickelt werden. Eine Auswertung mittels 
der Relaxationsmethode soll in einer folgenden Arbeit an praktischen Beispielen gezeigt werden. 
Eine genauere Kenntnis der instationaren Strémungszustande ist z. B. fiir die Beurteilung der 
Rutschgefahr von Dammen sowie fiir die Abschatzung der Ergiebigkeit von Brunnenanlagen 
von praktischer Bedeutung. 


2. Analyse der wirkenden Krafte und deren Gleichgewicht. Um zu einer methodisch ein- 
wandfreien Darstellung der Grundwasserstrémung zu gelangen, ist es wesentlich, zwischen 
Mikro- und Makrostrémung zu unterscheiden. Unter Mikrostrémung soll die tatsachlich auf- 
tretende Bewegung der Fliissigkeit in dem unregelmafigen Gefiige des Festkérpers verstanden 
werden. Diese Mikrostrémung ist nur insoweit von Interesse, als sie durch geeignete Mittel- 
wertshildungen durch eine Makrostrémung ersetzt werden kann. Um zu geeigneten Mittelwerts- 
bildungen zu gelangen, miissen die fiir die Untersuchung der Makrostrémung zugrunde geleg- 
ten Volum- und Oberflachenelemente einerseits groB genug gewahlt werden, damit eine Mitte- 
lung ohne unzulassig grofe Streuungen mdglich wird, andererseits klein genug gehalten 
werden, um die Abhangigkeit der durch Mittelung gewonnenen Gréfen vom Ort wieder- 
zugeben. Es mu dabei vorausgesetzt werden, das die Volumina der Einzelkérner des durch- 
strémten Festkiérpers tber das zugrunde gelegte Volumelement statistisch verteilt und 
geniigend klein sind gegentiber dem Gesamtvolumen des betrachteten Bereiches, so da die 
obigen Mittelwertsbildungen zulassig sind. 

Es ist wesentlich, eine genaue Unterscheidung zu treffen zwischen den Kraften, die auf die 
Fliissigkeit, und jenen, die auf den Festkérper wirken. 

Es wird angenommen: 1. Die Korner des Festkérpers, die an sich beliebige Formen aut- 
weisen kénnen, berihren einander praktisch punktformig!. 2. Der Festkérper befindet sich 
dauernd in Ruhe. 3. Die Flissigkeit sei praktisch inkompressibel und homogen. 4. Die in der 
Mikrostrémung auftretenden Kapillarkrafte werden nicht beriicksichtigt, es wird also der sog. 
Kapillarsaum aufBer acht gelassen. Desgleichen sollen andere Krafte molekularer Natur (z. B. 
elektrische) wie sie in den feindispersen Anteilen des Bodens auftreten kénnen *, nicht in 

* Betracht gezogen werden. 5. Im Inneren der Strémung soll die Flissigkeit die Poren des Fest- 


1 Hine Erweiterung auf allgemeine Festkérperskelette soll in einem nachfolgenden erganzenden Aufsatz 


durchgefiihrt werden. 
2 Vel. J. Kozeny, Hydraulik, S. 381, Wien 1953. 
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kérpers vollstindig erfillen. 6. Die Filterwirkung des Festkérpers sei geniigend grof, so dah 


eine ,,quasistationare Strémung“, bei der die Tragheitskrafte in der Flissigkeit gegeniiber den | 
anderen Kraften vernachlissigt werden kénnen, bereits bei kleinen Strémungsgeschwindig- — 


keiten auftritt. Bei gewisseni. a. sehr kurz dauernden Anlaufzustanden darf man die Tragheits- 
krafte nicht mehr vernachlassigen, (vgl. Abschn.5). 7. Eine Verdunstung der Flissigkeit wird 
nicht beriicksichtigt. 

Wir denken uns nun aus dem System Festkérper und Fliissigkeit (es werde im folgenden 
immer als ,,Gemisch“ bezeichnet) ein beliebiges endliches Volumen herausgegriffen, aber so, 
daB die Begrenzungsfliche kein Korn des Festkérpers durchschneidet. Es sollen nun fiir den 
auf diese Art abgegrenzten Flissigkeitskérper die Bewegungsgleichungen und fiir den zugehori- 
gen Teil des Festkérpers die Gleichgewichtsbedingungen aufgestellt werden. In der makro- 
skopischen Betrachtung teilen wir die auf den abgegrenzten Fliissigkeitskérper wirkenden 
Krafte ein in Oberflachenkrafte, Volumkrafte und Quasivolumkrafte. Die Oberflachenkrafte 
wirken an der Begrenzungsflache des abgegrenzten Volumens. Da zufolge der statistischen 
Verteilung der Kérner der Quergradient der Geschwindigkeit der Mikrostrémung in Richtung 
der auSeren Normalen des Flachenelementes mit gleicher Haufigkeit positiv und negativ sein 
wird, werden sich bei der makroskopischen Betrachtung die Schubspannungen iiber ein Flachen- 
element gerade aufheben. Fiir die Makrostrémung sind daher die an der Begrenzungsflache 
wirkenden Oberflachenkrafte reine Druckkrafte. Da die Begrenzungsflache so gelegt wurde, 
da8 sie kein Korn durchschneidet, wirkt nach Voraussetzung 1 dieser Druck auf die gesamte 
Schnittflache. Als Volumkraft kommt praktisch nur die Schwerkraft in Betracht. Die in der 
Mikrostrémung von den Kérnern auf die Fliissigkeit ausgeiitbten Krafte werden in der Makro- 
strémung als Quasivolumkrafte iitber das ganze Volumen verteilt gedacht. Entsprechend der 
Aufteilung dieser Krafte in der Mikrostrémung in Normaldricke und Reibungen sollen auch die 
Quasivolumkrafte in der Makrostrémung in Normal- und Reibungskrafte unterteilt werden. 
Im folgenden sei do der nach aufen positiv gezahlte Vektor des Oberflachenelements der 
Begrenzungsflache, dV ein Volumelement des Gemisches, n das auf die Volumeinheit des Ge- 
misches bezogene Porenvolumen1?, das i. a. eine Funktion des Ortes sein kann, aber nach 
Voraussetzung 2 keine Funktion der Zeit, g die Fallbeschleunigung, y,, das spezifische Gewicht 
der Fliissigkeit, p der in der Flissigkeit herrschende Uberdruck iiber den Atmospharendruck, 
U das Potential der Schwerkraft bezogen auf die Masseneinheit der Flissigkeit, ‘hi, die von der 
Hautreibung an den Kérnern herrihrende, auf die Flissigkeit wirkende Quasivolumkraft, 
bezogen auf die Volumeinheit des Gemisches, \7 der bekannte Vektoroperator. Mit diesen 
Bezeichnungen ergibt sich die Resultante der Oberflachenkrafte zu — |p do, die der Schwer- 


kraft zu — + n \/U dV, die von der Hautreibung herrihrende Quasivolumkraft zu ih Ry dV 


und die von den Normaldriicken der Kérner auf die Flissigkeit herrithrende Quasivolumkraft 
zu | (1 — n) \/p dV, wobei die Integrationen iiber die Oberflache bzw. das Volumen des heraus- 
gegriffenen Gemischkérpers zu erstrecken sind. 

Der Ausdruck i (1 — n) \/p 4V ergibt sich auf folgende Weise: in einem Volumelement dV 
der Makrostrémung mégen sich im Mittel N Kérner befinden. Die Resultante der von diesen 


N 
Kérnern auf die Flissigkeit ausgetitbten Druckkrafte ist dann ») f p do; wobei do; einen nach 
j=1 


. . - . I 
auSfen positiv gezahlten Vektor eines Oberflachenelementes des j-ten Kornes bedeutet. Nach 
Anwendung des Gaufschen Satzes erhalt man daraus 


Nee N 
=) Vp 4Vx; = Vp & Vx;; 
diz i= 

wobei dVx; das Volumelement des j-ten Kornes und Vx; das ganze Volumen des j-ten Kornes 
bedeutet. Die obige Gleichsetzung ist zulassig, da der makroskopische Wert \/p tber den 
Bereich dV als konstant angesehen werden kann. Ferner gilt die Beziehung 


* Bei fein verteilten Béden kommt hierfiir nur das sog. ,,spannungslose Porenvolumen“ in Frage 
das kleiner ist als das wahre Porenvolumen, vgl. Kozeny a.a.0., S.381. Das an die Kérner gebun- 
Bae hygroskopische Wasser soll in diesem Fall als zum System des Festkérpers gehérig angesehen 
werden, 
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Aus diesen Gleichungen folgt demnach 
N 
> f p do; =(1—n) VpdvV. 
jai 


Da bei Vernachlassigung der Tragheitskrafte die angegebenen Krafte fiir jedes beliebig 
herausgegriffene Volumen im Gleichgewicht stehen miissen, ergibt sich fiir die Flissigkeit die 
Beziehung 


aM |p ig | Zen VU dV + [ Rea + | a ay aves. 


fo} 


Durch Anwendung des Gaufschen Satzes auf das erste Glied erhalt man 


Wg Vp ene 5h VU+R, + (1 — n) VP] dV = 0% 


Da diese Gleichung fiir jeden beliebigen Integrationsbereich gelten mu8, so folgt schlieBlich 


Seen ihe U 4 P)- (1) 


Man sieht aus (1), da das Feld der Krafte ,, bei konstantem Porenvolumen n ein Potential- 


feld ist. 


Die auf den abgegrenzten Teil des Festkérpers wirkenden Krafte sollen ebenfalls in Ober- 
flachenkrafte, Volumkrafte und Quasivolumkrafte eingeteilt werden. Die Oberflachenkrafte 
erhalt man durch Mittelung der Beriihrungskrafte zwischen den Kérnern, die beiderseits eines 
Oberflachenelements von geeigneter GréBe liegen. Als Volumkraft kommt praktisch nur die 
auf die Kérner wirkende Schwerkraft in Frage. Die Reaktionskrafte der auf die Flissigkeit 
wirkenden Quasivolumkrafte (Normal- und Reibungskraft) wirken auf den Festkérper zuriick. 


Dai. a. von Korn zu Korn nicht nur Normaldricke, sondern auch Reibungs- und Kohasions- 
krafte ubertragen werden, bilden die tiber ein Oberflachenelement gemittelten Korn-zu-Korn- 
Krafte ein allgemeines Spannungsfeld, das durch einen symmetrischen Spannungsaffinor ©x 
gekennzeichnet werden soll. Es ist dann! ©x- dp der durch das Oberflachenelement do iiber- 
tragene Kraftvektor. { x: do bedeutet die Resultante der auf die Begrenzungsflache des 
herausgeschnittenen Festkérperteiles wirkenden Oberflachenkrafte (Korn-zu-Korn-Krafte). 
Die Resultante der auf den herausgeschnittenen Festkérper wirkenden Schwerkrafte ist 

LB) VK 

8 
das vom Ort abhangen kann, nach Vor. 2 aber nicht von der Zeit. Die beiden Anteile der 
Quasivolumkrafte sind bis auf das Vorzeichen dieselben, wie die auf die Flissigkeit wirkenden. 
Die Gleichgewichtsbedingung fir den herausgeschnittenen Festkérperteil lautet demnach 


[ Sx- ao lbecgerecae VuUdV — [ Re av — | a —n)VpdV=0. 


\/ U dV, worin yx das mittlere spezifische Gewicht des Kornmaterials bedeutet?, 


Durch Umformung des ersten Gliedes mit Hilfe des Gaufschen Satzes erhalt man 


[|v-& care yu Re (1 n) Vp|dv = 0. 


Da diese Gleichung fiir jedes beliebige Volumen gelten muf, so folgt 


V:Gx = Sve TU LR, + (1 —n) Vp. (2) 
Aus (1) und (2) ergibt sich schlieBlich 
Se. ep eee Yo) VU | Wipe (3) 
& 


1 Es bedeuten im folgenden das Zeichen - das skalare Produkt, -- das doppeltskalare Produkt und - 
das Zeichen>» das dyadische Produkt. 
2 Bei feinverteilten Béden wire unter yx das mittlere spezifische Gewicht von Kornmaterial samt 


dem gebundenen Wasser zu verstehen (vgl. FuBnote 1, 5. 74). os 
6 
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Fir den Spezialfall, daB yx und n ortsunabhangig sind, lat sich (3) auf die Form bringen 
Vx beers pe p| 6} =o, (3a) 
& 


worin (© den Einheitsaffinor bedeutet. 


3. Das verallgemeinerte Filtergesetz. Es soll nun ein allgemeiner Zusammenhang zwischen 
der Reibungskraft ,, und der Bewegung der Fliissigkeit hergestellt werden. Zu diesem Zweck 
fiihren wir den Begriff der Filtergeschwindigkeit ein. Fir ihre Definition denken wir uns aus 
dem durchstrémten Festkérper ein kleines Tetraeder (Abb. 1) 
herausgeschnitten, wobei die Achsen x, y,z ein rechtwinkeliges 
Koordinatensystem bilden. Wir bezeichnen mit dD den nach 
auBen positiv gezdhlten Vektor des schragen Flachenelements. 
Seine vektoriellen Komponenten dO,, dD,, dD, nach den drei 
Achsenrichtungen sind zugleich die nach innen positiv ge- 
zahlten Vektoren der in den drei Koordinatenebenen liegen- 
den Begrenzungsflachen. Wir bezeichnen mit dq, dt das 
wahrend des Zeitelements dt durch die Flache d hindurch- 
tretende Fliissigkeitsvolumen, analog seien dq, dt, dq, dt, 
dq, dt die durch die in den Koordinatenebenen liegenden 
Flachen wahrend dt eintretenden Flissigkeitsvolumina. Nach 
Voraussetzung 3 gilt dann 


Abb. 1. Zur Einfiithrung des Vektors 


der Filtergeschwindigkeit. 
dq, = dq, + dqy + dq, . 


Bezeichnen nun v,, vy, v, die auf die Flacheneinheit bezogenen Werte von dq,, dq,, dq,, so kann 
man auch schreiben 


dqy = |€Dz| v2 + |dDy| vy + |dO,| ¥. =0- dO, 2 


wobei ) den aus v,, vy, v, gebildeten Vektor bedeutet. Da sich also dieser Vektor wie ein Ge- 
schwindigkeitsvektor verhalt, soll er als Filtergeschwindigkeit bezeichnet werden. Man sieht, 
dafi eine derartige Definition einer Filtergeschwindigkeit nur bei inkompressiblen Medien 
méglich ist. 

Wir beschranken uns nun im folgenden auf solche Falle, bei denen ein linearer Zusammen- 
hang zwischen der Reibungskraft h,, und der Filtergeschwindigkeit » besteht. Der allge- 
meinste lineare Zusammenhang zwischen den Vektorenfeldern t,, und ) wird vermittelt durch 
einen Affinor 2{. Diesen Zusammenhang schreiben wir in der Form 


Jip = — Wd. (5) 


Damit ist die Méglichkeit der Beriicksichtigung einer Anisotropie des Bodens (z. B. Schichtung) 
gegeben. Aus (1) und (5) folgt 


7 oe EPs eat (6) 
oder 
LRA sagt 
| AG a Wb 
mit 
one 2 
eg = (7) 


YW soll als Widerstandsaffinor bezeichnet werden. Eine andere Form von (6) erhalt man durch 
linksseitige Multiplikation mit dem Affinor 

W-1 — K (8) 
Es ergibt sich 

U 
pee Qe eee 

& sei als Durchlassigkeitsaffinor bezeichnet; er wird im allgemeinsten Fall ein unsymmetrischer 
ortsabhangiger Affinor sein, Wenn langs einer Flache zwei Festkérper aus verschiedenem 
Material aneinandergrenzen, andert sich dort der Durchlassigkeitsaffinor & unstetig. Wir 
fassen § deshalb als eine stiickweise stetige Funktion des Ortes auf. 
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& habe in bezug auf ein beliebiges cartesisches Koordinatensystem (x, y, z) die Matrizen- 


darstellung 
bes kyy | i 
W=[hys  kyy i) (10) 
et a 
Wir setzen 
—+2=h (11) 
& Vw 


und nennen h die Standrohrspiegelhéhe, da sie die Spiegelhéhe in einem an der betreffenden 
Stelle errichteten Standrohr, gemessen von einem willkiirlichen Nullniveau (U = 0) angibt 


(Abb. 2). Aus (9) und (11) ergibt sich 


Z 
b= —&-Vh. (12) 
Fihrt man den Vektor des Standrohrspiegelgefalles 
FS =- Gh (12a) 
ein, so schreibt sich (12) 
b=R- J. (12b) 


Dies ist als eine Verallgemeinerung des bekannten 
Darcyschen Filtergesetzes anzusehen. 

Wenn im Festkérper keine Richtungsanisotropie vor- 
liegt, kann man § = k © setzen, dann spezialisiert sich 
(12) zu der Form 


b=—kVh=kJ, (12c) 


wobei k der in der Hydraulik als Durchlassigkeit be- 
zeichnete Skalar ist 1. Bezeichnet man mit v den Be- 


Abb. 2. Zur Definition der Standrohrspiegelhéhe. 


trag von ) und mit J den Betrag von J, dann erhalt man aus (12c) die bekannte Form 
des Darcy-Gesetzes 


een ie (12d) 


Aus (12c) ersieht man, daB in dem Spezialfall eines ortsunabhingigen k die GréBe h die Rolle 
eines Geschwindigkeitspotentials ibernimmt. 


4, Die Gleichung fiir die Fliissigkeitshewegung. Um zu einer Differentialgleichung fiir die 
Flissigkeitsbewegung zu gelangen, gehen wir aus vom verallgemeinerten Filtergesetz in der 
Form (12). Da das Vektorfeld der Filtergeschwindigkeit » quellfrei ist, gilt 


ioe weal (13) 
Aus (12) und (13) folgt 
Vi (&- VA) =O. (14) 
Dies ist die Differentialgleichung, der die Standrohrspiegelhéhe h geniigen muf. Eine einfache 
Umformung von (14) liefert 


(Vest) (Vb) PRY, V) k=O. (15) 
Fihrt man die in (15) angegebenen Operationen aus, so erhalt man die Differentialgleichung 


Oke | Okyx , Oh,x\ Oh Okay , Okyy | Okyy\ Oh , (Oke , Ohyy | Okgs\ Oh 
(3 ; rane a) ae Creer Verge dy | \ ax + dy | Oz) dz 


Ox 


orh 0rh 07h ; OPh orh 
- hixn 55 =f YY Oy? k.. Oz2 (kez y t ky x) ax dy ‘e (ky: k, y) dy 0 


Orh 
a (kez x at ky z) anos = 0. 


1 Vgl. z. B. Ph. Forchheimer, Hydraulik, S. 53, 3. Aufl., Leipzig und Berlin 1930. J. Kozeny, a, a. O. 
S. 389. 
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Fir den praktisch wichtigen Spezialfall eines ortsunabhangigen Durchlassigkeitsaffinors erhalt 
man 


O7h Ozh o*h o*h O7h Chae 
Kise a3 + kyy Oy b hss 0x2 + (Fix y + Key x) Ox dy + (ky + key) dy Oz T (hee + hes) Oz Ox Ue 
(16a) 


Wenn eine Richtungsanisotropie der Durchlissigkeit des Bodens vorliegt, so wird sie in prak- 
tischen Fallen meist durch einen symmetrischen Durchlassigkeitsaffinor § wiedergegeben 
werden kénnen. Das bedeutet, daB es drei aufeinander senkrecht stehende Richtungen gibt, in 
denen das Druckgefalle mit der Filtergeschwindigkeit richtungsgleich ist. Ist § auBerdem 
ortsunabhangig, so kann man eine generelle Hauptachsentransformation durchfihren. Wir 
legen in die drei Hauptachsenrichtungen ein neues Koordinatensystem (§, 77, ¢), das i. a. gegen 
die Achsen x, y, z verdreht sein wird. Im neuen Koordinatensystem erhalt der Durchlassigkeits- 
affinor die Form 


R& =| 0 Lee a0 Bs (17) 


Die Hauptachsenrichtungen werden sich in praktischen Fallen bereits aus dem Gefiige des 
Festkérpers erkennen lassen. 

Bei bekannten Hauptachsenrichtungen &, 7, € kann man nun zur experimentellen Ermitt- 
lung der Werte kz;, k,,, k-- einen in Richtung der Hauptachsen herausgeschnittenen Quader 
des Festkoérpers in drei aufeinanderfolgenden Versuchen in den drei Richtungen €, 7, ¢ durch 
strémen lassen. Man miBt die in der Zeiteinheit durch die Flacheneinheit der betreffenden Quer- 
schnitte durchstrémenden Flissigkeitsmengen (sie sind nach Definition gleich den Komponenten 
Ve, UV», Ve der Filtergeschwindigkeit )) und bestimmt die zugehérigen Standrohrspiegelgefalle 
J:, J,, J;. Es gelten dann die Beziehungen 


Vg = hee Ie, My = heyy Ins 1 = Bee Ses ey) 
aus denen sich die Werte ke;, k,,, und k-- bestimmen. Da der Durchlassigkeitsaffinor im neuen 
Koordinatensystem (&, 7, ¢) die Form (17) hat, geht (16a) tiber in 


OPh 07h Oh 
Kin On? r Kee oc? 0. (16b) 


55 Qe2 I 


Um die Gleichung (16b) auf die Form der Laplaceschen Differentialgleichung zu bringen, 
fiihren wir dimensionslose Werte fiir die Komponenten des Durchlissigkeitsaffinors ein durch 
die Beziehungen 


k é eres k® i lg ae eo 
312 ae V TIL EE 9 Ky y ay ) Kun Kin ? ke ar \ Kun ke, (19) 
mit Kyy = ke:k,, ke. Ky ist der dritte Skalar des Durchlassigkeitsaffinors. 


Durch eine reine Dehnung des Koordinatensystems (£, 7, €) vermittels der Transformations- 
gleichungen 


E=VRRE, 1 =Ving, C= Yi? (20) 
geht (16b) itiber in die Laplace-Gleichung 
: oth ect eas (166) 


=o I aA Tt aan 

gz | 07? | gee 
Im verzerrten Bereich verschwindet demnach die Anisotropie. In einer nachfolgenden Arbeit 
soll dies naher ausgefiihrt werden. 


5. Die Randbedingungen und Anfangswerte fiir die Standrohrspiegelhéhe. In praktischen 
Fallen kommen folgende Randbedingungen ! in Frage: 

1, Die Randflache ist undurchlassig. Sie stellt dann eine Stromlinienflache dar. 

2. Die Randflache bildet eine Grenzflache zwischen dem Strémungsbereich und der fest- 
kérperfreien Fliissigkeit mit statischer Druckverteilung. Hier wird auf die Randfliche die 
jeweilige hydrostatische Druckverteilung aufgepragt. 


* Vgl. H. F, Rossbach, Ing.-Arch, 12, (1941), S. 222; G. Hamel, Z. angew. Math, Mech, 14 (1934) S$. 131. 
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3. Die Randflaiche wird gebildet durch eine Fliche, langs der der Uberdruck in der Fliissig- 
keit tberall verschwindet (p = 0). Sie soll als Grundwasserspiegel bezeichnet werden. Bei 
instationaren Strémungen wird dieser mit der Zeit veranderlich sein. 

4. Die Randflache ist eine Sickerflache (Grenzflache zwischen Festkérper und freier 
Atmosphare). Es ist dann in ihr ebenfalls dauernd p = 0. 

Ks ist zweckmafig, bei der Berechnung des Strémungsverlaufes von der Standrohrspiegel- 
ee auszugehen und es sollen daher die vier Randbedingungen fiir die Funktion h formuliert 
werden: 

Randbedingung 1: Der Vektor der Filtergeschwindigkeit ) mu stets normal auf 
dem Vektor der Flachennormale der undurchlissigen Randflache stehen. Es sei 11 ein belie- 
biger Vektor senkrecht zur Randfliche, seine Komponenten seien n,, ny, n, baw. ne, n 
Ks gilt dann aes) ‘i 


"? Ne. 


Ta — 0) 21 
oder mit (12) SC 
He (eV hk) 0", (22) 
In Komponenten nach x, y, z lautet (22) 
toh Setae Gh yess dh ah ah ah 
Ny (hes Om kay dy on 5) Ny (Ass Ox ae yy ay =r ky. z) 
ai, ewes) ah 
a (tee 5, by 5 hes.) =0, (22a) 
in Komponenten nach &, 7, ¢ 
oh oh oh 
Ng kee FE + 1, ky» an + ng kee 3 = 0. (22b) 


Randbedingung 2: die Druckverteilung lings der Randflaiche ist die jeweilige 
hydrostatische p = p;, also bekannt. Aus (11) erhalt man dann die Werte von h langs der 
Randflache: 


Grundwassersplege/ 


| ee > + = =konst. (23) 


Randbedingung3 und 
4: langs des Grundwasserspiegels 
und langs einer Sickerflache ist 
p = 0 und daher nach (11) 


U 


a 
o 
oO 


(24) 


Die angegebenen Randbedingun- 
gen haben sowohl fiir denstationaren 
als auch fiir den instationaren Stré- 
mungszustand Giiltigkeit. Beim 
instationaéren Zustand tritt jedoch 
die Besonderheit auf, daB die Grund- 
wasserspiegelflache (Randbedin- 
gung 3) selbst eine unbekannte Funk- 
tion der Zeit ist. Damit ist auch der 
Bereich, in dem die Differential- Abb. 3. Zur Ableitung der Bewegung des Grundwasserspiegels. 
gleichung fiir h zu lésen ist, eine. 

Funktion der Zeit. Hingegen muB die Grundwasserspiegelflache zu Beginn der instationaren 
Bewegung gegeben sein (Anfangsbedingung). Im stationadren Fall hingegen ist die Form der 
Spiegelflache unbekannt und muf bei der Lésung mitbestimmt werden. 

Um fiir den instationaren Fall die Bewegung des Grundwasserspiegels rechnerisch zu ver- 
folgen, betrachten wir nach Abb. 3 ein Volumelement, das sich in vertikaler Richtung von 
einer undurchlassigen Flache bis zum Grundwasserspiegel erstreckt. Wir legen ein kartesisches 
Koordinatensystem (x, y, 2) zugrunde, wobei die z-Achse vertikal nach oben zeigt. Den Koor- 
dinatenursprung legen wir der Ubersicht halber so tief, daB die ganze undurchliassige Rand- 
flache im betrachteten Strémungsbercich iiber der («, y)-Ebene zu liegen kommt. 
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Die vorgegebene undurchlassige Randflache habe die Gleichung 2, = 2,(«, ¥), die gesuchte — 
Grundwasserspiegelflache sei z, = z(x,y ,t). Wir setzen beide Funktionen als stiickweise — 
stetig und mindestens einmal nach x und y stiickweise differentiierbar voraus. Das betrachtete _ 
Volumelement befinde sich an einer Stelle x, y, an der beide Flachen stetig und differentiierbar | 
seien. Es soll nun fiir das betrachtete Volumelement eine Bilanz der ein- und austretenden 
Flissigkeitsvolumina aufgestellt werden. Das in der x-Richtung wahrend dt eintretende 


Flissigkeitsvolumen betragt 


das austretende 


dx 


dy dt role + Foy t)de 


Analoge Ausdriicke ergeben sich fiir die y-Richtung. Zur Zeit t befindet sich in dem betrachte- _ 
ten Volumelement das Flissigkeitsvolumen 


z,(t) 
Vo = daay, | n(x, ¥, 2) dz, (25) 


worin n(x, y,2z) dasi. a. ortsabhangige Porenvolumen bedeutet. Wahrend dt vermehrt sich 


demnach V,, um 

Tea = in(x,.95-2 cid dy dt (25a) 
at eka NE Fa a 
Diese Vermehrung mu8 gleich sein dem Uberschu8 des in das Volumelement eintretenden 


Flissigkeitsvolumens gegeniiber dem austretenden. Es besteht somit die Beziehung 


% 2, dx . 02, dx 
Hex 2 TS 
‘ d 
dy s(x — F595 5 t)de — dy | ve(= + Fs 9 #1) de 

_ bx, de @z, dx 
Se eso} TO on 

aay ae 

DT Giy D a Tg Ds 

+ dx | ny(my pia i) de — dx i vy(% ¥ Tee t) ds 
pe ee 
2 ag, 2 ta age 
+¥.t 
= n(x, y, m4) MEY aa dy . 


Nach Ausfiihrung der obigen Integrationen und Division durch dx dy ergibt sich 
2,(t) 


Lean nad Aa(4%, y) Abies 
ie r 5) de + a Ys Ht) “eS t 94%, ys a5) 2 


Oz(x, y,t Oz(x, ¥, t 
—t(% V5 Zo, t) a — vy(X, ¥, %, t) Oaa(%, Y, t) 


02_(%, ¥, t) 
ae UCR. (26) 


SA EP 


Bei Beachtung der Kontinuitatsgleichung (13) ergibt sich fiir den ersten Term links 


%,(t) 3,(t) 
Ov, dvy dv, 
; | (3: ea =| es | Oz dz = v(x, Y> %29 t) a. v(x, Yo F1> t) : (27) 


Da die undurchlassige Randflache 2,(x, y) stets Stromlinienflache sein mu, gilt fiir sie die 
Gleichung (21). Als Normalvektor wahlen wir 


a7 Oz, Oz, 
n anf — Grey + es. 
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Damit liefert (21): 
) Oz 
a) ,( x, Yo 219 t) ae =F V4(%, Vo 2415 t) By +o v,( 2, Vo 245 t) (es (28) 


Dies ist demnach die Bedingung dafiir, daB die Flache z,(x, y) eine Stromlinienfliche bildet. 
Aus (26), (27) und (28) ergibt sich 


Oz9(x, ¥, t) 1 02,(x, y, t) 0z_(x, ¥, t) 
ae OES 0, Ko 2V'y Bast) — V4 iy Ys Zyy t)— a Uys Yo Za» 2) 2 rhe (29) 
Dies stellt eine Differentialgleichung fiir die Bewegung des Grundwasserspiegels dar. Da zur 
Zeit t = 0 die Lage des Grundwasserspiegels vorgegeben ist, so kann mittels (29) die neue Lage 
nach der Zeit dt bestimmt werden, wenn b(x, ¥, %, 0) bekannt ist. 


In vielen Fallen denkt man sich die Grundwasserstrémung zur Zeit t = 0 durch eine plétz- 
liche Stérung des Gleichgewichtszustandes ausgelist. Dann ist zur Zeit t = 0 streng genom- 
men noch iiberall ») = 0. Da zufolge (5) dann auch die Reibungskraft $h,, noch Null ist, konnen 
die auf den Fliissigkeitskérper wirkenden Krifte nicht im Gleichgewicht sein. Es werden daher 
Beschleunigungen auftreten, die zunachst nicht vernachlassighar sind. Doch wird sehr bald 
ein quasistationarer Strémungszustand erreicht, da der Festkérper laut Voraussetzung 6 der 
Flissigkeit schon bei kleinen Geschwindigkeiten einen geniigend grofen Reibungswiderstand 
entgegensetzt. Wir kénnen daher diese ,,Anlaufzeit“* im Vergleich zu dem ZeitmaB, in dem 
dann der quasistationare Strémungsvorgang ablauft, unberiicksichtigt lassen; es ist daher 
berechtigt, zur Zeit t = 0 bereits die Geschwindigkeitsverteilung der quasistationaren Str6- 
mung in die Rechnung einzufiihren. Um diese mit v(x, y, z, 0) bezeichnete Geschwindigkeits- 
verteilung zu erhalten lést man (16) bzw. (16c) mit den zugehérigen Randbedingungen. Ist 
h(x, y, 2,0) auf diese Weise ermittelt, so kann man aus (12) die Geschwindigkeitsverteilung 
D(x, y, 2,0) finden. Man kann dann schrittweise dasselbe Verfahren fir die aufeinanderfolgen- 
den Zeitelemente, die man praktisch durch endliche Zeitintervalle approximieren wird, wieder- 
holen und erhalt so den zeitlichen Verlauf der instationaren Strémung. 

Man erkennt, daB die Strémung stationar wird, sobald der Klammerausdruck auf der rech- 
ten Seite von (29) verschwindet. Dies bedeutet, analog zu (28), daB dann der Grundwasser- 
spiegel zu einer Stromlinienflache geworden ist. Dies ist also gleichzeitig die Randbedingung 
fiir den Grundwasserspiegel im stationaren Fall. 

Es ist also durch die Differentialgleichung (16) bzw. (16c) in Verbindung mit den Rand- und 
Anfangsbedingungen und den Gleichungen (11), (12) und (29) das Problem der instationdren 
Grundwasserstrémung vdéllig bestimmt. 

Man kann in der Beziehung (29) die Geschwindigkeitskomponenten vermittels (10) und (12) 
auch durch die GréSe h ausdriicken. Man erhalt dann 


ay _ [1 ah ah ah ah ah , , ah\ dz, 
Gees \n (ae ze Redias 5 he) (Fs Ox hey ay thes z) Ox 
ah ah ah\ az, 

eee fe pee ce (30 

( %* 0% ee oy r ") 5) alae, der Spiegelfliche ( ) 


Dies gilt, ebenso wie (29), fiir einen beliebigen, i. a. unsymmetrischen, ortsabhangigen Durch- 
lassigkeitsaffinor. Wenn der Durchlassigkeitsaffinor im betrachteten Strémungsbereich orts- 
unabhangig und symmetrisch ist, geht durch die generelle Hauptachsentransformation, die auf 
die Form (18) fiihrt, die Beziehung (30) itiber in 


Oz, 1 oh k Oh O25 dh 0z.\ 
FOE Ox = On OY 


lk, 
tt ala yom) | CO eee) 


langs der Spiegelflache 
Es ist auch méglich, in diesem Fall die Gleichung fiir die Bewegung des Grundwasserspiegels 
direkt im Hauptachsensystem (£, 7, ¢) abzuleiten. Es mu hierbei nur vorausgesetzt werden, 
daB® die Richtung einer der Achsen, z. B. die der ¢-Achse, die Grundwasserspiegelflache im 
betrachteten Strémungsbereich itiberall schneidet. Gleichung (30a) nimmt dann mit den 
analogen Bezeichnungen die Form an: 


aly bin tama dh ah ay ah ay 
= Kee ae — Agee ae — Fun 5 : 
ot n(&, 7, Ce) | a EE EL Os On On lings der Spiegelfliche 


(30b) 
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Ist der Festkérper homogen und isotrop, dann eriibrigt sich eine Koordinatentransformation 
und es spezialisiert sich (30) zu 

OZ k {0h oh Oz, Oh Az} (30c) 
at re | Oz Ox Ox dy OY tangs der Spiegelfliche 


Die bisherigen theoretischen Untersuchungen iiber die Grundwasserspiegelbewegung ver- 
wendeten stets den Dupuitschen Ansatz 1. Bei diesem geht man von der vereinfachenden An- 
nahme aus, daB das Standrohrspiegelgefalle fiir alle vertikal unter der betrachteten Stelle des 
Grundwasserspiegels liegenden Punkte gleich dem Spiegelgefalle an dieser Stelle sei. Diese 
Annahme trifft nur dann naherungsweise zu, wenn alle Stromlinien annahernd horizontal ver- 
laufen. AuBerdem geht beim Ubergang in die stationére Strémung nach dem Dupuitschen 
Ansatz die Spiegelflache nicht in eine Stromlinienflache tiber. Die hier entwickelten Beziehun- 
gen (30) bis (30c) hingegen sind von solchen Voraussetzungen iiber die Strémungsrichtung 
vollig frei und geben, wie schon erwahnt, den Ubergang zur stationadren Strémung richtig wieder. 

Grenzen zwei verschiedene Festkérper 1 und 2 lings ciner Flache aneinander, so wird der 
Durchlissigkeitsaffinor beim Uberschreiten dieser Grenzflache eine unstetige Anderung 
erfahren. Das hat zur Folge, daB die Stromlinien an der Unstetigkeitsflache gebrochen werden. 
Dieses ,,Brechungsgesetz soll nun ermittelt werden. Wir betrachten einen beliebigen Punkt 
der Trennflache und bezeichnen mit §, bzw. %, den Durchlassigkeitsaffinor im Material 1 baw. 
2 an dieser Stelle. Mit den aus Abb. 4 ersichtlichen Bezeichnungen ergeben sich gemaB (12) die 
Beziehungen 


0, REV Gs vs =— R cea Vad Fu (31) 


Aus Griinden der Kontinuitéat kann sich die Komponente der 
Filtergeschwindigkeit senkrecht zur Trennflache beim Ubergang 
von 1 in 2 nicht andern. Es gilt daher 


ae es Oa (32) 
Abb. 4. Die Brechung der Stromlinien , > 
an einer Unstetigkeitsflache (beider- (1, — Einheitsvektor senkrecht zur Trennflache). Da h, ebenso 
seits anisotroper Boden). f E 3 . ‘ 2 ¥ 
wie p, eine eindeutige Funktion des Ortes sein muBb, mub langs 
der Trennflache h, = h, sein. Es miissen daher auch die Projektionen von \/h, und \/A, auf 
die Tangentialebene der Trennflache im betreffenden Punkt iibereinstimmen. Demnach gilt 


Vay — 1 (Vain) = Vag —1 (Va * tn). (33) 
Um zu einer Formulierung des ,,Brechungsgesetzes“‘ zu kommen, setzen wir 
Y= Oe by (34) 


und bezeichnen § als ,,Brechungsaffinor’. 
Wegen der Linearitat des Darcyschen Gesetzes ist anzunehmen, daB B nur von §t,, %, und 
n, nicht aber von », abhangt. 
Aus (32) und (34) folgt 
- [(G — By dy] = 0". (35) 
Da (35) fiir beliebige Vektoren bp, Giiltigkeit haben muB, ist © — % ein planarer Affinor, der in 
der Form 
© —B=(,q) + ¢. 1) (36) 
geschrieben werden kann *. Hierin bedeuten § und t zwei beliebige feste aufeinander senkrechte 
Einheitsvektoren in der Tangentialebene senkrecht zu 11; q und t sind zwei noch zu bestimmende 
Vektoren. Aus (33) mit (31) und (34) ergibt sich nach einfachen Umformungen 


(Stes Bi RE) 5 Oy = ER BL Oe ath (37) 
Aus (37) erkennt man, daB K-!- B — K- ein is ee Affinor ist, der in der Form 
Sat che Oa (tt) (38) 


geschrieben werden kann. Hierin ist p ein noch zu bestimmender Vektor. Aus (36) und (38) 
mit © = (8, 8) + (t,t) + (1, n) folgt 


B=K- (1p) + 2=(8,8 —gt+(tt—)+(a,n), (39) 


. Kozeny, a. a. O. S. 433. 
Lagally, Vektorrechnung, S, 48, 2. Aufl. Leipzig 1934, 


ps 
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worin 


x == RX, 3 Cet: (40) 
bedeutet. 


Wenn man alle vorkommenden Affinoren im Dreibein 8, t, 1 darstellt (Indizes-1,:2, 3), 80 
kann man (39) dazu benutzen, die Komponenten von % zu ermitteln. Zu diesem Zweck 
schreibt man (39) in Matrizenform an und vergleicht gliedweise. Dann erhalt man aus der 
dritten Zeile die Komponenten von ). Damit ergeben sich die Komponenten q und t und 


schlieBlich der Affinor 8 in der Form 


k@ k() k2 
13 13 
a k(2) Is) lye a k(2) 32 13s a} (1 33) \ 
oe k(2) Kea) k(2) (41) 
Loy — 23 Loy ss Se le Ihe el ees 1 —_— I ‘ ; : 
\ k 2) oa k(2) 32 23 He§2) ( 33) 
ll 


Hierin bedeuten die ae bzw. 1;; die Komponenten von St, bzw. im Dreibein 8, t, 11. Damit ist 
das allgemeine Brechungsgesetz gefunden. Aus (31) und (34) kann man auch den Zusammen- 
hang zwischen \/h, und \/h, ermitteln. Es ergibt sich 


Viz = (81+ BK) > Vy. (42) 
Im Falle der Isotropie beider Festkérper gilt 
Rh, Ss - und Ro hy © « (43) 
Man erhalt dann durch Spezialisierung von (40) 
IB wae 
1 —_2(5 
ps hs, (40a) 
und nach (41) 
ky 
( = 0 0 
% =| 0 ue 0 |: (SU) sicetetighetafliehes (eidemneits. otzoper) Haden, 
ky ky > k,). 


0 
Legt man den Vektor ¢ in die durch », und n aufgespannte Ebene, so erhalt man aus (34) 


und (41a) 


lp a US ee LS etal Fo (44) 
ky 
Der Vektor v, liegt demnach ebenfalls in der durch », und nt aufgespannten Ebene (Einfalls- 
ebene). Sind a, und a, die Winkel von b, und b, gegen 11, so erhalt das Brechungsgesetz der 
Stromlinien fiir den Fall der Isotropie beider Festkérper die Form 


tg Oo, ky 
ae et (45) 
: 2. Oth Ab OP ae : Boge ae 
Bezeichnet man mit Pee ea die Komponenten von \/h im Dreibein 8, t, 11, so erhalt man 
nm 
analog die Beziehungen 
dhy , Oh, — ah, Ah, hr, hy 
as: (Beate ert Bis ee one >) 


was ebenfalls auf das Brechungsgesetz (45) fiihrt (Abb. 5). 

Kin Sonderfall liegt vor, wenn der durchstrémte Festkérper an einen festkérperfreien 
Flissigkeitsbereich mit hydrostatischer Druckverteilung grenzt (Randbedingungen 2). Diesen 
Fall kénnte man durch Spezialisierung von (42) gewinnen; es ist jedoch einfacher, ihn direkt 
zu untersuchen. Da bei hydrostatischer Druckverteilung die Standrohrspiegelhohe in dem 
betreffenden Flissigkeitsbereich konstant ist (h, = konst., \/h, = 0), ist wegen der Stetigkeit 
yon h die Projektion von \/h, auf die Tangentialebene der Grenzflache ttberall Null. Es ist 
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daher \/h, itberall senkrecht auf die Grenzflache gerichtet. Die Richtung der Stromlinie | 


ergibt sich dann aus der Beziehung (31): 
be = —R,° V hy. 
Man erkennt, da® im Fall eines isotropen Festkérpers auch die Stromlinien senkrecht auf die 
Grenzflaiche stehen. 
6. Die Korn-zu-Korn-Kriifte. Um einen Zusammenhang zwischen dem Spannungsaffinor 


© x der Korn-zu-Korn-Krafte und der Geschwindigkeitsverteilung bzw. der Standrohrspiegel- 
hohe herzustellen, gehen wir von (3) aus. Bei Verwendung von (11) erhalt man die Beziehung 


ace eee 
We : ©x = (WK ve a Ve at Vw Vh 3 (47) 

Driickt man in (47) die GréBe \/h nach (12) durch » aus, dann erhalt man 
Vs Sx = et et ye Gee oe (48) 


\/ : Gx bedeutet die auf die Volumeinheit des Gemisches bezogene Resultante der auf den 
Festkérper wirkenden Korn-zu-Korn-Krafte. Man erkennt dies sofort aus der Definitions- 
gleichung: 

Je =m $8888 

V0 

Aus (48) erkennt man den Zusammenhang zwischen \/ - ©x und der Filtergeschwindigkeit b. 
Fir » = 0 wird nur dem um den Auftrieb der Kérner verminderten Gewicht des Festkérpers 
das Gleichgewicht gehalten. Da durch (47) oder (48) nur die Divergenz des Spannungsaffinors 
bestimmt ist, laBt sich daraus ©x nicht eindeutig bestimmen. Es liegt hier ein statisch un- 
bestimmtes Problem vor, das erst gelést werden kann, wenn ein Zusammenhang mit den 
Deformationen des Festkérpers hergestellt wird. Da dieser Zusammenhang wesentlich von 
der physikalischen Beschaffenheit des Festkérpers bzw. des Gemisches abhangt, wird er fiir 
jede Bodenart eigens ermittelt werden missen. Wir wollen hier auf dieses Problem, das zur 
Beurteilung der Rutschungsgefahr von groBer Bedeutung ist, nicht naher eingehen. Es ware 
eine Aufgabe der Bodenmechanik, hieriiber allgemeinere Untersuchungen anzustellen. 

Eine Integration der hier fiir die Grundwasserbewegung abgeleiteten Differentialgleichun- 
gen wird auf rein analytischem Wege nur in einfachen Fallen mit tragbarem Rechenaufwand 
méglich sein. Hingegen bietet die Relaxationsmethode die Méglichkeit, auch komplizierte 
Falle (z. B. mit graphisch vorgegebenen, analytisch schwer zu formulierenden Randern und 
instationarer Strémung) zu behandeln. 

In einer nachfolgenden Arbeit soll auf diese Lésungsmethode an Hand von praktischen 
Beispielen naher eingegangen werden. 


7. Zusammenfassung. Nach einer eingehenden Analyse der wirkenden Krafte wurden die 
Gleichungen fir die Grundwasserstrémung und die Korn-zu-Korn-Krafte in verallgemeinerter 
Form abgeleitet. AnschlieBend wurden die miglichen Randbedingungen formuliert und die 
Bewegung des Grundwasserspiegels bei instationdren Strémungszustanden rechnerisch verfolgt. 
In einer nachfolgenden Arbeit soll die Anwendung der entwickelten Gleichungen fiir die Lésung 
praktischer Aufgaben gezeigt werden. 


(Eingegangen am 11, Marz 1954), 


Anschrift der Verfasser: Prof. Dr. G. Heinrich und Dr. K. Desoyer, Wien IV., Karlsplatz 13, 
Technische Hochschule, 
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Uber die Dimensionierung diinner Kreisplatten unter exzentrisch 
aufgebrachten Einzellasten’*. 


Von Kh, Nasitta. 


1. Einleitung. Es soll das Problem der durch eine Einzellast exzentrisch belasteten und am 
Rand frei aufgelagerten Kreisplatte behandelt werden. Das Ziel der Arbeit ist, die von 
E. Reifner 1 1935 abgeleitete Lésung dieses Problems durch cine einfache Naherungslésung zu 
ersetzen von der aus die schnelle und hinreichend genaue Dimensionierung einer exzentrisch 
durch Einzellast belasteten und am Rande frei aufgelagerten Kreisplatte méglich ist. 


2. Die Parameter in der Platienekene. Betrachtet man diejenige linear gebrochene kom- 
plexe Funktion, die das Innere des Einheitskreises einer w-Ebene auf die rechte Halbebene 
einer z-Ebene abbildet, so werden die in der w-Ebene konzentrischen Kreise r = konst. und 
0 <r<1 in eine um x = b exzentrische Kreisschar in die z-Ebene iiberfiihrt (Abb. 1). Die 


ty 


Z 


Abb. 1. Konforme Abbildung. 


Linien g = konst. bilden ebenso wie in der w-Ebene auch in der z-Ebene die Orthogonal- 
trajektorien zur ersten Schar, weil jede linear gebrochene Abbildung konform ist. Bringt man 
die Punkte J, IJ, III aus der ersten Ebene nach I’, II’, III’ = oo der zweiten Ebene, so ist 
damit die Abbildung eindeutig gegeben, und man erhalt 


Te OTs ee ert (1) 
Mit mw =re'? undz = x +7y entsteht hieraus die Parameterdaistellung 
x=bi—, yao? (2) 
wobei 
B=1+7r? —2rcosg 
sein soll. 


Weiterhin werden wir mit der reellen Darstellung (2) in einer reellen (x, y)-Ebene arbeiten, 
die zugleich die Plattenebene sein soll; der Parameter r = ry < 1 soll ein fiir allemal den 
Plattenrand beschreiben, wahrend r — 0 der ausgezeichnete Punkt des beniitzten Systems ist, 
in dem die singularen Beanspruchungen der Platte aufgebracht werden. Um einfachere Ver- 
standlichkeit zu erméglichen, mégen die hier benutzten Parameter als ,,exzentrosymmetrische 
Koordinaten“ bezeichnet sein. 


* Gekiirzte Fassung einer von der Technischen Universitat Berlin-Charlottenburg genehmigten 
Dissertation, Der Verfasser ist den beiden Berichtern, Herrn Professor Dr.-Ing. I. Szabé und Herrn 
Professor Dr, phil. E. Mohr sehr zu Dank verpflichtet. te 

ee E. Reipner, Math. Ann. 111 (1935) S. 777; diese Loésung des Problems ist in Abschnitt 10 durch 


Forme! (38) wiederge geben. 


86 Nasitta: Uber die Dimensionierung diinner Kreisplatten usw. Ingenieur-Archiv 


Damit im nachsten Abschnitt die Querkrafte und Schnittmomente der Platte fiir das vor- 
liegende spezielle System berechnet werden kénnen, muf jetzt der gesamte Parameterraum, in 
den die Platte eingebettet wird, betrachtet werden. Die Mittelflache der Platte liegt vor der 
Beanspruchung in der (x, y)-Ebene, die durch exzentrosymmetrische Koordinaten beschrieben 
wird (Abb. 2). Nach Aufbringen der Last entsteht in w-Richtung eine Durchbiegung w(r, @). 
Danach ist die Durchbiegungsflache durch folgenden Ortsvektor gegeben: 


u(r, p) = {a(rs 9) x(t G) » W(r, P)} > (3) 
wo x und y durch (2) bestimmt sind. Fihrt 
man die Abkiirzungen 


(8) +B) 


Ox\2 dy \? 
umd Ve (55) fis (55) =Ur (4) 


ein, so sind die Koeffizienten der ersten Grund- 
form der krummen Flache 1 


i) B= (t,t) O 4; BS (0g) = ee 


Abb. 2. Lage der Platte und der Einzellast im (x, y, w)-Raum. (Cs (to 3) == |/ —- We : (5a) 


3. Schnittmomente und Querkriafte der diinnen Platte fiir den Fall exzentrosymmetrischer 
Koordinaten. Fir kartesische Koordinaten werden die SchnittgréBen in Abhangigkeit von der 
Durchbiegung durch folgende Formeln gegeben *: 


0Aw 
My = — K (we + Wyy)>~ Qe) = —K—a_, | (6) 
0Aw 
My) = —K (Wey, + v Dy) > Qy) = —K Oy” 2 | ; 
M xy) = 1K (1 = v) Wey 9 
3 
wobei y die Querkontraktionszahl, E den Elastizitatsmodul und K = pies die Platten- 


steifigkeit bedeuten. 

Da bei einer diinnen Platte samtliche SchnittgréBen langs orthogonaler Parameterlinien 
allein durch die geometrische Gestalt dieser auf der Durchbiegungsflache verlaufenden Kurven 
geliefert werden, liegt es nahe zu versuchen, die vom Parametersystem abhangigen Schnitt- 
momente und Querkrafte durch absolute Invarianten der Parameterlinien darzustellen. 
Beachtet man, da die Theorie der diinnen Platte nur fiir soleche Durchbiegungsflachen gilt, 
deren Tangentialebenen kleine Winkel mit der (x, y)-Ebene einschlieBen, so lassen sich in der 
Tat die unter (6) gegebenen Ausdriicke im eben erwahnten Sinne invariant formulieren. 

Nach der soeben ausgesprochenen Voraussetzung fiir die Durchbiegungsflache hat offenbar 
die Normale n der Durchbiegungsflache dem Betrage nach kleine Komponenten in der (x, y)- 
Ebene. Damit ist das Quadrat des Betrages von 1 im wesentlichen durch nj) gegeben, so daB 
fiir den Fall der diinnen Platte nj), njy gegeniiber nj,,) vernachlassighar klein sind, also 


| = nly + nly + thy ~ they (7) 
gilt. 

Falls die Flache der Durchbiegung einer diinnen Platte mit beliebigen orthogonalen Para- 

metern u,v durch z(u,v) = {x(u, v), y(u,v), w(u, v)} dargestellt wird, so errechnet man ® 


aus der Bedingung (7) 
(n+ Yu) + wa ~ (a + yu)» 
(x8 + y2) + wh = (ad + 92) | ® 
Wenn E, G die Koeffizienten der ersten Grundform der Durchbiegungsflache sind, L, M, N 
diejemgen der zweiten Grundform ‘ so gestatten die unter (6) gegebenen Krafte und Momente 


1 Vel. B. Baule, Die Mathematik des Naturforschers und Ingenieurs, Bd. VII. Differentialzeo- 
metrie, S. 34ff., 2, Aufl. Leipzig 1945. : 

2 Vel. C. B. Biezeno u. R. Grammel, Technische Dynamik, Bd. I, S. 486 ff, 2. Aufl. Berlin-Géttingen- 
Heidelberg 1953. 

3 Vel. B. Baule, a. a. O. S. 39 u. 52. 

= MWralle tits lynn ene, O) Cy ae) aig 52% 
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bei Beriicksichtigung der in (8) formulierten Vernachliassigungen folgende invariante Dar- 
stellung, wobei « durch uw und y durch v ersetzt sind: 


8, N 0A 
Mw) = — Kg +»), se Sr rae | 
N hs 0A 
Mwy —=— K ( + yp z) Gl Qi) eS K ae Ss (9) 
; M 
M. 0) KG (da B) a 
(uv) ( ) VEG | 


Hierin bedeutet ds, das Linienelement der Linien v = konst., also ds, — VE du ; entsprechen- 
des gilt fiir ds,, also ds, = VG dv. Durch die Formeln (9) ist die Theorie der diinnen Platte 
durch Differentialinvarianten belicbiger orthogonaler Parameterlinien formuliert. Fiir exzen- 
trosymmetrische Parameter erkennt man nun, da die Beziehungen (5a) auf Grund von (8) 
ubergehen in 


BUs) 6 = V, WV=hC—-F =U. (5) 


Im Falle exzentrosymmetrischer Koordinaten ergibt sich fiir die Koeffizienten der zweiten 
Grundform der durch (3) gegebenen Durchbiegungsflache 


pee) 1 r 0U 1 oU 
pe ee 5 
W luv ( eee en aap *, 
(ts Fa) ey r oU lpels 
M i. = arr. (* Die ae dp Or or ) Nae 
ee) a pe nls 
Nee eperoue! | 
W per (to to Tiel See Oe Ca 


Um auch die Querkrafte nach (9) im vorliegenden System zu kennen, verbleibt noch die Auf- 
gabe, den invarianten Differentialoperator 4 in exzentrosymmetrischen Koordinaten anzu- 
geben. Die durch (1) gelieferte analytische Ab- y 
bildung laBt als solche Potentiallinien invariant ; 
d.h., wenn in der w-Ebene (Abb. 1) fiir eine 
Funktion AF = 0 gilt, so behalt diese Beziehung 
auch in der z-Ebene ihre Giiltigkeit, woraus 
folgt, da sich der A-Operator in exzentro- 
symmetrischen Koordinaten von dem in Polar- 
koordinaten nur durch einen Faktor unter- 
scheiden kann. Da sich der MaBtensor der 
Polarkoordinaten von dem der exzentrosymme- 
trischen Koordinaten nur durch den in (4) 
definierten Faktor U unterscheidet, gilt 1: 
L fet. tym. y- L.oFw 
G Toh. 3 Fe an) 


Damit sind durch (9), wenn man hierin u=r 
undv = @ setzt, unter Zuhilfenahme von (10) 
und (11) sowie mit (5) alle Querkrafte und Schnittmomente in exzentrosymmetrischen Koor- 


dinaten gegeben (Abb. 3). 


Abb. 3. Volumenelement der Platte mit Schnittkriften. 


4, Die Lésungen der homogenen Plattengleichung AAw = 0 fiir exzentrosymmetrische 
Koordinaten. Die Durchbiegung w einer durch Einzellast belasteten diinnen Platte muf 
bekanntlich der homogenen Differentialgleichung JAw = 0 geniigen. Wie man wei 2, lassen 
sich alle Lésungen der Bipotentialgleichung AAF = 0 fiir den Fall zwei freier Rander des 
Problems in Polarkoordinaten 0, w darstellen durch 


F(g, y) = % 0 file, ) + & @ cos p- fale, yp) + ¢s *fs(Q> ¥) - (12a) 
1 Vel. M. Lagally, Vorlesungen iiber Vektorrechnung, S, 325, Leipzig 1945. 
2 Vel. C. B. Biezeno u. R. Grammel, a. a. O. 5. 134 ff. 
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Hierin bedeuten die f;(i = 1, 2, 3) Potentialfunktionen — es gilt also Af; = 0 —, und die ¢; 
sind beliebige konstante Faktoren. Da es sich bei JA um einen invarianten Differentialopera- 
tor handelt, muB der eben ausgesprochene Satz giiltig bleiben, wenn einerseits AA und anderer- 
seits F(0, y) auf exzentrosymmetrische Koordinaten umgeschrieben werden. Wir denken uns 
also im folgenden AA fiir exzentrosymmetrische Koordinaten formuliert und befassen uns 
damit, F(9, y) auf G(r, ~) umzuschreiben. Bei Verwendung der Formeln (2) erhalt man aus 0? 
und @ cos p 


xp yt =" (1 +72 + 2rcosq) und «=o (1—r), 


Wenn man die Funktionen f; in neue Potentialfunktionen 8;(r, p) ibergehen lat, so entsteht 


nun aus (12a) 
1 > a A ie . 
G(r, ~) = 7 [r? (c, b? g, — cy bg, + ¢3 gs) + r cosy (2c, b? g, — 2 Ca Ba) 
+ (¢, 5? ra + Cy bg, + C3 83) | . (12b) 


Fat man die Klammern in diesem Ausdruck zu neuen Potentialfunktionen C; g; zusammen, so 
ist gezeigt: Alle Lésungen von AAG(r, y) = 0 in exzentrosymmetrischen Koordinaten lassen 
sich mit beliebigen C; und Ag; = 0 (i = 1, 2, 3) in der Form darstellen 


1 
G(r, p) ==zlGi gilt, y) + Cyr cos p: galt, ~) + Cy Ba(r, @)] - (12) 


Mit Hilfe dieser grundlegenden Aussage, iiber die im exzentrosymmetrischen System zur Ver- 
fiigung stehenden Lésungen der Bipotentialgleichung, ist nun eine vollkommene Ubersicht 
iiber die einfachsten Lésungsanteile gegeben. 

Zur Konstruktion dieser Lésungen ist festzustellen, daB wegen (11) im Falle exzentrosym- 
metrischer Koordinaten dieselben Potentialfunktionen wie im Falle der Polarkoordinaten 
verwendet werden kénnen. Geht man nun von dem Gesichtspunkt aus, nur die einfachsten 
elementaren Potentialfunktionen heranzuziehen, so wird man folgende Lésungen von Ag; =0 
ansetzen kénnen: 


see ae (13) 
fio Inte (14) 
£. = TF cOs.Q (15) 
1 
& = — cosy. (16) 
Durch Kombination der Funktionen (13) bis (16) mit den Koeffizienten ve - r cos @ und = 


aus (12) erhalt man zehn voneinander verschiedene Partikularlésungen der homogenen Platten- 
gleichung. Mit diesen bekommt man mit beliebigen c;, wenn man beriicksichtigt, daB nach den 
Grundlagen der Theorie w und die Neigung der Tangentialebene fiir alle den Plattenbereich 
beschreibenden Parameter beschrankt bleiben miissen und weiter w am Rande verschwinden 
soll, 


i |p g 2 2 
w(r,Q) = ae r tn a ++ €g (rg — 7?) + eg r (78 — r®) cosy + ce, (ré — 1?) cos? gy]. (17) 


Mit (17) ist die einfachste Lésung der homogenen Plattengleichung hergestellt, die geeignet ist, 
die Durchbiegung einer am Rande aufgelagerten diinnen Platte zu beschreiben. 


5. Singularitét der Lésung. Die Lésung (17) der homogenen Bipotentialgleichung soll nun 
henutzt werden, um fiir das vorliegende Problem eine einfachste Naherungslésung zu kon- 
struieren, derart, da fiir kleiner werdende Exzentrizitat ¢ (Abb. 1) die Naherungslésung immer 
besser wird und schlieBlich fiir den konzentrischen Fall ¢ = 0 mit der bekannten, symmetri- 
schen Lésung zusammenfallt; diese Tendenz der Naherung ist die fiir die Bemessung erwiinschte. 

Fiir r = 0 erhalt man aus der Grundlésung (17) 


w(0, ~) = 16 (cg + c4 cos? y) ; (18) 
d. h. es liegtinr = 0 eine unendliche Vieldeutigkeit vor. Hieraus erhellt, daB, wenn e, + 0 (17) 
nicht ohne weiteres zur Beschreibung einer Durchbiegungsflache benutzt werden kann. Es 


liegt jedoch nahe, die Stelle der wesentlichen Singularitat in r = 0 vorerst aus der Platte heraus- 
geschnitten zu denken, indem ein beliebig kleiner Kreis mit von Null verschiedenem Halb- 
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messer r, als innerer Rand der Platte eingefithrt wird. Wir werden nun die Grundlésung (17) 
den duBeren Randbedingungen der freien Auflagerung 1 am Rande r = ry anpassen, und dann 
wird es durch Interpretation der am inneren Rande entstehenden Kraft- und Momentenver- 
teilungen und anschlieBender Superposition méglich sein, die stérenden Krafte und Momente 
am inneren Rande der Platte in erster Naherung zu beseitigen und die fiir den Lastfall der Ein- 
zellast geforderte Beanspruchung allein zuriickzubehalten. 


6. Aufere Randbedingungen. Die im Falle einer freien Auflagerung zu befriedigenden 
Randbedingungen lauten {w}, =” = 0) und {Muy}, =r, =0. Wegen der ersten Auflager- 
bedingung ist w, = w,, = 0 am duferen Rande, und damit wird nach (9), (10) und (5) die 
zweite der Auflagerbedingungen 


file! 
{My}, =To = K) U (20,, 


Da fiir die Lésung (17) {w},_,, = 0 gilt, mu® diese nur noch der Randbedingung (19) an- 
gepaBt werden. Setzt man (17) in (19) ein, so entsteht 


LS a La 
Opi on es Cy Ors ter, 


=0. (19) 


Ary c, —Arpc, cosp — 419 (1 — 13) cp — 816 c3 cos p + 8 13 cy cos? — 4 rp €y(1 — 72) cos? p 

— (1 — 76) (1 — ») (ro ey — 2 9 ep — 2 reg cos — 2 1, ce, cos? gy) = 0. (20a) 
Da (20a) fiir aJle m gelten soll, miissen die Koeffizienten von 1, cos g und cos? p identisch ver- 
schwinden. Daraus flieBen drei Gleichungen zur Bestimmung der c,;; mit 


a = 4—(1 —y»v)(1—*7), Bp = (1 — 7) (1 +7) (21) 
erhalt man 
OC, 2,0 c= 10); | 
TG ioe OiCg see.) 5 (20b) 
Aréc, —Bey=0. | 


Eine Konstante ist noch frei verfiighar; wahlt man c,, als vorerst noch beliebig, so folgt aus 
(20b) 
a 2 8 r? 

nedie ate CS ie Maurer ee (20) 
Mit den unter (20) gegebenen Konstanten c; wird fernerhin weitergerechnet, und ohne beson- 
deren Zusatz sind stets die c; der freien Auflagerung gemeint. Damit ist den Auflagerbedin- 
gungen des hier behandelten Problems exakt geniigt, und wir wenden uns nun den Lastverhalt- 
nissen zu. 


7. Krafte und Momente am inneren Rand. Voraussetzungsgemaf ist die Platte frei von 
Flachenlasten, aber es befinden sich Lasten in Form von Querkraften und Momenten am 
inneren Rand r, der Platte. Um die Deutung und Vergleichbarkeit dieser Querkrafte und 
Momente zu erleichtern, leiten wir aus den nach (9) gegebenen SchnittgréBen der Dimension 
nach wahre Krafte und Momente folgendermafen ab: 


dM(,) = Mo) ds, = Mr) dy y 

IMG g) = M¢@ g) ds, = mq) dp , (22) 

dK,,) = Qo) ds, = ki) dy . 
Hieraus erhalt man unter Beachtung von (4) und (5) 

256 20 7 —26rT 
my == May, ~™,) =z Mens key =z Om: 

Mit Benutzung der Plattenlésung (17) errechnet man mit (9), (10), (11) und (5) die folgenden 
Momente und Kriafte, wobei sofort die Glieder, die mit r,—> 0 nach Null gehen, also die Poten- 
zen erster und héherer Ordnung von r, wegen der beliebigen Kleinheit von r, weggelassen 


4 = Es sei darauf hingewiesen, da fiir den einfachen Fall der starren Kinspannung c, cS 0 wird, wodurch 
die Singularitat verschwindet und man eine auferst einfache und handliche Lésung erhalt, die in bipolaren 
Koordinaten schon 1926 von E. Melan angegeben worden ist; fiir die hier benutzten Koordinaten wird 

1 PKA I 


i° 
Weingespannt = ee TPAD a + (r2 — a) : 


0 
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werden: 

2 2 

Toes K [470.9 £08 Th ae K2 c,7g (1 + ») cosg sin? , 
b Tg b 

5 : 
a c,7o (1— v) sin 2 p 
Mir )— K 2b Te 2 (24 a) 
K cos 29 9 COS /~ 
kr) aa b2 (e. ai C4 "9 ree €4T9 Te : 


\ 


Fiir eine Abschatzung der Wirkungen ist es nun zweckmabig, Krafte und Momente ver- 
schiedener GréSenordnungen beziiglich r, zu unterscheiden. Dann erhalt man 


( ( 
my = mi + mi}, Meg) =m), oy = ey + hG) + kG 


mit 
wv (al : 
7 hee +”) cos sin’ gp, j= —Ko, 
(24) 
2 2 
(2) C47 Y Cos2Q (te) C479 COS P 
2G Rae ea emia y By eae 
2 . 2 
(2) C479 (1— ») sin 2 (Qos C479 Cos 2 p 
Tg) oa K 2b roe kG) = K Bb? Ore? 


Man erkennt nun leicht folgendes: Die konstanten Krafte ke} am inneren Rand r, bilden 


zusammen eine resultierende Kraft, deren GréBe und Richtungssinn durch die noch freie Kon- 
stante c, festgelegt werden kann. Diese Wirkung allein ist also die gewiinschte Wirkung einer 
Einzellast; denn es ist von vornherein klar und 1a4Bt sich auch im konzentrischen Fall leicht 
nachweisen, daS ein Ring konstanter Querkrafte an einem beliebig kleinen Loch angreifend 
dieselbe Wirkung hat, wie eine Einzellast, die anstelle der Querkrafte auf einer ungelochten 
Platte in der Mitte des vorher vorhandenen beliebig kleinen Loches angreift und die der GréBe 
nach gleich der Summe der Querkrafte ist. Die auBer ke} angreifenden Krafte und Momente 
bezeichnen wir ein fiir allemal als Stérsysteme, und zwar wegen ihrer verschiedenen 
GréBenordnungen die mit dem oberen Index (1) als Nebenstérsystem und jene mit oberen 
Index (2) als Hauptstérsystem. Wie aus (24) ersichtlich ist, sind die beiden Stérsysteme 
,,Gleichgewichtssysteme“ 1, Da nach Festlegung von ki} keine Konstante in der Lésung (17) 
mehr verfiigbar ist, kénnen die Stérungen nur noch durch Uberlagerung zu Null gemacht 
werden. Das Stérsystem (1) kann wegen r, beliebig klein zu einem endlichen singularen Moment 
inr =0 zusammengefaBht werden. Dazu macht man sich klar, da ein singulares Moment, wenn 
iiberhaupt, nur um die Achse g = + baw. a z wirken kann, da die durch (17) dargestellte 
Durchbiegungsflache symmetrisch um die Achse p = 0 bazw.a verlauft. Die Krafte kG) liefern 


Momente dM beziiglich der Achse py = > bzw. = z und man erhalt, wenn s,, der wahre Halb- 
messer des beliebig kleinen Loches um r = 0 mit dem Parameterwert r, ist, 
dM = ki} s, cosp dy. (25a) 


Beachtet man nun, daf s,. = i| Ob. = f ju dr ist und wegen (4) und (2) yu in der Umgebung 
0 a 

von Null naherungsweise 2 6 ist, so erhalt man aus (25a) durch Integration von p = 0 bis 2 a 

fiir das resultierende Moment des ,,Gleichgewichtssystems“ von Kraften ky, 


Fi ESCA fas 


MS (25) 


In gleicher Weise wirken die Momente m/;} aus (24), deren Resultierende in Richtung y = 0 


bzw. a Null ist, und in Richtung g = > bzw. zm ergibt sich durch Integration 


1 In der vorliegenden Arbeit sind yon der iiblichen Bedeutung des Wortes », Gleichgewichtss ystem‘ 
abweichend solche Systems von Kraften oder Momenten gemeint, deren Resultierende bzw. deren 
Integralmittelwert iiber m Nullist. 
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Mithin ist das singulare Moment in r = 0, als das sich das Nebenstérsystem S() = {ki} » my} 
in seiner Wirkung deuten laBt, gegeben durch M() = M + M, und man erhilt 
Ker a 
M0) =—2 2% (3 —»). (27) 


Die GréSe der Nebenstérung kann nun bei jedem Zahlenbeispiel nach (27) ermittelt werden. 
Man kann jedoch sagen, daB die Wirkung des Stérsystems S(!) sehr klein gegeniiber der- 
jenigen von S@) — {ki} » mi} mira) sein wird. 
Stellt man die Krafte kG} und Momente m}}}, mir’ y) schematisch nur an ihren Maximalstellen 
4 


J 3 ; 
dar, d.h. alsog =0, + » Ty > % bzw. 9 = Fg Gm Ait, so ergibt sich endgiltig fiir die 


Gesamtheit der am inneren Rand 
der Platte auftretenden Kriafte und 
Momente das in Abb. 4 dargestellte 


Bild. -¢-F 
; A FL) ap 
8. Wirkung desHauptstérsystems. APY) 
Da S( als endliches Moment ge- 
deutet werden konnte, erkennt man, os frnt3}(p-0) 
dafS nur das Vorhandensein von S(2) ? fees 
die unendliche Vieldeutigkeit der rs) 
Lésung (17) in r =0 bewirken kann, eS, 
da ein endliches Moment eine solche ¢ ; 
Wirkung auf eine Platte bekanntlich min (p= $0, 4 
nicht hat. 7 bea 
Um das Hauptstérsystem S(?) = pH 
eliminieren zu kénnen, kommt es Abb. 4. Krafte und Momente am inneren Rand der Platte. 


nun im weiteren darauf an, allein 

seine Wirkung abzuschatzen. Zu diesem Zweck soll folgende einfache Durchbiegungsflache, 
die ein ahnliches ,,Stérsystem® und dieselbe Singularitaét enthalt, zum Vergleich herangezogen 
werden: 


w= ao*In- + dy (0? ~ eo] cos2y, (28) 
: : : : IB 2 Sey 5 
wobei 0, y zentrische PoJarkoordinaten sind und a = — ta K und dy = See sollen. 


Mit diesen Konstanten a und dp, die weiterhin die hier definierte Bedeutung beibehalten sollen, 
erfiillt w aus (28) die Randbedingungen einer freien Plattenauflagerung, wenn 0 = Q, den Platten- 
rand beschreibt. AufBerdem hat w in 9 = 0 die gewinschte algebraische Singularitat, und da 
der Faktor vor cos 2 p in (28) genau die Lésung einer frei aufgelagerten und im Mittelpunkt 
durch die Einzellast — P belasteten Kreisplatte ist, liefert (28) fiir eine Schnittlinie y = konst. 
eine solche Durchbiegungskurve, wie sie auf der Durchbiegungsflaiche einer Platte fiir alle w 
entstehen wiirde, die allein durch die Einzellast — P cos 2 y = Ps belastet und am Rande frei 
aufgelagert ist. Insbesondere haben wir auf der Linie y = 0 bzw. eine solche Durchbiegungs- 
linie, als ob eine scheinbare Einzellast Ps von der GréBe — P in 9 = 0 auf die Platte wirken 
wurde, — ,,scheinbare Einzellast in 9 = 0° deshalb, weil wir ja wegen der Singularitat von (28) 
wiederum wie bei (17) den Punkt 9 = 0 ausschlieBen und uns einen beliebig kleinen Kreis mit 
_dem Halbmesser 9, >0 um 0 = 0 geschlagen denken miissen. Wie bei (17) werden wir nun die 
am inneren Rand 0, der Platte wirklich auftretenden Krafte und Momente berechnen. Es 
ergibt sich analog zu (24a) 


Mo) = K—*—“cos 2y, 
2 
My») = —K eee sin 2 Y> (29) 
2 
a ky = — K2%% cos2y. 


AuBerdem liefert (28) nicht wie das bei (17) der Fall war, AAw = 0, sondern man erhalt mit 
AAw = p/K fir die flachenhafte Belastung, durch die die Durchbiegungsflache nach (28) 


7s 
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beansprucht wird, 


p= K Tf cos 2 y. (30) 


Wir wollen zunachst die Belastung (30) auSer acht lassen und zuerst das System (29), das 
wir durch S, kennzeichnen, mit dem Hauptstérsystem S() aus (24) vergleichen. Um fir beide 
Systeme die Lage gegeniither 9 baw. r = 0 und die Gré®e der Krafte und Momente vergleichbar 
zu machen, gehen wir in S() aus (24) auf wahre Radien und Winkel itber; wie schon bei (25a) 
gezeigt, erhalt man fiir den Halbmesser des beliebig kleinen Loches in der Umgebung von Null 
s,. =2br,. Wir setzen nun ohne Einschrankung der Allgemeinheit voraus, daB s,, = @, ist. 
Beachtet man weiter, daB in der Umgebung von r = 0 die Winkel py der exzentrosymmetri- 
schen Koordinaten wahre Winkel sind, also p = y, so errechnet man unter Zuhilfenahme von 
(29) und (24) leicht 


2 
(2) — | Ca To 
Mir) == F d, A MM (0) ° 
rail) ces | a8) on (31a) 
(r@) 2d, 02 (oy) ? 


Hieraus erkennt man, da8 das Hauptstérsystem S() aus (17) dem System S, aus (28) ahnlich 
ist, da sich Krafte und Momente nur durch einen konstanten gemeinsamen Faktor unter- 
scheiden, also 


SQ) Sa (31) 


Wir wenden uns nun der Wirkung der flachenhaften Belastung p in (30) zu. Diese Belastung 
ist tiber die ganze Platte als,,Gleichgewichtssystem“ verteilt derart, daB sie im ersten Oktanten 
positive Werte, in den nachsten beiden negative 
Werte annimmt und in den zwei darauffolgenden 
Oktanten wiederum positive Werte zeigt usw., wie 
man leicht aus (30) unter Beriicksichtigung, daB a 
negativ ist, abliest (Abb. 5). Um die Wirkung von 
p beziiglich seiner Verteilung tiber y abzuschatzen, 
machen wir uns klar, da wegen des viermaligen 
Vorzeichenwechsels der Belastung p die Wirkung nur 
ein Bruchteil derjenigen ist, wenn p einsinnig bezig- 
lich y tiber die ganze Flache verteilt ware. Man kann 
eine obere Schranke fiir diesen Bruchteil der Wirkung 


3 Sak 3 
angeben, wenn man sich paralfel zur Linie y = raid 


Ge ; 
bzw. 7 einen Balken aus der Platte herausgeschnitten 


Abb. 5. Belastung zur Durchbiegungsfliche (28).  denkt (Abb. 5). Erstens ist damit ein Schnitt ge- 
ringsten Vorzeichenwechsels beziiglich der Last durch 
die Platte gelegt, und zweitens werden durch die Betrachtung eines Balkens zu groBe Durch- 
biegungen geliefert; beides sind also Umstande, durch die wir auf der sicheren Seite bleiben. 
Nun verhalten sich aber die maximalen Durchbiegungen eines Balkens, der einmal durch 
eine einsinnige konstante Last belastet wird und das andere Mal der Summe nach durch 
dieselbe konstante Last, jedoch derart, daB zu gleichen Teilen von der Mitte aus betrachtet 
auf der einen Seite positive und auf der anderen Seite negative Belastungen vorliegen (Abb. 5), 
wie 16 zu 1. Also wird die Wirkung von p héchstens 1/16 derjenigen sein, wenn p als einsinnige 
Last tber die Plattenflache beziiglich wy verteilt ware. 
Wir befassen uns nun mit der Verteilung von p tiber den Radius 9. Dazu bilden wir aus 
(30) eine wahre Kraft und erhalten 


dst = pdF = — 16 K cos2 ydyd(Ino). (30a) 


Die Anderung der Kraft St ttber den Radius 0 geht also logarithmisch vor sich, woraus wir 
erkennen, dafs in der Umgebung von 9 = 0 die Wirkung von p vernachlassighar klein gegenitiber 
derjenigen des Systems S, aus (29) ist. Wir kénnen demnach fiir die Wirkung p iiber 0 an- 
nehmen, daf} diese in 9 = 0 praktisch Null ist und fiir gréGBer werdende 9 anwichst, jedoch in 
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der Nahe des Randes bis zu Null wieder abnehmen wird, da eine Last in der Mitte einer Kreis- 
platte angebracht erheblich gréBere Wirkung, d. h. Durchbiegung, zeitigt, wie wenn dieselbe 
Last in der Nahe des Randes aufgebracht werden wiirde. Diese Umstinde legen es nun nahe, 
die wirklich auftretende ,,Gleichgewichtslast“ p mit einer Gewichtsfunktion g zu multiplizieren 
und somit eine neue ,,Ersatzlast‘‘ P =8&p zu bilden; die so erzeugte Funktion p soll solche 
Belastungsordinaten aufweisen, wie sie in erster Naherung der tatsdchlichen Wirkung von p 
als tiber den ganzen Radius der Platte gleichmafig verteilte Last entsprechen, so daB man sich 
durch Integration von p tiber einen Quadranten unter Beriicksichtigung der Wirkung von p 
als Gleichgewichtssystem tiber y ein Bild ither die tatsachliche Wirkung von p machen kann, 
wenn man die durch die Integration entstehende Summe als gleichmaBig verteilte einsinnige 
Last aufgebracht denkt. Nach dem eben Gesagten muB p in 9 = 0 verschwinden, woraus folgt, 
daB g wie 9° gegen Null gehen muB; wenn weiter p fiir @ = Q) Null sein soll, so muB g am 
Rande 9 = Q, ebenfalls verschwinden, so da® wir fiir g eine Parabel vierter Ordnung ansetzen 
werden, von der Form 


g = Co (a — 0). (2) 

In (32) ist C eine noch freie Konstante, mit deren Hilfe itber die GréBe des Maximalwertes 
von g zwischen 9 = 0 und 9 = Qy verfiigt werden kann. Mit (32) ist die Verteilung der Wirkun- 
gen von p tiber 9 festgelegt, und wir kénnen nun die in (32) noch freie Konstante C dazu ver- 
wenden, die weiter oben besprochene obere Schranke der tatsachlichen Wirkung von p beziiglich 
y zu erfassen. Wir hatten festgestellt, daB wegen der Verteilung von p iber y die tatsachliche 
Wirkung von p hichstens 1/16 derjenigen sein kann, wenn p iiber die ganze Flache einsinnig ver- 
teilt ware. Also werden wir das Maximum vong zu 1/16 wahlen, wodurch danng und somit auch p 


eindeutig festgelegt sind. Integrieren wir dann p im Quadranten von + bis + auf und 


denken uns die so erhaltene Summe als vierten Teil einer konstanten einsinnigen Flachen- 
belastung, so haben wir in erster Naherung eine obere Schranke, fiir die aus dem Vorhanden- 
sein der Flachenbelastung p aus (30) zu erwartenden maximalen Durchbiegungen. Mit 


max = 1/16 bestimmt man C aus (32), bildet mittels (30) p = gp, integriert von — + bis 


% 4 
uns 1/8 der scheinbaren Einzellast —P, von der weiter oben gesagt wurde, da sie auf der 
Linie y = 0 bzw. a in voller Héhe wirkt, auf die ganze Platte als konstante Last verteilt zu 
denken. Der Vergleich P einmal als Einzellast in 9 = 0 angreifend und das andere Mal P/8 
als konstante Last iiber die ganze Platte verteilt, liefert, da die maximale Durchbiegung in 
o =0 von P/8 als konstante Last rund 5% derjenigen von P als Einzellast aufgebracht, 
betragt. Damit kénnen wir nun abschatzend feststellen, daB die Wirkung der Flachenlast p 
aus (28) im ungiinstigsten Fall héchstens 5% derjenigen des Systems S,, wie es durch (29) 
gegeben ist, betragen kann. Wir werden also die Flachenlast p vernachlassigen kénnen und 
kommen zu dem SchluB, daB® mit sehr guter Naherung allein das System S, die Form der 
Durchbiegungsflache (28) erzeugt. 

Nun wollen wir von den Ergebnissen, die im konzentrischen Fall fiir S, gefunden wurden, 
auf den exzentrischen Fall schlieBen. Wenn man sich fir einen Augenblick als Beobachter auf 
der nach (28) gegebenen Durchhbiegungsflache befindlich denkt, auBerhalb des Mittelpunktes 
der Platte, und sich auf der Platte bewegt, so erkennt man, daf eine reine Entfernungsanderung 
vom Mittelpunkt die Gré8e der in 9 = 0 scheinbar wirkenden Einzellast Ps nicht andert, da 
man bei einer reinen Entfernungsanderung auf einer Linie y = konst. voranschreitet, und 
Ps = — Pcos2y vom Winkel, nicht aber vom Radius abhangt. Wenn sich jedoch vom 
Beobachter aus gesehen eine Drehung,des Mittelpunktes vollzieht, d. h. wenn man sich quer 
zu den Linien y = konst. bewegt, so wird sich eine Anderung der scheinbaren Einzellast Lid 
ergeben. Fir das nun Folgende mége sich der Beobachter auf der Linie y = 0 an einem 
beliebigen aber festen Punkt 9 =konst. aufhalten, und statt seiner bewege sich jetzt der Punkt 
o = 0, der Angriffspunkt von S, bzw. Ps aus dem Mittelpunkt lings der Linie y =a heraus. 
Da sich der Beobachter auf der Linie y = 0 aufhalt, entspricht dieser Vorgang einer reinen 
Entfernungsanderung zwischen S, und dem Beobachter. Da also eine reine Entfernungsande- 
rung von S, bzw. Ps vorliegt, befindet sich der Beobachter, auch wenn S, nicht in der Mitte 
der Platte, sondern exzentrisch angreift, langs der ganzen Linie y = 0 bzw. a, also fiir alle 0 
auf einer Durchbiegungslinie, die in erster Naherung durch eine konstante scheinbare Einzel- 


auf und multipliziert die so erhaltene Summe mit vier: Man erhalt a P; d. h. wir haben 
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last Ps erzeugt gedacht werden kann. Kurz gesagt haben wir damit das Ergebnis gewonnen, 
daB das Kraft- und Momentensystem S,, auch wenn es exzentrisch auf der Linie y = 0 baw. 2 
angreift, wie eine scheinbare Einzellast wirkt; qualitativ gesprochen erzeugt S, langs der 
Hauptdurchbiegungslinie, wie wir die Linie p = 0 bzw. a fernerhin nennen wollen, eine Einzel- 
lastdurchbiegungskurve. Da nun aber nach (31) S@) ahnlich zu S, ist, sind auch die durch die 
Systeme erzeugten Durchbiegungslinien ahnlich und wir kommen zu dem Ergebnis: Das 
Hauptstérsystem S() aus der Plattenlésung (17) wirkt allein in erster Naherung langs der ganzen 
Hauptdurchbiegungslinie wie eine scheinbare Einzellast im Quellpunkt r = 0 angreifend; d. h. 
die Wirkung von S() ist entlang der Linie p = 0 baw. a der gewiinschten durch ki} aus (24) 
hervorgerufenen Punktlastwirkung abnlich, d. h. bis auf einen konstanten Faktor gleich. 


9. Elimination des Hauptstérsystems. Die zunachst erzielten Aussagen kénnen sich also 
nur auf die Hauptdurchbiegungslinie beziehen. Zusammenfassend wissen wir bisher, daB die 
nach (17) fir y = 0 baw. x gegebene Durchbiegungskurve qualitativ derjenigen einer in r = 0 
angreifenden Einzellast naherungsweise entspricht; jedoch ist die GroBe dieser Einzellast noch 
nicht bekannt, da sich die Einzellastdurchbiegungskurve aus (17) aus zwei Komponenten 
zusammensetzt: erstens Wirkung von k{}) und zweitens Wirkung von S(?), die derjenigen von 
k\) auf der Hauptdurchbiegungslinie entgegengerichtet ist (Abb. 4). Wenn mit w, allein die 
Wirkung von ks} bezeichnet wird, so gilt bei Vernachlassigung der Wirkung des Nebenstér- 
systems S(!) wegen der Ahnlichkeit der Wirkungen von kG} und S(?) lings der Hauptdurch- 
biegungslinie folgende Ahnlichkeitsrelation: 


(33a) 


Hieraus folgt 
wp(ote 
Wp (r, ) se (0-3) w ( , ) ‘ (33) 


Falls also wy in r = 0 bekaunt ist, kann nach (33) die Wirkung von k(} allein durch Multipli- 
kation der Urlésung (17) mit einem Faktor ermittelt werden. Weiter gilt ganz allgemein aus 
Uberlagerungsgriinden bei Vernachlassigung der Wirkung des Nebenstérsystems, wenn wsi) 
die Wirkung des Hauptstérsystems bezeichnet, 


w(0, 0) = w,(0 , 0) + wsi(0 , 0) , 
10 (0,3) = ep (0.5) + ese (0.5). fe) 


Nun ist aber die Wirkung bzw. die Durchbiegung von ki), also der Einzellast, in der Um- 
gebung von r = 0 vom Winkel unabhangig, und es gilt w,(0, 0) = w,| 0, + - Gemaf dem 
Verlauf der algebraischen Singularitaét in r = 0 ist die Wirkung von S(@) auf p = 0 derjenigen 
aul @ = + entgegengesetzt, also wg2)(0,0) = — wgi){ 0, z): Ferner ist, wie ‘schon einmal er- 


wahnt, die Wirkung des Hauptstérsystems bei py = 0 derjenigen von ki’) entgegengerichtet, so 
da8 wir schlieBlich aus (34a) erhalten 


w(0,0) = w,(0,0) — wes2)(0, 0) , 


100.) = w9(0, 0) + 1059(0, 0). (34) 


Mit (34) und (33) entsteht nun fiir die Wirkung von ki) allein langs der Hauptdurch- 
biegungslinie 


nfo (2) SG 
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Nach (34) kann man weiter die GréBe des Anteils bestimmen, den die Hauptstérung gegeniiber 
der Wirkung w, ausmacht, und der vermége (35) in erster Naherung eliminiert worden ist. 
Nach (34) wird 
100 
— ; = wg) % Von wy. (35a) 
142 
C4 
Um die Lésung fiir die Hauptdurchbiegungslinie nach (17) explizit angeben zu kénnen, 
mu} abschlieBend noch die Konstante c, errechnet werden. Dazu bildet man 


100 wiz) 
Wp 


20 


EN dp ee, (36a) 


om 
Diese Summe muB nun in iiblicher Weise gleich der gewiinschten Einzellast P mit umgekehrten 
Vorzeichen sein, woraus folgt 
Piereals 
eas 


Wenn man weiter aus rechentechnischen Griinden anstelle der Konstanten c; nach (20) 


(36) 


¢=— (20c) 


einfithrt, so erhalt man schlieflich in erster Naherung gemaB (35) fir die Hauptdurchbiegungs- 
linie einer durch die Einzellast P exzentrisch beanspruchten und frei aufgelagerten Kreisplatte 
in exzentrosymmetrischen Parametern 


Ronee eter omit ee a ee 
wp (r. -s =TnKG@ +e) yada Zt — 1°) [ee + £65 2 (37) 


Fiihrt man in (37) r = oa @ ein, wobei R, den wahren Plattenradius bedeutet, und geht mit 
0 


ry —> 0 unter Beachtung, daB, die c; und 6 Funktionen von ry sind, so entsteht die bekannte 
Lésung fiir den konzentrischen Fall. 


10. Numerische Ergebnisse. Jetzt soll zwischen der unter (37) hergestellten Naherungs- 
lésung und der von E. Reifner angegebenen exakten Lésung ein numerischer Vergleich der 
Hauptdurchbiegungslinien durchgefiihrt werden. Die Lésung von E. Reifner einer am Rande 
frei aufgelagerten und durch eine Einzellast belasteten Kreisplatte in Polarkoordinaten lautet 


a—2z 


2 
wt (p, y= PEE Real |e — a? | logs 


| 


+ (1 a’) (1—9") 


az 
Yaz 


1+» d 
te w dt l—vy 
we) | tests 


(38) 


Hierin ist z = (0/R,) e' ¥ zusetzen, a Ry ist der Abstand vom Mittelpunkt, in dem die Einzellast P 
aufgebracht ist. 

Es werden fir den Vergleich von (37) und (38) drei Falle verschiedener Exzentrizitat ¢ 
betrachtet: ¢ =1/4 R,, e=1/2 Ry, €= 3/4 Ry ; hieraus errechnet sich ry, der Parameter des 
Plattenbegrenzungskreises, nach der MaBgabe ry) = ¢/Ry, was durch elementargeometrische 
Uberlegungen an Abb. 1 leicht zuerkennenist. Diese drei Falle sollen der Reihe nach mit a, f, y 
bezeichnet werden. Es werden durchgehend als konstant angenommen: erstens die Einzellast 
zu P = 2a K/b;, zweitens der Halbmesser der Platte mit Ry = 4 und drittens die Quer- 
kontraktionszahl zu y = 1/3. Die Formel (35a) liefert dann als prozentualen Anteil, den die 
Hauptstérung S@) gegeniiber der Wirkung von k({/} ausmacht, 


&) Ws2) = 4,6% von wy, | 
fp) Ws(2) = 19,5% von Wp > (39) 
v) ws) = 48,6% von w,. | 


Um soviel Prozent, wie in (39) angegeben, ware also die Hauptdurchbiegungslinie gemaB der 
Urlésung (17) falsch, wenn S@) nicht eliminiert worden ware. 


Ingenieur-Archiyv 
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Die Wirkung des durch die Nebenstérung hervorgerufenen singuléren Momentes gemaf (27) 
wird ein sehr geringer Bruchteil der Wirkung von S(@) sein, also in Prozenten ein geringer 


Bruchteil derjenigen Wirkung, wie sie durch (39) angegeben ist. 
Zum Vergleich der Funktionen (37) und (38) wurde der in (38) auftretende Integrand in 
eine Reihe bis zum siebenten Glied entwickelt und aufintegriert; dann wurde weiter der Real- 


teil der nichtanalytischen komplexen Funktion gebildet, wobei man erhalt: 


2 
e ogee 


§ cosyw+a 
PER el a 5 R2 R 
eOle y= gamle lag ee mete eee ar eee 
1 2 as cos y+a Re 
Pa a ea ar cael 2 e : 1 ilu Se I 8 hea ert 
(l—a V(t — Ke} eer 2(1l—a V(t R) |» + RG +H) cosy] (38a) 

GO) en. (a 9) (eee 
zm re Es) cos 2 y + R740) cos 3y + Ri (9-9) cos 4y 
ero ka Gene Cai eve el 

| FAYE TEESE Dae TRB ss By) oe Oar Ts em cos 7p + fe 


Hierin ist nun @ mit ry zu identifizieren und y =z baw. 0 zu setzen. Beide Funktionen sind 
an den Stellen 9 = 0, 1, 2, 3, 4 berechnet worden. Es wurden die in Tabelle I aufgefiihrten und 
durch Abb. 6 dargestellten Funktionswerte gefunden. 


Tabelle I. 
" R 5 Max. Abweich 
3 yg =nbzw.y = 0 ¢ = 0bew.yp =a (37) gepeultber (Ge) crete 
Fall nach Se oa a 2 auf an den Stellen: 
Formel = We Oss Que — Qe = Ore: 0 0 
ee 3|/ e=2|) @=1] e=0)| e=1] @ aise Ale: ae 3,9 =0|9=3,9=7 
(37) 0 | 1,30~.2,52 | 3,334 3,19 | 2,57 | 1,75 | 0,87 | 0 p 
o | Z g | {ie ORT TS 
38 OQ. | 1/31 12.52-453,99 | 3219") 258° 99, 77) 0n88 40 fo | 7o 
(37): file—20 4 BAW O40) 953 9:20 eT Gonadal And 0.56 0 | 
3 ae? , , 5,19 1,5° 
P 38). hy Ox i) 4828) 2,860). 2,52. 2,24 a1 77 | 21591) (0,50 0 7o 7o 
(37) 0 | 1,20 | 1,43 | 1,30 | 1,06 | 0,80 | 0,52 | 0,25 0 : 
» | , , 13,89 8,82 
(38) 0 | 1,01 | 1,32 | 1,30] 1,13 | 0,88 | 0,59 | 0,29 0 yo | 18.8% 


Abb, 6. Vergleich der Hauptdurchbiegungslinien. ——— Naherungslésung, — — — exakte Lésung. 


Wie man erkennt, ist die Ubereinstimmung der Durchbiegungskurven fiir kleine Exzentrizi- 
taten sehr gut; von Exzentrizitaéten gréBer 1/2 Ry ab diirfte jedoch die Brauchbarkeit der 
Lésung (37) in Frage gestellt sein. Andererseits ist zu bemerken, da® auch die praktische Ver- 
wendbarkeit der Lésung (38) nach E. ReiBner immer geringer wird, je gréBer die Exzentrizitat ist 
da die Reihe in (38a) um so schlechter konvergiert, je gréBer a = ry ist. Fir den Fall y ans. 
um die Genauigkeit der Tabelle I zu erreichen, sechs Glieder der Reihe beriicksichtigt ‘recdem 
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11. Extrapolation von den Ergebnissen fiir die Hauptdurchbiegungslinie auf die ganze 
Durchbiegungsfliche; Hebung der algebraischen Singularitat. Won der Ahnlichkeitsrelation 
(33a) wissen wir, daB sie, falls die Exzentrizitat nicht zu groB ist, mit sehr guter Naherung 
fir p = 0 baw. gilt; dartiber hinaus nehmen wir an, daB fiir alle Linien gy = konst. (Abb. 1) 
des exzentrosymmetrischen Bezugssystems diese Ahnlichkeitsrelation bestehen mége. Dann 
wird aus der so erweiterten Ahnlichkeitsrelation (33a) 

wplr, ~) __ wp(0, @) 

w(r,~)  w(0,q) uy 
Mit diesem Ansatz fallt die gesuchte Funktion w,(r, p) langs der Hauptdurchbiegungslinie 
mit der bereits hergestellten Naherungslésung (37) zusammen, und aufBerdem enthalt w,(T, P) 
aus (40a) fiir r = 0 keine unendliche Vieldeutigkeit mehr. Dies leuchtet ein, wenn wir aus 


(40a) bilden 


Ww p(0, wp) 
wp(T, Y) = w(0, y) w(r, p) Cy 


In (40) mu8 noch w,(0,~) auf der rechten Seite eliminiert werden, damit w,(r, y) allein 
durch Bestandteile der Urlésung (17) ausgedriickt ist. Nun gilt analog zu (34a) ganz allgemein 
aus Uberlagerungsgriinden fiir alle p 


w(0, Y) ay w,(0, ~) =F ws(2)(0, 7) ’ 
IU TU 


w (0 P+ 5 =, (0, Vor a) + Ws(2) (0. @ + 3): 


wobei wiederum wg) die Durchhangung in r = 0 sein soll, die allein durch die Wirkung des 
Hauptstérsystems entsteht. Ebenso wie schon einmal im Abschnitt 9 ist nun zu beriicksichtigen, 


} 
| 
[ (41a) 
J 

da w,(0, p) = w,(0,y + > , und daf weiter wse)(0, p) = — wg) ¢ Pp 7) wegen des Ver- 


laufs der algebraischen Singularitat tiber ¢y ist, so daB aus (41a) entsteht 


w(0, p) = w,(0, y) + wse(9, ~) , 


———-— aS 


(0.9 + -F) = wol0, g) — wesel0, 9). | oe 
Mit (41) und (40) erhalt man fiir die Durchbiegungsflache w, 
w(0, ~) + (0. Cis Zl 
WAT PY = 2 wl0, @) w(r, @) - (42a) 


Diese Darstellung kann noch vereinfacht werden, denn nach (17) ist 
w(0, ») + (0, p+ Z| = 1? [2 cy + 4] = w(0, 0) + w (° aki (42b) 


Damit wird aus (42a) 


1wp(r, 9) => (0, 0) + w (0. 5) a (42) 


Die so konstruierte Funktion w,(r,q) werden wir nun nachtraglich auf ihre physikalische 
Bedeutung hin zu interpretieren haben, damit die Berechtigung der Voraussetzung, die zur 
Konstruktion von w,(r, 7) fihrte, bestatigt werden kann. 

Mit (17) wird aus (42), wenn wir wieder ¢; nach (20c) einfiithren und c, gemaB (36) einsetzen, 


PU? 2 2 2Gr 4, 


‘ r D\ es bn eal a 
w,(r, ~) = TR He oe oa r? In > + (r2 — r?) (¢, + €, cos? p + ¢,rcosq)|. (43) 


12. Physikalische Deutung der nach (43) gegebenen Funktion. Durch Einsetzen von (43) 
in (19) bestatigt man leicht, daB mit (43) die Durchbiegungsflache einer am Rande frei auf- 
gelagerten Kreisplatte gegeben ist. Wir wenden uns jetzt den Belastungsverhaltnissen zu, wie 
sie durch (43) geliefert werden. Zuerst stellen wir nochmals fest, da (43) im gesamten Platten- 
bereich wegen des angebrachten Zusatzfaktors 1/(c, + c, cos’ y) singularitatenfrei ist, und 
also kein Punkt der Platte ausgeschlossen oder gesondert betrachtet zu werden braucht. Da 
die Auflagerbedingungen durch den Zusatzfaktor nicht gestért worden sind, ist anzunehmen, 
daB (43) im Gegensatz zu (17) nicht die homogene Plattengleichung 44w = 0 erfiillt. Flachen- 
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hafte Belastungen p(r, y) werden bekanntlich fiir die Theorie der diinnen Platte durch die 
inhomogene Plattengleichung wie folgt geliefert: 


KAAw,(r, ¢) = p(t 9) » (44) 
wobei K die Plattensteifigkeit aus (6) ist. Bevor jedoch die flachenhafte Last p(r, y) behandelt 
wird, soll (43) auf singulare Lasten untersucht werden; denn laut Konstruktion miissen wir in 
r — 0 eine Einzellast annehmen. Demnach bilden wir mit (43) an einer beliebigen Stelle 
r = konst. die Summe aller Querkrafte von 0 bis 27 und bezeichnen die resultierende ent- 
gegengerichtete Kraft mit Pp(r). Mit (9), (5), (4) und (11) erhalt man 


OS Oy ae ky 

P(r) = 2 4 p — 45 

a) 2V cg-(cg + ¢,) ae ( ) 

Aus (45) lesen wir ab, da®B Pz keine Funktion von r ist, woraus folgt, daB die Flachenlast P(r, y) 

der Platte auf jedem beliebigen Ring von der Breite Ar die resultierende Null hat, wegen 

Pa(r + Ar) — P(r) = 0. Weiter erkennen wir aus (45), daB in r = 0 die Einzellast Pp 

wirkt. Das ist aber gerade der gewiinschte Belastungsfall. Der Faktor x aus (45) ist deutbar 
als 


w(0,0) + u( a 


4 
V w(0, 0) uo. =) 


d. h. man dividiere das arithmetische Mittel der Funktionswerte der Urlésung (17) in r = 0 


5 (45a) 


fur gy — 0- und -g = + durch das geometrische Mittel derselben Funktionswerte. Wir wissen 


im vornherein, daf das geometrische Mittel stets kleiner als das arithmetische Mittel ist, und 
zwar unterscheiden sich beide Mittelwerte um so weniger, je dichter die Bezugszahlen liegen. 


Da die GréBe des Unterschiedes von w(0, 0) und w{0,—)durch die GréBe des Stérsystems 


S() festgelegt wird, erkennt man, daf x sich fir kleine Exzentrizitaten wenig von 1 unter- 
scheiden wird, da das Hauptstérsystem S@) mit verschwindenden Exzentrizitaten wie ri gegen 
Null geht. Fiir die in Abschnitt 10 behandelten Fille a, f, y ergibt sich z. B. fiir x: 

Hai Uns % = 1,019, Hy, = 1,140. (45b) 
Mithin ist fiir nicht zu groBe Exzentrizitaten, auf die der Anwendungsbereich der Funktion (43) 
ohnehin beschrankt ist, die nach (43) gelieferte Durchbiegungsflache einer frei aufgelagerten 
Platte mit hinreichender Genauigkeit durch die Einzellast P belastet. 

Die Wirkung der nach (44) zum gewiinschten Belastungsfall zusatzlich auftretenden ,, Gleich- 
gewichtsbelastung’’ muB von vornherein als gering angenommen werden, da die Hauptdurch- 
biegungslinie aus (43), die auf Lasten am ehesten reagiert, mit der exakten Hauptdurchbie- 
gungslinie nach ReiBner gemaf (38) sehr gut ibereinstimmt. Um diese Tatsache auch numerisch 


zu belegen, wurde die ,,Gleichgewichtslast“‘ (44) von y = 0 bis = uber die Plattenflache auf- 


integriert, wobei entsprechend dem Prinzip von de Saint-Venant um den Punkt r = 0 ein Kreis, 
dessen Flacheninhalt 5% der gesamten Plattenfliche ausmachte, bei der Integration weg- 
gelassen wurde, wodurch also nur solche Zonen der ,,Gleichgewichtsbelastung“ erfaBt sind, 
von denen man annehmen kann, daB sie auch als ,,Gleichgewichtssystem“ fahig sind, merkbare 
Durchbiegungen zu erzeugen. Die so entstehende resultierende Kraft W der », Gleichgewichts- 
belastung” wurde fiir die drei in Abschnitt 10 behandelten Falle berechnet zu 


W,=0,002 P, W,=0,093P, W,=0,646 P. (46a) 

Die Werte (46a) dividiert durch die zwischen y = 0 und => eingeschlossenen Flachen- 
inhalte liefern die folgenden mittleren Spannungen om: 

Om, = 9,00019 P, On,g = 90,0056 P, on, = 0,0286 P. (46) 


In Abb. 7 sind die aus der ,,Gleichgewichtsbelastung* p(r,y) resultierenden mittleren 
Spannungen o,,/P, wie sie nach (46) geliefert werden, in Abhangigkeit vom Parameter r, dar- 
gestellt. Miteingezeichnet wurde die mittlere Spannung o,/P = 1/m R?, wie sie bei gleich- 
maBiger Verteilung der Einzellast P iiber die ganze Plattenflache auftreten wiirde, und man 
erkennt, da die Wirkung von oy im Vergleich zu derjenigen von P als Einzellast im ungiinstigsten 
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Fall etwa 40% betragen kann, wie ein einfacher Vergleich im konzentrischen Fall lehrt, daB bis 
zu Ty = 0,5 die durch die »,Gleichgewichtsbelastung auftretende Stérung des Einzellast- 
falles vernachlassigbar klein ist, zumal die hier durchgefihrte Abschatzung sehr grob ist und 
weit tiber den tatsachlich auftretenden Stérungen liegen wird. 


Um die Wirkung von p(r,q) als auBere Last zu klaren, blieben gemaf} dem Prinzip von 
de Saint-Venant die Flachenlasten innerhalb eines kleinen Kreises um r — 0 unberiicksichtigt. 
Es ist jedoch méglich, wenn die Flachenlasten um r = 0 auch nicht als iuBere Lasten wirksam 
werden kénnen, da sie ebenso wie das Stérsystem S() aus (24) ein singulares Moment in r = 0 
erzeugen. Es la8t sich jedoch leicht bestatigen, da® dies nicht der Fall ist. 


Damit ist gefunden, da die vorausgesetzte erweiterte Ahnlichkeitsrelation (40a) mit 


groBer Naherung fiir nicht zu grofe Exzentrizitaten bis etwa ry = =; richtig ist und mit (43) 


eine recht gute Naherungslésung der Durchbiegungsfliche einer am Rande frei aufgelagerten 
und exzentrisch durch die Einzellast P belasteten Kreis- 
platte konstruiert worden ist. 


0 G25 OS GIS 10 Abb. 8, Durchbiegungslinien unter dem zur Hauptdurchbiegungslinie orthogonalen 


Abb. 7, Zunahme der Stérwirkung nach (44) Radius der Platte. 
mit wachsender Exzentrizitat. Naherungslésung, — — — exakte Lésung. 


13, Numerische Ergebnisse. Mit (43) ist nun die Durchbiegungsflache fir alle ~ gegeben, 
und so wie bei Funktion (37) ein Vergleich mit der exakten Lisung durchgefiihrt wurde, soll 
jetzt ein Vergleich zwischen der Reifnerschen Lisung (38) und der Naherungslésung (43) langs 
des zur Hauptdurchbiegungslinie orthogonalen Durchmessers der Platte angestellt werden. 

; Wegen der Symmetrie der Tabelle IL. 
Lésungen reicht es aus, nur einen 


Radius, der unter dem ortho- €n/ Ro 
gonalen Durchmesser entstehen- Fall 0 = . = 1 
den Durchbiegungslinie zu _ be- Pama 
trachten. Dazu ist in (38) bzw. 
it ‘a (43) 3,19 2,84 2,04 1,01 0 
(38a) wy — 2 zu setzen. Es Wwure- (38) 3,19 2,84 2,03 1,02 0 
ae es eee vee 8 (43) | 2,20 2,02 149 0,7 0 
cr nator 1, 4/3, (38) 2,24 2,03 1,47 0,74 0 
2, 4, oo. Es wurden nun wieder- 2 
um die in Abschnitt 10 behan- (43) | 1,06 1,00 0,77 0,41 0 
delten Beispiele zugrunde gelegt, ? (38) Us 1,03 0,75 0,38 0 


und man erhalt fiir die Durch- 4 
biegungen die in Tabelle II angegebenen Werte, deren gute Ubereinstimmung in Abb. 8 er- 


sichtlich ist. 


14, Dimensionierung exzentrisch durch Punktlast belasteter Kreisplatten. Bisher wurden 
nur Durchbiegungen behandelt, und es ist mit den gefundenen Ergebnissen méglich, an jedem 
Punkt einer frei aufgelagerten Kreisplatte, die exzentrisch durch eine Einzellast belastet ist, 
fir nicht zu gro%e Exzentrizitaten die Durchbiegung mit groBer Genauigkeit anzugeben. 
Weitaus wichtiger als die Dimensionierung nach den Verformungen diirfte jedoch die Bemes- 
sung nach Momenten sein. Wie bekannt, liefern die Punktlastlésungen dinner Platten unend- 
lich groBe Momente unter dem Lastangriffspunkt. Fir gewéhnlich geht man nun so vor, um 
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diese Lésungen dennoch zur Momentendimensionierung heranziehen zu kénnen, das man die 
GroBe des Momentes in einer solchen Entfernung s,, vom theoretischen Lastangriffspunkt 
betrachtet, wie sie dem tatsdchlichen Halbmesser des Laststempels entspricht. Andererseits 
muB gefordert werden, da® fiir die sinnvolle Verwendung einer Punktlastlésung der Halb- 
messer s,, des Laststempels im Verhaltnis zum Plattendurchmesser sehr klein sei. Deshalb 
kénnen wir uns bei der Berechnung des Dimensionierungsmomentes M,,) darauf beschranken, 
die Untersuchungen in der Umgebung des Punktes r = 0 zu fiihren, was durch einen Nullkreis 
iiber dem Funktionssymbol bzw. durch den Index ¢ bei r angedeutet werden wird. 


Um das Moment M,, in der Umgebung von r = 0 zu berechnen, ziehen wir die singulari- 
taitenfreie Lésung (43) heran. Diese gibt die Durchbiegungsfliche einer am Rande frei auf- 
gelagerten Kreisplatte, die nach (45) durch die Einzellast x P belastet ist und auBerdem zu- 
sitzlich durch ein ,,Gleichgewichtssystem’ von Flachenlasten beansprucht wird, dessen 
Wirkungen sich als vernachlassigbar klein gegeniiber derjenigen der Punktlast herausgestellt 
hatten. Wir bemerken zuerst, daB das durch eine Einzellast hervorgerufene Moment Me in der 


Umgebung von r = 0 mit hinreichender Genauigkeit auch im exzentrischen Fall keine Funk- 


tion vom Winkel ist, daim Kleinen Zylindersymmetrie herrscht. Mithin ist also Mo invariant 
gegen Integralmittelwerthildung iiber y, und es gilt _ 


| MP ds, = MP. (47a) 
f dsy 3} 
0 


Wenden wir die Integralmittelwertbildung auf das aus (43) zu bestimmende Moment M,, an, 


° 
so ist einerseits das gesuchte Moment Mi) gegeniiber dieser Operation invariant, und anderer- 


ce} 
seits werden diein M,,) noch enthaltenen geringfigigen aus der ,, Gleichgewichtszusatzbelastung 
herrithrenden Stérungen im wesentlichen eliminiert, eben deshalb, weil ein ,,Gleichgewichts- 
system“ iiber y um Null alternierende Spannungs- und Momentenzustande erzeugt. Dariber 
hinaus hat die Formel fir den Integralmittelwert des Schnittmomentes eine erheblich ein- 
fachere Gestalt als der Differentialoperator, der das Moment selbst liefert. Da die Dimen- 
sionierungsformeln fiir die Schnittmomente dasjenige Ergebnis jeder theoretischen Platten- 
untersuchung darstellt, auf das die Praxis letztlich am meisten Wert legt, soll in diesem Fall 
der Faktor x aus (45) nicht identisch gleich 1 gesetzt werden, sondern seine wenig von 1 ver- 
schiedenen Werte werden, was besonders fiir gréBere Exzentrizitaten eine Verbesserung bedeu- 
tet, exakt beriicksichtigt, was physikalisch soviel bedeutet, da die Platte nicht wie nach (43) 
durch xP belastet wird, sondern genau mit P. Um diesem Umstand Rechnung zu tragen, 
brauchen wir die Lésung (43) nur mit 1/x durchzumultiplizieren, oder einfacher, da x eine Kon- 
O° 


stante ist, wir multiplizieren nachtraglich M,,, mit 1/x. So entsteht schlieBlich die Formel fiir 
das gesuchte Dimensionierungsmoment 


pee ee Mo) 
ues 2 B jee é (47) 
0 


worin M,,) das nach (43) zu berechnende Moment und x der durch (45) bzw. (45a) gegebene 
Faktor ist. Geht man mit dem Differentialoperator nach (9) fiir das Moment M,, unter Zu- 


hilfenahme von (10) in (47) ein und beachtet (=) =1-+2r,cos » +.---, so erhalt man 
unter Weglassung der Glieder, die fir r — 0 nach Null gehen, nachdem U und V gemaB (4) 
ersetzt und aufintegriert wurde, 

Me = = EV eg) in ee vy Lee ce eel ep ces 

i) 8x Feit | rey pre ena witines = freee A 


0 Cy 


Nimmt man hierin nachtraglich x wieder gleich 1, so entsteht fiir kleine Exzentrizitaten die 
noch einfachere Formel 


MP = — 2 (1 +>)In2@—r2(1—») (142 a). (48a) 
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Der Vollstandigkeit halber und zum Vergleich ist auch das Dimensionierungsmoment fiir 
den Fall der eingespannten Platte berechnet worden. Man erhalt! aus (47) mit x = 1] 


aie P re 
elles tiga | ol Poe) : + 2+ (1+) r?}. (49) 


Die Abklingung der Momente bei wachsender Exzentrizitat ist fiir den frei aufgelagerten und 
eingespannten Fall, wie sie durch (48) und (49) geliefert wird, fiir zwei verschiedene Quer- 
kontraktionszahlen in Tabelle III angegeben und in Abbildung 9 dargestellt, wobei fir die 
Berechnung jeweils ein Parameter r, so gewahlt ist, daB alle Momente in einer Entfernung von 
r = 0 berechnet sind, die 3°% des Plattenradius betragt. 


Tabelle III. 


| | SET | ee ee 


| Formel | ha =_0 r= 5 n=5 maa ip = 1 
sok Lp shy 09,38 Oe ee (a RC TOMO 
Bh (49ye|  -7:35 7,09 6,25 4,40 0 
ante | (48) | 7,01 6,86 6,34 510 |: 6 
(49) 5,01 | 4,82 4,19 2,80 0 


15. Zusammenfassung. Behandelt wurde das Problem der durch eine Einzellast exzentrisch 
belasteten und am Rande frei aufgelagerten Kreisplatte. Es wurden dem Problem angepafbte 
Parameter, als exzentrosymmetrische Koordinaten bezeichnet, eingefiihrt. Durch invariante 
Schreibung der SchnittgréBen gelang es, A, 
diese sofort fiir exzentrosymmetrische Ko- Mr) 
ordinaten anzugeben. Unter Zugrundelegung 
eines bekannten Satzes itiber die Lésungen 
der homogenen Plattengleichung in Polar- 
koordinaten wurde durch Transformation 
eine Darstellung aller Lésungen der homo- 
genen Plattengleichung in exzentrosymme- 
trischen Koordinaten gewonnen. Durch Vor- 
gabe der vier einfachsten Potentialfunktionen 
wurde dann eine Funktion gefunden, die ge- 
eignet war, die Randbedingungen der freien 
Auflagerung zu befriedigen, jedoch im Ein- 
zellastangriffspunkt eine algebraische Sin- 
gularitat enthielt. Durch physikalische 
Deutung dieser Singularitat konnte eine 
Naherungslisung fiir die Hauptdurchhbie- 
gungslinie konstruiert werden mit der Ten- 
denz, mit wachsender Exzentrizitat an Ge- 9 025 O5 O15 10 
nauigkeit abzunehmen. Weiter wurde durch Abb. 9. Momentenabklingung mit wachsender Exzentrizitit. 
Verallgemeinerung der Methoden, die sich nach (48), — — — nach (49). 
zur Konstruktion der Naherungslésung fiir 
die Hauptdurchbiegungslinie bewahrt hatten, eine singularitatenfreie Lésung fiir die geé- 
samte Durchbiegungsflache aufgefunden. Die singularitatenfreie Lésung fiir alle p war dann 
geeignet, die zur Dimensionierung notwendigen Momente zu ermitteln. 


(Eingegangen am 19, Marz 1954). 


Anschrift des Verfassers: Dr.-Ing. Karlheinz Nasitta, Technische Universitat Berlin-Charlottenburg, 
Lehrstuhl fiir Mechanik Prof. Dr.-Ing. I. Szabo. 


1 Siehe FuBnote 1 von S. 89. 
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Uber spezielle Systeme von Tragfliigelgittern IL. 


(Theorie der diinnen Profile.) 


Von. K. Nickel. 


1. Einleitung. In einer I. Abhandlung 1 wurde iiber Systeme von Fligelgittern berichtet, 
die sich aus unendlich vielen ,,tragenden Linien“* aufbauen. Die vorliegende Arbeit behandelt 
die ebene Strémung durch endlich viele Gitter aus diinnen Profilen (Abb. 1 und 6). Die beiden 
Teile I und II zusammen erweitern die Ergebnisse des Verfassers iiber Systeme von Einzel- 
fligeln ? auf Systeme von Fliigelgittern. 

Die Voraussetzungen tiber das strémende Medium sind dieselben wie in I: es wird reibungs- 
freie, stationdre und inkompressible Potentialstrémung angenommen. Dabei kann die letzte 
Einschrankung in gewissem Umfange fortfallen, da man eine kompressible Unterschallstrémung 
mittels der Prandtlschen Regel auf eine aquivalente inkompressible Strémung abbilden kann °. 
Um die Formelausdriicke iibersichtlich zu erhalten, wird jedoch im folgenden auf die Aus- 
fihrung dieser Transformation verzichtet. 

Es werden nicht samtliche méglichen Systeme von Profilgittern behandelt, sondern nur die 
spezielle Anordnung, daf alle Gitter parallele Gitterfronten besitzen; zudem ist dies auch der 
einzige in den Anwendungen auftretende Fall. Dariiber hinaus miissen zur Vereinfachung der 
mathematischen Behandlung noch gleiche Gitterteilungen t = |h| (Abb. 1) vorausgesetzt 
werden. Diese weitere Einschrankung ist allerdings in der Praxis schon nicht mehr immer 
erfillt (z. B. Leitgitter mit vorgelagerter Stiitzenreihe). SchlieBlich besitzt ein Teil der mit- 
geteilten Ergebnisse sogar nur unter der noch viel einschneidenderen Bedingung Giiltigkeit, 
daB die Gitter alle senkrecht (Gleichdruckgitter) oder waagrecht stehen und die Schaufeln 
hintereinander oder ,,auf Licke“« angeordnet sind. Die im folgenden angewandte Methode ist 
dieselbe wie frither 2, d. h. man ersetzt die Profile durch eine Wirbel-, bzw. Quell-Senken- 
Belegung auf ihrer Sehne. Damit kénnen dann auch nur kleine Profildicken und -Wélbungen 
zugelassen werden. 

Fir diese, z. T. schwerwiegenden, Einschrankungen und Vernachlassigungen tauscht man 
dann allerdings den Vorteil ein, geschlossene Ausdriicke fiir Druckverteilung, Auftrieb und 
Moment anschreiben zu kénnen. Weiter lassen sich mit gewissen Modifikationen die Rechen- 
verfahren von Riegels und Truckenbrodt samt ihren ein fir allemal berechneten Konstanten 
ibernehmen. Fir gewisse Gittersysteme aus Streckengittern sind iiberdies die erhaltenen Er- 
gebnisse exakt giiltig, steht die verwendete Singularitatenmethode also in Ubereinstimmung 
mit der konformen Abbildung. 

Die vorliegende Abhandlung ist wie folgt gegliedert: In Ziff. 2 werden die Grundgleichungen 
aufgestellt und daraus die Stérungen im Unendlichen berechnet. Es werden die oben erwahnten 
Vernachlassigungen in die Grundgleichungen eingefiihrt und weiter die beiden Typen des ,,ge- 
wohnlichen* und des ,,alternierenden* Gitters in Sonderfallen als ebene Strahlstrémungen ge- 
deutet. Ziff. 3 enthalt die wichtigsten Formeln fiir die Behandlung der ersten und der zweiten 
Grundaufgabe im Falle eines einzelnen Gitters (n = 1), Ziff. 4, dasselbe fiir ein System aus 
n Profilgittern (7 > 1), auBerdem werden dort einige wichtige Folgerungen fir Skelettprofile 
zusammengestellt. In Ziff. 5 schlieBlich werden zuerst die exakten Lésungen im Falle eines 
speziellen Doppelgitters aus Streckenprofilen mitgeteilt. AnschlieBend werden die Formeln 
aus Ziff.3 angewandt, um die Truckenbrodtsche und Riegelssche Theorie auf Profilgitter zu 
iibertragen, wobei noch gepriift wird, wieweit die Giltigkeit einer linearen Theorie fiir endlich 
dicke Profile reicht. 


2. Vorbemerkungen. Die folgenden Betrachtungen sind allgemein fiir Systeme aus beliebig 
vielen Profilgittern giltig. Um die Ubersichtlichkeit zu wahren, werden sie jedoch nur fiir den 
Fall eines einzigen Gitters (n = 1) formuliert. 


1K, Nickel, Ing.-Arch, 22 (1954), S. 108,im folgenden kurz als I zitiert. 

2 K. Nickel, Ing.-Arch, 20 (1952), S. 363. 

* Man vergleiche dazu P. Byrd und M, Huggins, Ing.-Arch. 21 (1953), S. 191, sowie H. Schlichting, 
Z. Flugwissensch. 1 (1953), S. 109. 
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a) Ein Profilgitter in der komplexen z-Ebene sei durch seine Profilkontur C und die Gitter- 
konstante k = a/h = z/(h, + ih,) gegeben (Abb. 1). Man denkt sich die Strémung U +i V 
in der z-Ebene zusammengesetzt aus einer parallelen Grundstrémung wu, (cosa + isin a), 
der noch eine Stérstrémung u + iv ttberlagert ist. Diese Stérstrémung, die also den Einflub 
des Fligelgitters wiedergibt, sei erzeugt durch eine Wirbelbelegung y(z) auf den Konturen C 
und berechnet sich fiir die Stelle z aus der Formel 1 


u(z) = tr(z) = 2x {70 ctg k (2 —C) de. (1) 


Die Starke y der Wirbel unterliegt dabei der Bedingung, da die Profilkontur nicht durch- 
strémt werden darf, daB also die Strémung 


U(z)-+ 1 V(z) = u,, (cosa + isin a) + u(z) + iv(z) 


fiir alle Punkte z auf C die Richtung der Tangente von C besitzen mu. Sei etwa y(z) der 
Richtungswinkel der Profiltangente im Punkte z gegeniiber der reellen Achse (Abb. 1), so 
lautet diese Umstréimungsbedingung fiir alle z auf C 


u,, sin [x — y(z)] = cos y(z) Im ie ¢ vi cig bz = 0) ath + 


[ik . eae ee (2) 
+ 1) Re le G10) tg k (= —¢) ae} 


Da hier der Aufpunkt z auf dem Integrationsweg C liegt, sind die beiden Integrale singular. Sie 
sind als Cauchysche Hauptwerte zu nehmen, was durch das Symbol ¢ gekennzeichnet werden soll. 
b) Im Gegensatz zur Umstrémung 
eines Einzelprofils verursacht ein Pro- 
‘filgitter im Unendlichen eine endliche 


Ablenkung der Grundstrémung. Um 
diese zu bestimmen, geniigt es offenbar, C 
die Stérungskomponente fiir z = + ov, 


also allein in Richtung der reellen 
Achse zu betrachten, da iiber die Lage 
des Achsenkreuzes zum Gitter nichts - 
vorausgesetzt war. Es gilt fiir alle h, + 0 


lim ctgkz=+1 


z—+>-- co 


und also in (1) 


j : oi BN) 
Damit findet man fiir h, = 0 = f 
IRE ne g Ny 
Gy fo) 


[U +iV],.. =u.cos a 


, : Pooh, ae 
Ae et Abb. 1. Profilgitter. 
+1 (u- sin a —-- 2 he a) (3) ; ct 


mit der Profilzirkulation [ = | y(Z) df oder in Worten: 
c 


Die induzierte Zusatzgeschwindigkeit des Gitters im Unendlichen ist von der Gestalt des 
| Gitterprofils unabhangig und hangt nur ab von der Profilzirkulation und der Gitterkonstanten. 

Formuliert man dies Resultat gleich im allgemeinen Falle fiir ein System aus n Schaufel- 
gittern, so ergibt sich: 


1 Die Integration iiber die geschlossene Kontur C mége hier und im folgenden immer im mathematisch 
positiven Sinne erfolgen. 
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Jedes System aus n Schaufelgittern mit den Zirkulationen J’, und der einheitlichen Gitter- 
konstanten k hat im Unendlichen dieselbe Wirkung wie ein einziges Gitter aus Einzelwirbeln 


nm 

mit der Zirkulation /’ = 5) J’, und der Gitterkonstanten k. Die Strémungskomponenten weit | 
y=l1 
vor bzw. weit hinter dem System sind durch Gleichung (3) gegeben. | 
Im Gegensatz zur Umstrémung bei einem einzelnen Profil (oder auch bei endlich vielen), 
wo im Unendlichen die Parallelstrémung u. (cos a + i sin x) herrscht, kann man hier u~ also 
nur als ,,mittlere** Anstrémgeschwindigkeit im Unendlichen und ebenso a nur als ,mittleren“ 

Profilanstellwinkel deuten, sie hangen nach (3) itber die Formeln 
woe (ULSAN yee (Vege ae 2 Us OS VS eS (4) 
eindeutig mit den Geschwindigkeitskomponenten vor bzw. hinter der Gitterfront zusammen. 
Dieser Sachverhalt ist bei der Wahl eines Bezugsanstellwinkels und eines Bezugsstaudrucks zu | 
beachten. Die folgenden Resultate sind einheitlich auf u.. und « bezogen. | 
c) Zur Berechnung einer gesuchten Gitterstrémung dient die Integralbedingung (2). Je 
nachdem, ob darin die Wirbelstarke y(z) oder die Kontur C gegeben ist, spricht man von der 

ersten oder zweiten Hauptaufgabe ': 


Erste Hauptaufgabe. Gegeben: Zu- und Abstrémung [also nach b) auch ux und a] 
sowie y(z). 
Gesucht: k und C. 


Zweite Hauptaufgabe. Gegeben: Zustrémung, k und C, 
Gesucht: y(z) und Abstrémung. 


Beide Aufgaben bieten in dieser Allgemeinheit groBe Lésungschwierigkeiten. Beim ersten 
Problem sind in (2) die Integrale langs der Kontur C zu erstrecken, die jedoch gerade erst be- 
stimmt werden soll. Beim zweiten Problem mu in Gleichung (2) eine Fredholmsche Integral- 
gleichung erster Art aufgelést werden, fiir die ein lésender Kern im allgemeinen Falle bis heute 
nicht bekannt ist. Es ist daher angebracht, die Gleichung (2) zu vereinfachen. Hierzu werde 
die Kontur C fiir alle Profile durch eine gerade Strecke < a,b >, die Profil,,sehne‘ ersetzt. 
Damit geht (2) tiber in 

b 


uc sin [x = y(x)] = cow y(x) Im Fare) Cte hin 2) az\ kis | 
ao : (2a) 
+ sin y(x) Re Fee y(é) ctg k (x — é) as| a 


| 


wenn man o. B. d. A. annimmt, da die Strecke < a,b > auf der x-Achse gelegen sei. Im 
folgenden wird also eine spezielle Lage des Achsenkreuzes gewahlt, bei der die reelle Achse 
eine der Profilsehnen verlangert. 
Diese Vernachlassigung entspricht 
as) =. a ee aaa genau dem Vorgehen bei Einzelprofilen, 
* wo man das Profil in eine Skelettlinie 
und einen Profiltropfen aufspaltet 
(Abb. 2) und das Skelett durch eine 
Wirbelbelegung, den Tropfen durch eine Quell-Senken-Belegung auf der Profilsehne ersetzt2. 
Fir die Berechnung eines Profilskeletts geht damit die erste Aufgabe tiber in die Aus- 
wertung eines singularen Integrals, die zweite Aufgabe in eine gegeniiber friiher stark ver- 
einfachte Fredholmsche Integralgleichung erster Art. Fiir die Berechnung des Tropfens liegen 
die Verhaltnisse gerade umgekehrt. Durch die obige Vereinfachung hat man zwar schon bei 
den Grundgleichungen einen gewissen Fehler gemacht, doch zeigen Kontrollrechnungen, dab 
der Fehler im Endresultat fiir diinne, wenig gekrimmte Einzelprofile geniigend klein ist und 
da man ziemlich alle praktisch vorkommenden Falle mit dieser Methode behandeln darf. 


Profilkontur = Skelettlinie i Tropten 


+—— | ——_» 


Abb. 2, Aufbau cines Profils aus Skelettlinie und Profiltropfen. 


‘Vel. H. Schlichting, FuBnote 3 von 8.102, dochist in dieser Arbeit bei der Definition der ersten Haupt- 
aufgabe ein Irrtum unterlaufen. 

® Vel. F. Riegels, Ing.-Arch. 16 (1948), S. 373 und Ing.-Arch. 17 (1949), S. 94; E. Truckenbrodt, 
Ing.-Arch. 19 (1951), S. 365. 
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Sehr viel ungiinstiger liegen dagegen die Verhiltnisse bei Profilgittern. Einerseits hat man 
in der Praxis gerade auch an stark gekriimmten Schaufelprofilen groBes Interesse, fiir die ein 
Ersetzen des Profils durch eine Wirbel- und Quell-Senken-Belegung auf der Profilsehne recht 
gewaltsam erscheint. Andererseits ist aber auch schon hei wenig gekriimmten Profilen eine 
direkte Ubernahme der Riegelsschen Theorie nicht mehr méoglich, sobald die Schaufelabstande 
klein sind, da dann die gegenseitige Beeinflussung der Profile zu grof wird, wie in 5 d) gezeigt 
wird. 

Im folgenden werden nur Gitter aus wenig gekriimmten Skelettprofilen (Dicke Null) unter- 
sucht und nur im Anhang 5 c) und 5 d) wird noch kurz auf Profiltropfen eingegangen. Fir die 
Berechnung stark gekriimmter, 


dicht stehender Schaufeln kann die ee a Ce 
nachstehend entwickelte Theorie in Sept 
jedoch immerhin als erste Naherung Prof im i 
fiir ein Iterationsverfahren auf der offenenen = SS 
Basis von Gleichung (2) dienen. METS SB Rink ily 
d) Entsprechend, wie es in I — 
S. 116ff. geschehen ist, laBt sich co tee 


auch hier als einfachstes System von 
Schaufelgittern das ,,alternierende‘ 
Gitter einfiihren. Die Gleichung (1) 
ersetzt sich dann unter den obigen 
Vereinfachungen durch 
u(x) — tv(x) eye 
b Prot int); LLL 


eae 7(S) Luftstrah/ 
2 re dé, (la) zwischen Pe TE ire ee 


parallelen 
: MOT TTT: 
und Gleichung (2) modifiziert sich i 
entsprechend. eed 
Die Hauptanwendung des alter- Skelettprotil symmetrisches Profi] 


nierenden Gitters liegt im senkrech- 
ten Fall (h, = 0). In Abb. 3 sind 
die beiden Sonderfalle des symmetrischen und des Skelettprofils eingezeichnet, je einmal in 
gewohnlicher und einmal in alternierender Gitteranordnung. Man kann diese beiden Falle, 
wie skizziert, deuten als Profile im ebenen Kanal zwischen festen Wanden einerseits und 
(wenigstens in linearer Naherung) bei freier Strahlobrflaiche andererseits. In diesen Fallen 
stellt also die Gitteranordnung der Profile nichts anderes dar als einen Kunstgriff zur Er- 
zeugung des gewiinschten Strémungszustandes. 


Abb. 3. Anwendungsbeispiele fiir gewéhnliches und alternierendes Gitter. 


3. Der Falln = 1. a) Die erste Hauptaufgahe. Der Hauptpunkt dieser Aufgabe 
ist die Berechnung des Integrals 
b 
k 
SGN) etek (x — 6) dé 


a 


bei bekannter Zirkulation y(£). Die Behandlung dieses Problems wird durch die folgenden 
Entwicklungsformeln sehr erleichtert. Fiir gewohnliche Gitter gilt 


1 


S g sink Shain (E— a) sin? ea (E- ze k : - es "eos (:- a (x—&) d& 


2) 
= b a+b | | (5a) 
a cheatin iy x slim ne apa a Ay mi Sd Ya) Ba 
== SH g27 sin Ai 5 sin (= 9 eas et ip aad aber! ( 
| : >) 
Ie I 3et 8) 
fiir m = nae mit den Konstanten gg = —1,q, = 1/2 und 22 =—Q7G-@ 12d 


firA >1 eel 8. 110). Dabei sind die beiden Arten von Funktionen, die in der geschweiften 
8 
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; ' : 1 
Klammer getrennt werden, je symmetrisch oder antimetrisch zum Mittelpunkt 7 (a +b) der 
Strecke: <a, b>. 


Entsprechend findet man fiir alternierende Gitter 


b aids (ee a+b } 
“f /sin k (b —€) sin k(é — a) sin? ™—1k (: —— ‘ i ; i ae aneecs 


fed 


a+b 


: ab ae 5 )| 
Beer) ean 
- te Jeos ke x — 


ag oageand (Ug ea | 
fiir m hes 

Manchmal noch zweckmiafiger sind die folgenden Entwicklungsformeln fiir gewéhnliche 
Gitter: 


b sin 2my | 
k 
ee pel! o : . ctg k (x — &) dé = 
¢ te —sintk-—* cos’y sin (2m — 1) yp | 
cos2m@p 
(6a) 
W — sin? k * cos? gy cos (2m — 1) | 
(meer, 2 
mit sink/x — * qa sin k 2 = cos ~, entsprechend fiir € und y, sowie 0 y, yp <7; 
sowie fiir alternierende Gitter: 
si lgmeans * cost y sin 2 my | dé 
sin k (x — &) 
sin (2m — 1)yp | 
_| pi sine & =" cost y cos 2 mo | (6b) 


{ 
cos (2m — 1) | 


| 
(tigen Fo, 


b) Die zweite Hauptaufgabe. Die folgenden Ableitungen gelten nur fir die 
beiden Sonderfalle des senkrechten (hj = 0) und des waagrechten (h, = 0) Gitters, da die 
strenge Lésung der Fredholmschen Integralgleichung (2a) bei bekannter Funktion (x) bis 
heute nur fiir diese beiden Gitteranordnungen gelungen ist. Der Grund dafir liegt darin, daB 
nur fiir diese beiden Falle die induzierten Zusatzgeschwindigkeiten u + iv auf der Profilsehne 
<a,b> reine Abwartsgeschwindigkeiten werden, daB also gilt 


Se ae ctg k(x — dele 
(x) = Ref Ge v(@) etg k(x —8) :} =0 


a 


fir alle Zirkulationen y(€). Damit geht (2a) iber in 
é b 

; df (x k 

= sin ne cous ale : ae VIE) ots Fil eee (2b) 

wenn mit f(x) die Ordinate des betrachteten Skelettprofils bezeichnet wird. Wegen der 

Linearitat von (2b) 1aBt sich die zugehérige Lésung y(&) aufspalten in y(&) = y,(&) + V9(&). 


wenn y, die Lésung von (2b) fiir f’(x) = 0, y) entsprechend fiir ~« = 0 bedeutet. Fiir kleine 
Winkel « hangt der zweite Term fast nicht von « ab, so daB die Zirkulationsverteilung y sich in 
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linearer Naherung also aus einem zu sina ~ a proportionalen und einem anstellwinkelunab- 
hangigen Glied zusammensetzt. Die Auflésung von (2b) ist! unter sehr allgemeinen Voraus- 
setzungen méglich und die Lésung ist eindeutig bestimmt, wenn man noch die AbfluBbedingung 


y(b) = 0 fordert. 
c) f’(x) = 0 (Plattengitter). Man findet * aus (2b) die Zirkulationsverteilung 


2u. sina 7/sin k (b — é) 
Val) = pine Yee eo. (7) 
2 


cos k 


Nach der Formel von Kutta-Joukowski ist der zugehérige Auftrieb einer einzelnen Schaufel aus 


dem Gitterverband (9 = Luftdichte) 


b b—a 
tg k 
ae is fr df =2nsina fut, (b— a) — (8) 
een 
a DA 
Dieser Ausdruck unterscheidet sich nur oy 
durch den Faktor 
a 
tek oe waagrechtes Gitter 
hy i an ee ip 
bY Plame 
b—a Se). ae CT Sa (Sy ¥0) 


k 


von dem entsprechenden fiir ein Einzel- 
profil (Abb. 4). Das Moment um den 
Nullpunkt wird (vgl, I S. 119) 


b 
a+b 


M.. = 0x fro f di = Ay be 


a 


| 
| 
| 
| 
| eae 
| 
| 
| 


+ 2 ctg k 


b—a b—a 
5 “log cos b> | (9) 


und also die Lage des Auftriebsmittel- 
punktes fiir y,, der mit dem Neutral- 
punkt xy identisch ist 


Mx a+b 1 


| 


Miia. ee D senkrechtes bitter 
b—a b—a Abb. 4. Beeinflussung des Auftriebsanstiegs durch Gitteranordnung, 
xX ctg k 5} log cosk ach (10) Korrekturfaktor z. 
In Abb. 5 ist xy fir 6 — a = 1 in Abhangigkeit vom Teilungsverhaltnis t/l = |h|/(b — a) auf- 
getragen. 


Die Gleichuagen (7) bis (10) gelten fiir f’(x) = 0, d. h. fiir ein Plattengitter. Da in diesem 
Falle (2) und (2b) tibereinstimmen, sind die angegebenen Formeln nicht nur Naherungen als 
Lésungen einer linearen Theorie, sondern geltene x ak t. Sie stehen also in Ubereinstimmung 
mit den Ergebnissen der konformen Abbildung *. 


ad) x = 0. Die Lésung der Integralgleichung (2b) wird hier zu 


b ey 
We ks ae sin k (b — €) df (x) s. : _b—a sin k (% — a) 
vol) = ; V5 (E—a are ete Bee e os | Vee ge ace) 


eA sink a ) 


s 


Dazu gehért der Nullauftrieb eines Einzelprofils 


b b 


’ 2 our [ df(x)1/sin k(x — a) 2 
ae 


cos k 
a a 
1K, Nickel, Math. Z. 58 (1953), S. 49. 
2 Vol. fiir (7) und (8) ae Die hydrodynamischen Grundlagen des Fluges, S. 93—97. Braun- 


schweig 1917. 
8* 
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und das Moment um den Nullpunkt beim Anstellwinkel a = 0 


J 2 cos k(x —2 5) + |sin k (b — x) sin k (x — a) 
= EVE de = 20 us | eal 
My= ous [rlé) eds =2 out [GO| tog Fe: 
a a ~% x 
a+b sin k (x — a) 13 
baa Vimeo—a | 13) 
2 cos k ous 
t 
d 
senkrechtes Gitter a waagrechtes Gitter 


Abb. 5. Lage des Neutralpunktes. 


4, Der Fall n > 1. Die n Profilgitter aus Skelettlinien mégen alle dieselbe Gitterkonstante 
k =a/h =a/(h, +h.) besitzen und untereinander naherungsweise parallel verlaufen. Man 


° ersetzt dieses System durch eine 
S ‘ie o  Wirbelbelegung auf den  Profil- 
cari ° Sa sehnen, die o. B. d. A. alle parallel 
a dpe —= zur reellen Achse angenommen sind. 
ol Die Sehne eines beliebigen Profils 
= sae yB ch 4 des y-ten Gitters mége sich von a, 
4 é cage eR eo nach b, (mit Im {a,} = Im {b,\) 
CB OSEO ‘i erstrecken und sei so bestimmt, daB 
So ee ON i NS sie das Profilskelett nach Lage und _ 
& . — h 2 — {/ Lange méglichst gut annahert (vgl. 
° ir se e Abb. 6). Damit ersetzt sich die Be- 


as ee " stimmungsgleichung (2a) durch das 
Abb. 6. System aus n Profilgittern (ausgezogen). Ersatzsystem (gestrichelt) : . i 
Wirbelbelegung y auf den ,,Profilschnen“. te Gleichungssystem 


by 
Uo Sin |x — ¥,(%,,)] = cos x,,(z,) Im ' _ >» ¢ y(C,) ete kb (z, —C,) dé; ae 


er ese | 
a sin %,(2) Re >¢ VAC) ctg k (2, 5 Gs dc, ; : (u vs 1, fe n) (2c) 
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Dabei und im folgenden bezieht sich der Index fs immer auf den y-ten Fliigel, so daB z. B. 
?,(2,) die Zirkulation an der Stelle z, ist, die demnach zwischen a, und 6, liegt. Der schwie- 
rigste Teil der ersten Hauptaufgabe besteht wieder in der Berechnung der Integralsumme 


m by 
a VAG) ctg k (2, a ¢,) dc, 
bei bekannten Zirkulationen y,. Wie in IS. 115/116 ausgefithrt ist, lassen sich fiir die y, Entwick- 
lungsformeln angeben, die analog zu den Gleichungen (5a) und (5b) gebaut sind. Nur erfordern 
sie naturgema§ erheblich mehr Schreibarbeit, und da sich auch die Entwicklungskoeffizienten 
i. a. nicht geschlossen angeben lassen, soll hier auf eine Wiedergabe verzichtet und sofort die 
zweite Hauptaufgabe behandelt werden. 

Wie im Falle n = 1 mu8 man auch hier k reell oder rein imaginar (Gitterfronten waagrecht 
bzw. senkrecht) vorschreiben. Dariiber hinaus miissen aus der Menge aller solcher Gitter- 
systeme noch diejenigen ausgesucht werden, bei denen die induzierten Geschwindigkeiten am 
Orte der Profilsehnen reine Abwindgeschwindigkeiten sind, da nur hierfiir die strenge Lisung 
des Integralgleichungssystems (2c) bekannt ist. Das ist z. B. der Fall, wenn alle Profilkonturen 
hintereinander liegen, doch ist auch eine Staffelung um h/2 erlaubt (man vgl. I S. 114). Unter 
diesen Einschrankungen ist dann 


| ik x ; | 
Re = >¢ y(cy) etg k (z,, —C,) =) = 


und (2c) geht iiber in 


GiGi k \ 
th [sin a — cos a TH) = TD lls) 18 bey 0) (24) 


Wie im Falle n =1 1a8t sich auch hier wieder fir die Zirkulationen die Aufteilung 
¥(C,) = Yavr(6,) + Yo,(C,) vornehmen, wobei das erste Glied y,, fir f, = 0 und 9, fiir « = 0 er- 
halten wird. 


b) f.(z,,) = 0 (Systeme aus Plattengittern). Das Integralgleichungssystem (2d) hat, falls 
man die Abstrémbedingung y,(b,) = 0 (vy =1,...,m) voraussetzt, fir fi = 0 die eindeutige 
Auflésung * 


2 ux» sina Plain bbe Cy) oe 
Vy Cy) = ee IT ak ac (Gh ule tae, (14) 
cosk S) eae x=1 
e=! 


Dazu gehért der Gesamtauftrieb (Summe der Auftriebe von je einer Schaufel aus jedem Gitter 
des Systems) 


bo ao 
i: : ghd 
: ; iy Cues = : 
A, = 0 Ux > [rats Nise = 2 7 SiN 0 -F- Uso (Oe aes, nf aw ae (15) 
v=l1 [3 k See Se st 
ay O} 
gel) 


Das Gesamt-Langsmoment 
m by 
M, SSD es Py i Vane) o dc, 
y=la, 


148t sich fir n >1 nicht in geschlossener Form durch elementare Funktionen darstellen. Aus- 
driicklich sei auch hier wieder bemerkt, daB diese Ergebnisse fiir die genannten Systeme aus 
Plattengittern exakte Giiltigkeit besitzen. 


1 Siehe FuBnote 1 von S. 107. 


110 Nickel: Uber spezielle Systeme von Tragfliigelgittern II. Ingenieur-Archiv 


c) « = 0. Die Auflésung von (2d) liefert fir « = 0 


2 ihe sin k (b x — Cy) df u(zu) x 
Y'0 Ge )= i ee k ( C =a S$ dx 


Usa 


i sin k ( cee 
etek (25 —=¢,) tek > : _ |/ ie k (ba — Zu) dz, (16) 


o=1 


x 


fuie (ssl, 6 bo pide 
Dazu sans der Gesamtauftrieb 


b 
n Nd ya oe 
= ra 2 9 us * dfu(Zz) é sin k (az — 2p) 
rs on 5 | Vowly) @y = — ge i a Ao dT he dz,. (17) 
4 bo ao = b A=1 
ce cosk > 9. 4 (Se 
g=l 


Auch hier laBt sich das Gesamtmoment nicht geschlossen angeben. 


d) Folgerungen. Beziiglich der folgenden Aussagen vergleiche man mit 1, wo die 
entsprechenden Folgerungen fiir ein System aus n Einzelfliigeln gezogen werden. Es sei hier 
noch einmal daran erinnert, daB an dieser Stelle nur waagrechte und senkrechte Profilgitter 
betrachtet werden, deren Schaufeln hintereinander oder (evtl. nur teilweise) ,,auf Liicke“ 
stehen. Zur Vereinfachung der Schreibweise werde noch sin « durch « und cos a durch | ersetzt, 
man betrachtet also nur kleine Anstellwinkel a. 

1, Auftriebsverteilung ? c,(z) = 2 y(z)/u., Gesamtauftrieb A und Moment M sind linear 
im Anstellwinkel a, es gilt also 


: dca(zy) 
Ca(2,) = €a,(2,) ai a Ah ? 


2. Es gilt allgemein 
dea(zy) 4 + sink (bz — zy) 
da n —\ UE heen morse, i 


3. Es gilt allgemein 


eb 
‘ tek > ae 
“t —2n et =27-x 
da * by ce 
ae 5 


vgl. Abb. 4, wenn man dort | = 6 — a durch 5) b, — ag ersetzt. Der Auftriebsanstieg aller 


e=t 
Gittersysteme mit gleicher Summe der Fliigeltiefen ist also gleich. 
4. Es gibt genau einen Profilneutralpunkt N (aerodynamic center), um den das Moment M 
anstellwinkelunabhangig wird, d. h.: bezogen auf N ist dc,,/da = 0. Fir n = 1 ist die Abszisse 
zy des Neutralpunkts 


iy ae ie 


1 Siehe FuBnote 2 von S. 102. 
* Die Verwendung von ca(z) zur Kennzeichnung der értlichen Auftriebsdichte und von cq als Beiwert 
des Gesamt auitriebs wird kaum zu Mi®verstandnissen AnlaB geben, In der Theorie der tragenden 


Linie wird ebenfalls seit langem schon ca als Beiwert des Gesamtauftriebs neben Cq(y) als 6rtlichem Auf- 
triebsbeiwert benutzt. 
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(Abb. 5), fiir n > 1 gibt es dafiir keinen geschlossenen Ausdruck aus elementaren Funktionen, 
doch gilt dann wenigstens noch die Abschatzung Re {a,} < Re {zn} =< Re{b,}, d.h., der 


Neutralpunkt liegt immer innerhalb des Gittersystems. Andere, genauere Abschatzungen 
lassen sich leicht von Fall zu Fall gewinnen. 


5. Fir f,(z,) > f,(b,) ist ebenfalls c,, > 0. Dabei tritt das Gleichheitszeichen bei c, nur 
fiir f,(z,) = f,(b,) ein. Oder: verschiebt man alle Profile senkrecht so, da JT A0;) = ist fiir 
y =1,...,m und legt man in Anblasrichtung! durch den hintersten Punkt b, des Fliigel- 
systems eine Gerade, so kann diese bei Gesamtauftrieb Null nie ausschlieBlich auf einer 
Seite des Systems liegen. 

6.n —1. Wenn sich die Tangente der Fliigelkontur vom vordersten bis zum hintersten 
Punkt immer gleichsinnig weiterdreht, so ist Cm, == 9 (Cm, bezogen auf den Neutralpunkt). 
Oder: Druckpunktfestigkeit ohne Wendepunkt im Profil ist nicht méglich. — Ob dieselbe 
Aussage auch noch fiir n > 1 gilt, konnte der Verfasser nicht nachweisen. 

7, Sonderfall der Streckenprofile mit demselben Anstellwinkel: Fiir diese ist ¢,,(2,) = Ca, 
= Cm, = 9. Der Gesamtauftrieb ist bei allen Anstellwinkeln derselbe, wie wenn die einzelnen 
Fligel zu eine m einzigen Fligel zusammengeriickt waren; der Druckmittelpunkt liegt im 
Neutralpunkt. Bei senkrechten Plattengittern (h, = 0) besitzt die Strémung weit hinter 
(+ oo) bzw. weit vor (— oo) dem Gitter die Komponenten 


[U +iV],~.= a % 4 eal + Yq = =| ; (18) 


Ein solches Gitter lenkt die ankommende Strémung also nie vollstandig in die Waagrechte, die 
Richtung der Platten, um. 

Die Beweise zu den ausgesprochenen Behauptungen ergeben sich entweder aus dem Vor- 
angegangenen, oder kénnen genau wie in meiner friheren Arbeit ? gefiihrt werden. Sie sollen 
deshalb hier nicht gebracht werden. 


5. Anwendungsbeispiele. a) Spezielles Plattengitter aus zwei Pro- 
fiteu,. Bs sein — 2,0, —0.h, — 1 a, =2---1,b,=5 +1,h th, = 27% (Abb. 7). Damit 
gilt nach (14) allgemein 


Sin Z (1 — x) Cof F (5 — x) 
Gin Z x Coj Z (2 — x 
_ 4sina | ir) VE ( ) 


A oer ere | 1 /Goi Fle — 1) Sin FG — ») 
LV CofZ x Sin Z (x — 2) 


ru ONS el 


l\ 
2 
I 

| 
/\ 
wn 


fir 2< 


Siehe Abb. 7, wo zum Vergleich noch c,(x) fiir den einfachen Doppelfligel mit dem Strecken- 
profil0 <x <1,2< % <5 eingezeichnet ist. Die Auftriebsdichte, die fir das Doppelgitter 
schon an der Profilvorderkante erheblich unter der entsprechenden des Doppelfliigels liegt, 
wird am Ende des Hinterprofils verschwindend klein, weil dort die Strémung durch das Vor- 
dergitter schon stark umgelenkt ist. 

Nach 4 d) 3) ist weiter c, = 2 Tg a sina = 1,9925 sina. 

Die Strémungskomponenten im Unendlichen vor und hinter dem Gitter sind nach (18) 


[U +iV]z. = Ux jeos ao + wae 
Das bedeutet, da® dieses keineswegs besonders engstehende Doppelgitter die V-Komponente 
bis auf den kleinen Rest von nur 1,9°/og (!) in die Waagrechte umlenkt. 

b) Als Beis piel fiir die Anwendung der Formeln (6a) und (6b) werde fiir das senkrechte 
Gitter (h, = 0) die erste Hauptaufgabe fiir Skelettlinien gelést. Das Ziel sei dabei die Gewin- 
nung von Summenformeln, wie sie von Truckenbrodt ® fiir das Einzelprofil aufgestellt wurden. 
Man geht dazu aus von Gleichung (2b), worin man offenbar ~ = 0 setzen darf, da aus (7) die 
zusatzliche Lésung von (2b) fiir « + 0 schon bekannt ist. Spaltet man in (2b) (€) und f(x) in 


| isin 5 


1 D—, h. in der Richtung, die durch u,, (cos « + isin a) festgelegt ist, 
2 Siehe FuBnote 2 von S. 102. 
3 Siehe FuBnote 1 von S, 104. 
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ihre symmetrischen (yt, f+) und antimetrischen (y~, f~) Bestandteile auf und setzt noch 


wo. 6 Gee UN ha RH eee 19 
Gin * (« 5 Sin 5 ag or ke (O55 @y Sx0t)s5 (19) 
so geht (2b) mit der Schreibweise y(x) = y(9), f(x) =f (@), he = |h| = t, 6 — a =liberin 
df(p) __ l : oy +(p) r(p) sin p dy r y (yp) sin y dy 20 
“de.” ~ Baus sno g cos y — cos p 9) Cans omteere eo 
0 
(man vel. dazuI S$. 111). Diese Formel stimmt, bis auf die Glieder 
ae 
Sin — — 
2 1 
r(y) =——- (21) 


é 1 
Ay yi + Sin? a cos? 


iiberein mit der entsprechenden Forme] fiir den Einzelfliigel und geht fiir t > oo mit r(p) > 1 
auch in diese ttber. Durch die Koordinatentransformation (19) hat man also die gestellte Auf- 
gabe, bis auf die Faktoren r(y) und r(w) auf eine 
bekannte Gleichung zurickgefithrt und -kann 
die dafiir entwickelten Methoden mitbenutzen. 

Doppel- ————_ _ Etwas iiberraschend ist dabei am Anfang die 
ee aay Aufspaltung von y in einen symmetrischen und 

Pe einen antimetrischen Teil, da man eine solche 
Aufspaltung beim Einzelprofil nicht kennt. Sie 
rubrt daher, da sich diese beiden Teile bei der 
° Transformation (19) verschieden verhalten und 
also auch getrennt behandelt werden miissen. 

Mit der Truckenbrodtschen Bezeichnungs- 
weise y(x) = f(x), g(x) = + y(x)/2 u. und den 
Doppel- dort definierten Konstanten c,,, kommt man 


ae schlieBlich auf die Summenformel 1 


O° 
1} 
° 


e) 


ooo 


N—1 


; N—1 
Vrs = Oey Nee Tm = ey Cmn &m- : (22) 


m=1 m1 
Dabei bedeutet 7, =—9(¢,) mito, = ay und 
entsprechend r, =r(@n), gn =8(Pn), die Ab- 
leitung nach @ bei y ist durch einen Punkt 


gekennzeichnet. Die Aufspaltung g,, = gi + gh 
geschieht nach der Vorschrift 


oS AS 5 ~ 
Sm aa &N—m> Sm ee tee §&N—m ° 


Abb. 7. Auftriebsverteilungen cq(x) eines senkrechten F e 
Doppelgitters aus Streckenprofilen und des entsprechenden Der Zusammenhang mit x ist durch 


Doppelprofils. 


¥n = ¥(%n). Bn = 8(%n) 


t vA eats mn 
+ — Ur Sin | Sin 7g O08 N 


a+b 
ats 


mit 


(23) 


gegeben. In Tabelle 1 findet man x, und r, fir N = 12 und die drei Teilungsverhaltnisse 
Hh. 05 0,75-und 0:5: 

; Fir symmetrische Zirkulationsverteilungen, also g = gt ]aBt sich noch die Integration von 
y geschlossen ausfiihren und liefert mit den Truckenbrodtschen Konstanten d 


LP yt 
yn =i Sidue ey Tm: (24) 
m=1 


Bei antimetrischen Zirkulationsverteilungen kann man y etwa durch graphische Integration 
der Funktionswerte aus Gleichung (22) erhalten. 


t Der Einfachheit balber nur fiir durchweg endliche Geschwindigkeitsverteilungen angegeben. 
2 Die Uberstreichung von y deutet daraui hin, da ¥ nicht exakt, sondern durch eine Summenformel 
berechnet wird. 
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Die (20) entsprechende Gleichung fiir alternierende Gitter lautet 


4 A ~ 
df (p) ; y*(w) sin p dp os y) sin p di 
Et y (py) r(p) sin yp dy 
sas sin r a CaN ee y 
dy 4 Us P (9) ¢ cos Pp — Cos p ¢ cos P — cos p ‘ (25) 
0 0 | 
Daraus folgen wieder die Summenformeln 
. N—1 N—1 
Ae > Cnn en > ey Cmn Sm Tn (26) 
m=] T= 
und 
sod N—1 
Os = ps a Em Min > (27) 


m=l1" 


wobei die integrierte Gleichung (27) diesmal nur fiir antimetrische Zirkulationsverteilungen 
giiltig ist. 


o nach Summentorme/ 


Einzelflugel 


~ 


Abb. 8. y(x) (rechte Bildhalfte) und y’(x) (links) fiir konstante Zirkulationsverteilung bei verschiedenen Teilungen, 
y(x) ist symmetrisch, y’(x) antimetrisch in die andere Bildhalfte fortzusetzen. 


Als Zahlenbeispiel sei der Fall g, = gi = 1 gewahlt. Es ergeben sich fir a = 0, b =1 
(Profillange), h, = t (Gitterteilung) 


yf (eas = log — ae fir gewohnliche 
Gin = « 
t 
und 
Tq (1 — x) 
F il ot a . 

y’(x) = — log —— fir alternierende 

2 ONE 


Gitter. Fiir t > oo streben diese beiden Ausdriicke gegen die bekannte Funktion 
f 1 l— «x 
ye) Slog 
In Abb. 8 (rechte Halfte) ist y(x) aufgetragen. Die Ubereinstimmung zwischen exakten Werten 
(Integration numerisch durchgefiihrt) und Summenformel (NV = 12) ist zufriedenstellend, nur 
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im Falle t/1 = 0,5 streuen die Werte in Profilmitte etwas starker. Die Ubereinstimmung bei y’ 
(ebenfalls Abb. 8, linke Hilfte) mit der Summenformel ist fiir gewohnliche Gitter nicht sehr 
gut. Bei alternierenden Gittern erhalt man nach (25) die Kurve fiir y’ exakt durch Verzerrung 
in x-Richtung nach (23), wie die Waagerechten durch die Punkte y, zeigen, da in dem spe- 
ziellen Falle g = 1 die Funktion g(y) mit g(x) itbereinstimmt. Die Ungenauigkeiten der Sum- 
menformel fiir Einzelprofile (¢ —-oo) werden hierbei durch die Gitteranordnung also nicht 
verstarkt. 

Ein gewisser Nachteil der beschriebenen Methode soll nicht unerwahnt bleiben: schon bei 
Einzelprofilen (t/l oo) stért fiir manche Zwecke die Tatsache, da die Rechenpunkte x, 
nicht diquidistant iiber die Profiltiefe verteilt sind, sondern sich am Vorder- und Hinterende 
haufen (Abb. 8 und Tabelle 1). Diese Eigenart, die auf die benutzte Gaufsche Quadratur 


abelle ie mind e7peNe— 2) 


a a ee a ee Ce ae jn MELO 11 


0,9830 | 0,9330 | 0.8536 | 0,7500 | 0,6294 | 0,5000 | 0,3706 | 0,2500 | 0,1464 | 0,0670 | 0,0170 


0,3282 | 0,3657 | 0,4468 | 0,6273 | 1,1662 | 3,6761 | 1,1662 | 0,6273 | 0,4468 | 0,3657 | 0,3282 


3 


x 
sais aS 1,0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 | 1,0000 
fl = 1,00 {%, | 0.9899 | 0,9585 | 0,9021 | 0,8132 | 0,6799 | 0,5000 | 0,3201 | 0,1868 | 0,0979 | 0,0415 | 0,0101 
: \ rn | 0,6011 | 0,6570 | 0,7671 | 0,9610 | 1.2587 | 1,4651 | 1,2587 | 0,9610 | 0,7671 | 0,6570 | 0,6011 
sae {xn ] 0.9920 | 0,9668 | 0,9208 | 0,8446 | 0,7163 | 0,5000 | 0,2837 | 0,1554 | 0,0792 | 0,0332 | 0,0080 
aN \rn | 0.4785 | 0,5297 | 0,6366 | 0,8541 | 1,3266 | 1,9093 | 1,3266 | 0,8541 | 0,6366 | 0,5297 | 0,4785 
eee {x | 0,9945 | 0,9772 | 0,9451 | 0,8906 | 0,7889 | 0,5000 | 0,2111 | 0,1094 | 0,0549 | 0,0228 | 0,0055 

o 


zurickgeht, wird durch die Transformation (19) noch verstaérkt. Fir kleine Teilungen sind 
deshalb fast alle Punkte x, an den Profilenden konzentriert, wahrend in Profilmitte nur der 
Punkt x, tibrig bleibt. Z. B. liegen beim Teilungsverhaltnis t/1 = 0,5 auf der ganzen Strecke 
zwischen 11% und 89% der Profiltiefe / nur die drei Punkte x, = 0,2111, x, = 0,5000 und 
x, = 0,7889 (man vergleiche Abb. 8 und Tabelle 1). Dieses Verhalten ist ein Schénheits- 
fehler, der als Ausgleich fiir den einfachen Bau der angegebenen Formeln in Kauf genommen 
werden muB. 

c) Als weitere Anwendung fiir die Formeln (6a) soll noch die Theorie der Profiltropfen, 
wie sie von Riegels 1 gegeben wurde, vom Einzelprofil auf ein senkrechtes (h, = 0) Profilgitter 
tibertragen werden”. Zur Vereinfachung werde nur das symmetrische Profil ohne Anstellwinkel 
behandelt, ebenso mége der Fall des alternierenden Gitters unerwahnt bleiben. 

Ubertragt man die Riegelsschen Resultate sinngemaB, so findet man die Grundformel * 

b 


oder nach partieller Integration 
W i) kd 
i 


P= presal! te eG 10) ote he — 2) a]. (28) 


Dabei ist y(x) die Profilkontur und u. =V gesetzt, mit W bezeichnet Riegels die Strémungs- 
geschwindigkeit an der Kérperkontur. Man kommt auf (28) also durch eine Quell-Senken- 
Belegung der Strecke <a,b> mit der Starke dy(x)/dx. Schon hier sei bemerkt, daB diese 
Annahme, die sich bei Einzelprofilen gut bewahrt hat, bei Gitterprofilen zu Schwierigkeiten 
fiihrt, sobald die Gitterschaufeln eng stehen. 

Als Beispiel werde 


Gin = (b — x) (29) 


1 Siehe FuBnote I von S. 104 

g Zusatz bei der Korrektur: Dieselbe Aufgabe wird in einer inzwischen erschienenen Arbeit von 
L. Speidel, Ing.-Arch. 22 (1954), S. 295 behandelt. 

* Man vgl. auch Byrd una Huggins, Fufinote 3 von §, 102. 


XXIII. Band 1955 Nickel: Uber spezielle Systeme von Tragfliigelgittern IT. 1S 


gewahlt. Man findet damit fiir W/V nach (28) den geschlossenen Ausdruck 


V° + Gin? 2 us 


WwW 1 ME Dat | ae | 
— = =-5- hates =a a > + e oe 2 eee 14 . Ca b 
Pee rena sean eral ets (a7 
t a : 
: Te ae Ai). 
Cae it otal Bee et 
e? Go] = — Gin 2 AG ak (30) 


# Abb. 9 ist diese Profilkontur mit der zugehérigen Geschwindigkeitsverteilung eingezeichnet 
ur 


t h, 
— = z as I's (pe == 2 und E> eee — 0,018 882 ’ 


l b—a 


Um wieder zu einer Summenformel zu gelangen, spaltet man y(x) in einen symmetrischen 
(yt) und einen antimetrischen Bestandteil (y~) auf. Die Setzung (19) fiihrt dann (28) iiber in 


© nach Summentorme/ 
o mit Abspalting 
(volkommene Ubereinstimmung) 


14 
Tol df y*(y)r(y) sin y dy 
cos Y — Cos y 17 
0 
A 10 
1d are sin y dy Pike 


cos m — cos wy 


cos p — cos yp) 


Beim Ubergang zu einer Fourierreihe kommt Ee Ee ea 


man fiir die beiden ersten Integralglieder auf 
Summen mit den Riegelsschen Konstanten - Abb. 9. Profilkontur ae en ae ee nach (29) 
Amn. Anstatt jedoch die gegebene Profil- 
kontur y(x) in einen symmetrischen und einen antimetrischen Bestandteil aufzuspalten, ist es 
zweckmaBiger, diese Aufspaltung bei den A,,,, vorzunehmen durch 4,,, = At,, + Ann mit 
Pes p= Amn An—mn Und 2.45, = Amn — An—mn- Damit,.findet: man. schlieBlich die 


Summenformel 


0 
Pine f yi(y) r(y) sin yp _ (1) 
0 


W 1 | LS 44 sus 
V —— ae ae la,,| Th a ae poe AS n 2 VYmlm ae 24 ore 2 Ym ss 
12. Toll On 7 m=1 We 
2 N—1 
+ Pp Tn 25 AA,, n 2 Vm Tm |-+ (32) 
nm 


Die Punkte n sind wieder dieselben wie in b), ebenso ist r, von dort ibernommen. ¥ bedeutet 
wieder Ableitung nach gy, die Konstanten a, und c, sind dieselben wie bei Riegels. Die Koeffi- 


zienten At, A>, und AA,,, sind fiir N = 12 in den Tabellen 2 bis 4 zusammengestellt, 


mn? 


AA, , ist durch 


L200 Tn 
Sh 5 
N 1 J 
= 2 ae et ot (Mees m-+-n eA et ee ee : 
AAmn Hf 32 N be ae | a (m -+-n) ze _ a@(m—n) (e3) 
sin a sin ie 


erklart. 
Als Anwendung dieser Summenformel wurde die Geschwindigkeitsverteilung W/V fiir das 


- Beispiel der Profilkontur (29) bestimmt, sowie beim Joukowskiprofil 0015 die drei Teilungen 


t/l = 1,0; 0,75 und 0,5 mit t/l = co (Kinzelprofil) verglichen. Das Ergebnis zeigt Abb. 10. 


Wie man sieht, streuen die Punkte der Summenformel fiir die kleinste Teilung t/l = 0,5 in 


Profilmitte unzulassig stark. Ebenfalls eine groBe Abweichung findet man in Abb. 9. Der 
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Tabelle 2. At, N=12. 


es! 2 3 4 5 6 OES ees 9 10 11 
n=0 0 0 0 0 0 0 0 0 0 0 0 

1 | + 1,50000] —0,54197/ 0  |—0,04737/ 0 = 0,02312| 0 — 0,04737 0 — 0,54197] +- 1,50000 

2 | —0,54197| + 1,50000/—0,58933| 0 |—0,07049 0  |—0,07049} 0 — 0,58933|-+ 1,50000] — 0,54197 

3 0 — 0,58933| + 1,50000|—0,61245] 0 —0,11785} 0 — 0,61245| + 1,50000 | — 0,58933 0 

4 | —0,04737 0  |—0,61245] + 1,50000] — 0,65982 0 |—0,65982] + 1,50000| — 0,61245| 0 — 0,04737 

5 0 —0,07049} 0 ~~ | — 0,65982) + 1,50000] — 1,20178) + 1,50000| — 0,65982 0 — 0,07049 0 

6 | —0,02312 0 |—0,11785| 0  |—1,20178] + 3,00000]—1,20178] 0 — 0,11785 0 — 0,02312 

7 0 —0,07049] 0 —|—0,65982]+ 1,50000] — 1,20178]-+ 1,50000| — 0,65982 0 — 0,07049 0 

8 | —0,04737 0  |—0,61245|+ 1,50000]— 0,65982) 0 —|— 0,65982) + 1,50000! — 0,61245| 0 — 0,04737 

9 0 — 0,58933|-+ 1,50000]—0,61245| 0 —0,11785| 0 — 0,61245| + 1,50000|— 0,58933 0 

10 | —0,54197] + 1,50000|— 0,58933/ 0 |—0,07049] 0 |—0,07049 0 — 0,58933|+ 1,50000| — 0,54197 

11 | + 1,50000] —0,54197} 0 |—0,04737) 0 —0,02312; 0 — 0,04737 0 — 0,54197| + 1,50000 

12 0 0 0 0 0 0 0 0 0 0 0 

Tabelle 3. A,,,, N = 12. 

m 1 2 3 4 5 6 7 8 “) 10 11 
n=0 0 0 0 Ores lu. 20 0 0 0 0 0 0 

1 |-+ 1,50000|— 0,53874; 0 — 0,03867 0 0 0 + 0,03867 0 + 0,53874| — 1,50000 

2 |— 0,53874|+ 1,50000|— 0,57742 0 —0,03867| 0 + 0,03867 0 + 0,57742| — 1,50000| + 0,53874 

3 0 — 0,57742| + 1,50000| — 0,57742 0 0 0 +0 ,57742 | — 1,50000| + 0,57742} 0 

4 |—0,03867) 0 — 0,57742| + 1,50000] — 0,53874| 0 +- 0,53874] — 1,50000| + 0,57742| 0 | +4 0,03867 

5 0 —0,03867| 0 — 0,53874| + 1,50000| 0 — 1,50000| + 0,53874 0 + 003867; 0 

6 0 0 0 0 0 0 0 0 0 0 0 

7 0 + 0,03867} 0 + 0,53874| —1,50000} 0 + 1,50000| — 0,53874 0 —0,03867) 0 

8 1+ 0,03867| 0 + 0,57742| — 1,50000| + 0,53874| 0 — 0,53874 | + 1,50000 | — 0,57742 0 |—0,03867 

9 0 + 0,57742|— 1,50000| + 0,57742 0 0 0 — 0,57742| + 1,50000] — 0,57742| 0 

10 |+ 0,53874|— 1,50000|+ 0,57742 0 + 0,03867| 0 — 0,03867| 0 — 0,57742| + 1,50000|— 0,53874 

11 |—1,50000]+ 0,53874| 0 4. 0,03867 0" Maas 0 0 — 0,03867 0 — 0,53874) + 1,50000 

12 0 0 0 0 0 0 0 0 0 0 0 

Tabelle 4. 4A,,,, N = 12. 

m 1 2 3 4 Tee oe: rae 8 9 10 11 
n=0 0 0 0 Cn ewilne 0 0 0 0 0 0 

1 0 + 0,13565 0 + 0,06296 0 + 0,05322 0 + 0,06296 0 + 0,13565| 0 

2 |—0,10904] 0 ++ 0,21809 0 + 0,10904 0 + 0,10904; 0 + 0,21809] 0  |—0,10904 

3 0 — 0,14539 0 +. 0,25183 0 + 0,14539 0 + 0,25183| 0 |—0,14539/ 0 

4 |—0,01687| 0 — 0,12591 0 + 0,23496 0 + 0,23496 0 —0,12591| 0  |—0,01687 

5 0 — 0,00974 0 —0,06296} 0 + 0,19861 0 — 0,06296 0 |—0,00974| 0 

6 0 0 0 0 0 0 0 0 0 0 0 

q 0 — 0,00974 0 — 0,06296 0 | + 0,19861 0 —0,06296| 0  |—0,00974| 0 

8 |—0,01687; 0 — 0,12591 0 ++ 0,23496 0 + 0,23496 0 —0,12591] 0  |—0,01687 

9 0 — 0,14539 0 + 0,25183 0 + 0,14539 0 + 0,25183 0 |—0,14539} 0 

10 }—0,10904} 0 + 0,21809 0 =|+0,10904| 0 + 0,10904) 0 + 0,21809] 0 — 0,10904 

11 0 4 0,13565 0 | 4 0,06296 0 + 0,05322 0 + 0,06296 0 + 0.13565] 0 

12 0 0 0 0 ‘aaa ae 0 0 0 0 0 
~ Se ° : 5 - : a E a+b 
Grund dafir liegt darin, daB die Funktion r(p) nach Gleichung (21) fiir t/l +0 bei x = 4 
also pm = es uber alle Grenzen wachst und auch schon bei t/] = 0,5 dort eine starke Spitze hat 


SPS 


© nach Summentorme! 


a mit Abspaltung 


Abb. 10. Joukowskiprofil 0015 mit Geschwindigkeitsverteilung 
fiir vier verschiedene Teilungen. 


(vgl. Tabclle 1). Diese Spitze tibertragt sich 
auch auf die Funktion yt(g) r(p) und es ist 
verstandlich, daf sich solch eine Funktion 
nicht mehr gut durch wenige Glieder einer 
Fourierreihe approximieren la8t. Um diese 
Schwierigkeit zu umgehen, empfiehlt sich 
eine Abspaltung derart, daB man y+(q) durch 
eine symmetrische Funktion e+(p) annahert, 
fiir die sich das Integral in Formel (28) 
geschlossen berechnenlaBt!. Die Annaherung 
braucht dabei aus dem oben erwadhnten 


Grunde nur bei x 5 (a + b) gut zu sein. 


Damit kennt man W/V fiir diese Ersatzfunk- 
tion exakt (wenigstens im Rahmen der ge- 
schildertenTheorie) und braucht dieSummen- 


formel nur noch auf die Differenzfunktion 
y — et anzuwenden. 


* Kine solche Abspaltung 1aBt sich selbstverstandlich auch bei den Skelettprofilen des Abschnitts b) 


durchfiihren, wobei etwa die angegebene Lésung fiir g = 1 als Abspaltungsfunktion verwendet werden 
kann, 
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ee Punkte ,,mit Abspaltung‘ in Abb. 9 und 10 sind auf diese Weise gewonnen und zwar 
wurde 


Goj (x25) 


are COs 


et(x) = (34) 


2 


are COS 


mit der Dicke d = 2 y(2) gewahlt. Dazu gehért nach 1S.114 eine konstante Zusatzgeschwin- 
digkeit und man findet 


—|. (35) 


Mit diesem Ergebnis lassen sich auch die Geschwindigkeitsverteilungen von Abb. 10 fiir 
das Joukowskiprofil 0015 in ihrem Verlauf verstehen. Nach (31) wird namlich die antimetrische 
Geschwindigkeitsverteilung, die durch y~ erzeugt ist, durch die Gitteranordnung nur wenig 


Gebiet mit Widerspruch 


zur Kontinurtatsgleichung 


Abb. 11, Erklarung im Text. 


geandert. Weiter wird y* fir nicht zu kleine Gitterteilungen durch e* nach (34) sehr gut an- 
genahert und die dazu gehérige symmetrische Geschwindigkeitsverteilung ist nach (35) 


(bis auf den Faktor —— 
1+ y2 
Wert vergréBert. Dies ist aber gerade der Sachverhalt von Abb. 10, wo die drei W/V-Kurven 


fir t/] < coim wesentlichen durch Parallelverschiebung aus der Kurve t/l = oo fiir das Ein- 
zelprofil hervorgehen. 

d) Nun hat jedoch die in c) geschilderte Theorie fir Profiltropfen in Gitteranordnung noch 
einen bedenklichen Mangel. Als lineare Theorie in y(x) bestimmt sie eine Geschwindigkeit an 
der Kérperoberflache, die von dem unwesentlichen Faktor yl + y" abgesehen, nur proportio- 
nal zur Profildicke d wiachst. Dieses Verhalten steht aber in Widerspruch zur Wirklichkeit, 
da fir d/l— t/l, also bei verschwindender Breite der Schaufelkanale, diese Geschwindigkeit 
aus Kontinuitatsgriinden gegen Unendlich streben muff. Eine lineare Theorie kann also nur 
als Naherung fiir gro®e Teilungen t/] angesehen werden und muf bei kleinen Teilungen not- 


konstant und gegeniiber dem Einzelprofil um einen konstanten 
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wendigerweise versagen. Zum Beweis, daB dieses Versagen bei den itblichen Profildicken d/l 
schon bei recht groBen Teilungen t/l eintritt, werde die maximale Geschwindigkeit langs der 
Profilkontur mit der Kontinuitatsgleichung abgeschatzt. Zwar liegt der Anwendung der 
Kontinuitatsgleichung das Modell der eindimensionalen Strémung, der ,otromroéhre**, zu- 
grunde, da jedoch eine Potentialstrémung ihre Maximalgeschwindigkeit immer auf dem 


Rand annimmt, kann die héchste Umstrémungsgeschwindigkeit p max eines Profils niemals 
kleiner sein als der Wert, den man nach der Kontinuitatsgleichung berechnet. Es gilt also 
fiir jedes senkrechte Fliigelgitter der Teilung t aus symmetrischen Profilen mit der Dicke d 
und der Tiefe / die Ungleichung 


W ty 1 
er 


Setzt man (35) in diese Ungleichung ein, so ergibt sich als Grenzkurve 


TEm AU : 
wee oe arc sin —————__ (37) 


(vgl. Abb. 11). Im ganzen Gebiet links von dieser Kurve steht die vorstehend geschilderte 
Theorie im Widerspruch zu der grundlegenden Kontinuitatsgleichung, zumindest fiir das 
spezielle Gitter aus den Profilen nach Gleichung (34). .Sie kann also z. B. bei dem Dicken- 
verhaltnis d/l = 0,15 schon bei dem Teilungsverhaltnis t/] = 0,75 keine sicheren Resultate 
mehr liefern. Die Geschwindigkeitsverteilung des Joukowskiprofils 0015 fir t/] = 0,5 und 
t/t = 0,75 in Abb. 10 liegt schon jenseits dieser Giltigkeitsgrenze und besitzt damit nur noch 
Vergleichswert. 

Der Grund fiir dieses Versagen der erweiterten Riegelsschen Theorie ist offenbar darin 
zu suchen, daB fir kleine Teilungen eine zu y'(x) proportionale Quell-Senken-Verteilung die 
Wirklichkeit nicht mehr genau genug wiedergibt und da hier die gegenseitige Beeinflussung 
der Gitterprofile schon in der Quell-Senken-Verteilung beriicksichtigt werden miBte. 


(Eingegangen am 29. Marz 1954.) 


Anschrift des Verfassers: Dr. Karl Nickel, Instituto Aerotecnico, Cérdoba, Argentinien. 


= 
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Kin graphisches Verfahren zur Lésung von linearen Differential- 
gleichungen zweiter Ordnung. 


Von W. Vogel. 


Die Berechnung von Schwingungsvorgangen des Maschinenbaues und der Elektrotechnik 
erfordert haufig die Lésung von linearen Differentialgleichungen zweiter Ordnung. Als Bei- 
‘spiele seien genannt die Einschwingvorginge eines Schwingungsgebildes unter dem Einfluf 
einer veranderlichen auBeren Kraft und die Verzerrungen von nichtperiodischen Spannungs- 
ainderungen am Ausgang eines Vierpoles gegentiber dem Spannungsverlauf am Eingang (elek- 
tronische MeBgerate). 

Ist die Stérfunktion als einfache mathematische Funktion darstellbar, dann kann in 
vielen Fallen mit einem mehr oder minder groBen Aufwand an Rechenarbeit die Lésung ge- 
funden werden. In der Praxis treten aber auch Differentialgleichungen auf, deren von den 
unabhangigen Variabeln abhangige Funktionen empirisch ermittelt worden sind. Die Funk- 
tionsverlaufe sind dann fehlerbehaftet, und es geniigt, fiir die Lésung der Differentialgleichung 
ein Naherungsverfahren anzuwenden. Das folgende graphische Naherungsverfahren gestattet, 
bei geringem Zeitaufwand in einem Zeichnungsgang die gesuchte Lisung zu finden, wenn fiir 
einen Wert der unabhangigen Variabeln die Randbedingungen vorgegeben sind. 

Die Ausgangsgleichung sei: 

d*y 


Pil2) + Al) & + ¥ = a(x). (1) 


Sie wird ersetzt durch die Differenzengleichung 


ieee is 
1 (Ay, J 4 A 
a(x) = ¥ + falx)°-5 ( a iz) Efi(2)-— (2) 
oder 
g(x) _ Ay. F(x), fo(*) Ay, (fil%) f(x) 
op hack igs WA a Bie 2 Ax \ Ax of 
Setzt man 
( =H) a (3) 
dann wird schlieBlich 


Die graphische Lésung der Gleichung (4) zeigt Abb. 1. 
Die praktische Durchfihrung des Verfahrens erfolgt zweckmaBig auf Millimeterpapier 
unter Zuhilfenahme einer Zeichenmaschine. Das Lineal wird entsprechend der Neigung | bis 2 
(siehe Abb. 1) eingestellt und festgeklemmt. Mit 
dieser Neigung wird von Punkt 5 ausgehend 
eine Linie soweit gezogen, bis im Abstand a (3) 
der Punkt 6 gefunden ist. Nach Lésung der 
Klemmvorrichtung wird das Lineal derart ge- 
neigt, da die Punkte 2 und 6 verbunden 
werden kénnen. Von Punkt 2 bis 3 (4x-Ab- 
stand) wird die Linie ausgezogen. In dieser 
Lage wird das Lineal im Zeichenmaschinen- 
kopf festgeklemmt, da der so gefundene 
Neigungswinkel des Lineales zur Ausfiihrung 
der Hilfskonstruktion des darauffolgenden 
Zeichenschritts benétigt wird. Durch Wieder- 
holung des beschriebenen Verfahrens entsteht 
der gesuchte Kurvenzug. Abb. 1. Graphische Lésung der Differenzengleichung. 


ze. 


parallel 


8 


Azz 


120 Vogel: Ein graphisches Verfahren zur Lésung von linearen Differentialgleichungen. Ingenieur-Archiv 


Der Verlauf von f;(x) und a kann unter oder iiber das Zeichenschaubild gezeichnet werden, 
damit fiir jeden Zeichenschritt die entsprechenden Werte leicht gefunden werden. In vielen 
Fallen ist es zweckmafBig, die Funktion f,(x) von vornherein von der Funktion g aus waagrecht 
aufzutragen. 

Die Differenz Ax ist so klein zu wahlen, daB sie weniger als 1/, bis 1/, der kleinsten auf- 
tretenden Werte von a und f,(x) betragt, wenn gréBere Fehler vermieden werden sollen. 


1y 


Abb. 2, Lésungsbeispiel: Einschwingvorgang. 


Der MaBstab fiir die Funktionen a und f,(x) ist gleich dem Mafstab von x baw. Ax aus- 
zufiihren. Fir die y-Koordinate darf der MaSstab beliebig gewahlt werden; doch ist zu be- 
achten, daB das Rechenblatt nach oben und unten geniigend Raum fiir die Hilfskonstruktionen 
aufweist. 

Um die zu erwartende Genauigkeit des Naherungsverfahrens zu zeigen, seien zwei Beispiele 
angefithrt und mit der mathematischen Lésung verglichen. Fiir die Zahlenrechnungen wurde 


¥ 


LZ) f(a) 
Wales 2 


ae ot wo a 
0 7 a 3 4 ro 6 Vi 


Abb. 3. Lésungsbeispiel: Gleichdimensionale lineare Differentialgleichung. 


ein Rechenschieber mit 250 mm Basislinge verwendet. Das erste Beispiel betrifft einen Ein- 
schwingvorgang. Die Koeffizienten der Differentialgleichung seien konstant. Als Ausgangs- 
gleichung gelte ; 


dy dy 


Hierbei soll die Stérfunktion im Bereich x — 0 bis x = 4,4 den Wert 0 und im Bereich x — 4,4 
bis x oo den Verlauf g(x) = 1,6 sin (2 x — 8,8) haben. Die Randbedingungen fiir x = 0 
seien y = 1,1 und dy/dx = — 0,3; fiir den Beginn des zweiten Bereichs ergeben sie sich aus den 
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Werten y und dy/dx fiir x = 4,4 des ersten Bereiches. Die in Zahlen ausgedriickten Lisungen 
sind fiir den ersten Bereich: 


y = e—%325 * (1,1 cos 0,628 « + 0,901 sin 0,628 x) 


und fiir den zweiten Bereich 
y = — 0,2145 sin (2 x — 8,444) — e—%325 * (0.2376 cos 0,628 x + 2,43 sin 0,628 x). 


In Abb. 2 sind die durch das Naherungsverfahren gewonnene Kurve und einige durch Berech- 
nung aus den vorstehenden Formeln gefundene Werte als Punkte eingetragen. 

Das zweite Beispiel (Abb. 3) stellt einen Vergleich zwischen der Naherungslisung und der 
mathematischen Lésung einer gleichdimensionalen linearen Differentialgleichung zweiter 
Ordnung dar. 


Die Ausgangsgleichung lautet 
2 
0,05 25% 403 22% + y¥=3 402%. 


Werden als Randbedingungen fiir x = 1 die Werte y = 2 und dy/dx = 0,5 gewahlt, so ist 
die Lésung 


y =3 + 0,154 x — x—*5. [1,154 cos (In x*78) + 0,685 sin (In x3,708)) , 


0,05 x? 0.3 x 
Ax 

Aus den beiden Beispielen ist zuersehen, daf die Naherungslésung mit den mathematisch 
ermittelten Werten gut ibereinstimmt. 

Das Verfahren ist bei der Lésung von linearen Differentialgleichungen erster Ordnung 
selbstverstandlich ebenfalls anwendbar. AuBerdem sei darauf hingewiesen, dafs dieses gra- 
phische Naherungsverfahren nicht nur dann erfolgreich angewandt werden kann, wenn die 
verdnderlichen Koeffizienten und die Stérfunktion durch empirische Verfahren gewonnen 
wurden, sondern auch wenn die exakte mathematische Liésung einer Gleichung nur mit groBem 
Rechenaufwand zu finden ist. 


Der Verlauf von f,(x) = 0,3 x unda = 


ist in Abb. 3 eingetragen. 


(Eingegangen am 13. April 1954.) 
Anschrift des Verfassers: Dr. Wilhelm Vogel, Augsburg, Schillstr. 91. 
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Der ebene Verzerrungszustand des dickwandigen Rohres 
bei einem nichtlinearen Elastizitatsgesetz. 


Von F. Jindra. 


1. Einleitung. Es handelt sich hier um die Berechnung der Spannungsverteilung in einem 
dickwandigem Rohr unter gleichmaSigem Innen- und Auf endruck. Die elastischen Eigen- 
schaften des homogenen, isotropen Werkstoffs werden durch einen nichtlinearen Zusammen- 
hang zwischen Spannungen und Verzerrungen beschrieben. Das nichtlineare Elastizitatsgesetz 
fir kleine Verzerrungen kann allgemein in der folgenden Form geschrieben werden ?: 


oe = 3 K x(&) & +26 (yd) (ee — &) 
Oy = 3 K x(é) & + 26 y(yo) (€y — &) » 
0: = 3 K x(é&) & + 2G y(po) (&« — &) » 
Try =GY(Yo) Yeys Tye =EY(Po) Pyzr Tex = © (YO) Pex? 
wobei, bezogen auf ein rechtwinkliges (x, y, z)-Koordinatensystem, 6x, Oy, Oz) Tx y Tyzs Tex die 


Spannungskomponenten, €,, &y, Ez, Pxys Pyz> Yex die Verzerrungskomponenten, K und G den 
Kompressionsmodul und den Schubmodul, ¢) und y{ die invarianten VerzerrungsgréBen 


oy = = (ee + by + &2), | 


(1) 


(2) 


wash [(ex — &)? + (€y — &)? + (ez — €0)*] + Pay + Pye 4 Pay | 


bedeuten. Die Dehnungsfunktion~x (eg) und die Scherungsfunktion y (y}), 
die experimentell zu ermitteln sind, nehmen als Potenzreihen von é) und y@ die Gestalt an: 
H(€&) = 1 +x, +x eG +uze8+---, | (3) 
yy) =1+yevetyevitveyit---- J 
Im folgenden geben wir zuerst allgemein die Differentialgleichung fiir die radiale Verschiebung 
an (Ziff. 2). Beschranken wir uns auf bestimmte einfache Formen der Funktionen (3), so wird 
die Weiterbehandlung der aufgestellten Differentialgleichung durch Reihenentwicklung (Ziff. 3) 
und durch sukzessive Naherung (Ziff. 4) erméglicht. Die Zahlenrechnungen ergeben, dafi eine 
geringfiigige Abweichung des Elastizitatsgesetzes von der Linearitat eine recht erhebliche Ver- 
minderung des Wertes der Tangentialspannung am Innenrand verursacht. Die vorliegende 


Untersuchung kann somit als ein weiterer Beitrag zur Erklarung des Abbaues der Spannungs- 
spitzen angesehen werden. 


2. Die Differentialgleichung fiir die radiale Verschiebung. Wir betrachten einen Kérper von 
der Gestalt eines rotationssymmetrischen Hohlzylinders, der unter gleichmaBiger Belastung 
durch Innen- und Aufendruck steht. Das Koordinatensystem legen wir so, daB die z-Achse 
mit der Zylinderachse zusammenfallt. Wir setzen einen ebenen Verzerrungszustand voraus 2. 
Man spricht von einem ebenen Verzerrungszustand, wenn die Verschiebung in der z-Richtung 
verschwindet: 

es = 0 (4) 
Dieses Ebenbleiben der Querschnitte muB durch eine an den Stirnflachen des Zylinders gleich- 
zeitig auftretende Langsspannung o, erzwungen werden. Das Problem hangt somit nur vom 
Radius r ab. Wie das Kraftegleichgewicht an einem aus der Zylinderwand herausgeschnittenen 
Element im Abstand r vom Mittelpunkt (Abb. 1) zeigt, gilt zwischen Radialspannung o, und 
Tangentialspannung o, die folgende Beziehung: 


do, 
Og = 0; pr a (5) 
* H. Kauderer, Ing.-Arch, 17 (1949) S. 450; F. Jindra, Ing.-Arch, 22 (1954), S. 121. 
pe ex. praktisch wichtigere Fall des ebenen Spannungszustandes im dickwandigen Rohr bei einem 
nichtlinearen Elastizitatsgesetz ist von mir bereits behandelt worden. Die Arbeit erscheint demnichst in 
der Zeitschrift fiir angewandte Mathematik und Physik (ZAMP), 6 (1955). 
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Bezeichnet u die radiale Verschiebung, so ist die Formanderung vollstandig durch die Deh- 
nungen in radialer und tangentialer Richtung 


_ du =a 
of aps ? or Ss (6) 
gegeben. Dann wird nach (2) mit (4) und (6) 
1 1 fd 
wobe red 443) | 
5 ee Ok ae 8 [/du\2 d (7) 
Wo = foi (EF 8). Ew -[- E>) — 9 (a) : = (*)] . | 


Mit Hilfe der Spannungs-Dehnungsgleichungen (1) kann man die Spannungen durch die 
Dehnungen ausdriicken und erhalt mit (7) 


a, = Kx(e) (7 + 7) + 3 Crwd(27 — +). 
eK du u 2 aw u 
Oy = Kx(e) & =f = oa 3° Y(Yo) (3 —='2 7) 3 (8) 


0. = Kx(%) GF eng PWD Ge 


r 


Man geht mit den Ausdriicken (8) in (5) ein und findet 
d diene w d > du u 
3K az |e) (z; al + 26— y (po) (2 ae se >| 


ae ete 
+ 6G y(y) (z: as . = 


r 


Setzt man schlieBlich noch fiir die bezogene radiale 
Verschiebung 


Abb. 1. Das dickwandige Rohr. 


u 
i nes > (9) 
so folgt fiir 2(r) die Differentialgleichung 
d d d d d 
3K ar [*e) (: ane “)|+ 262 [v(vi 2 oe :) + 66 (yi) = = 0 (10) 
mit 
dz\2 d 
o=3("e +24). v= >|" (a) trad +a. (10a) 


Bezeichnet man mit r, und ry den auferen und den inneren Halbmesser, mit p, und p, den 
AuBen- und den Innendruck, so sind die Randbedingungen an den beiden Mantelflachen des 


Zylinders : 
(0;)r=r, seme ee OF, ai (0; ar. hice a Lie (11) 
Die Lésung der nichtlinearen Randwertaufgabe (10) mit (11) sté®t bei willkirlich vorge- 


gebenen Dehnungs- und Scherungsfunktionen auf uniiberwindliche mathematische Schwierig- 
keiten. Wir beschrinken uns von jetzt an auf solche Werkstoffe, bei den die Funktionen (3) 


in der Form ; : 

x(é) = 1, (yo) =1 + Ye Yo (12) 
darstellbar sind. Dies ist fiir viele metallische Werkstoffe der Fall!. Dann vereinfacht sich 
die Differential gleichung (10) fiir die bezogene radiale Verschiebung zu 


2 dz\2 d - d 
ean (Zh teersh bast +ilrge + [en (7) +5arz— +302 +3) =0, (13) 


wobei zur Abkiirzung die (negative) Stoffzahl 


7 3K +48 


a 


1 Vgl. H. Kauderer, a. a. O. S. 468. 
Qg* 
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eingefiihrt wird, die hier allein die Abweichung vom Hookeschen Gesetz kennzeichnet. Mit der_ 


Annahme (12) und wegen (9) nehmen die Spannungskomponenten (8) die Gestalt an: 


a, = K(22 7 | 56 (r+ ari) by GK +40)e((s +r FZ) +e(F)] 


dr 3 
d 4 d 1 dz\3 
oy — K(22 | ri) | 56 (: ied | y BK +46) a[2 — (5). (15) 


Ca eee 4\ oes} (5 Ko dena st 
OC araa — (22 +r FZ) — a a ) as +(s+r3)]- 


3. Lésung durch Reihenentwicklung. Die Differentialgleichung (13) laBt sich durch ge- | 


eignete Transformationen in eine Abelsche Differentialgleichung zweiter Art iiberfiihren. Da 
trotzdem die Gleichung nicht geschlossen gelést werden konnte, suchen wir die Lésung in 
Gestalt einer Reihenentwicklung. 


Wir wollen die Integralkurve von (13) in der Umgebung der Stelle r = r, durch eine Potenz- | 


reihe darstellen. Wir setzen daher in (13) 
sues Pee (16) 


Ta 


und machen fiir z(r) = z(1 | = 2(s) den Ansatz 


a 


lo ©) 
&(s) =D) 54 Shey iaaey S ieee 8? aS ea ae (17) | 
p=0 


Geht man mit (17) in (13) ein, entwickelt nach Potenzen von s und vergleicht die Faktoren | 


entsprechender Potenzen von s, so bestimmen sich die Koeffizienten c,, der Potenzreihenent- 
wicklung rekursiv aus dem folgenden Gleichungssystem: 


2 [a (2 ce — 2c, ce) + ¢3) + l]e, —a (4ci —S5e, ce +3 ci) c, —3¢,=0, (18a) 
6 [a (2 c2 — 2c, cy + cf) + 1] cg + 2 a (8 ce, c, — 4, cy — 20 cf + 16, cy — 4 6) cg — 8 ey 
“sa (13:6, — 11.e,)\c? ==-0., (18b) 
12 [a (2 2 — 2e,¢) + 3) + 1] eg + 3a (24 c,¢, — 12 eye) — 28 ci 4+ 22 c cy — 5c?) cg — 15, 
+ a (16 cF — 108 c, cy + 40 €g cy + 95 cf — 46 €, cy) cg — 12 ac? = 0, (18c) 


20 [a (2 cf — 2 cy ¢y + 6) + le, + 4 (32 cy ce, — 16 cg cy — 36 ci + 28 c, co — 6 €2) c, — 246, 
+ a (72 ¢3 cy — 36 cg cy + 96 c} — 412 cy cy + 148 cy cy + 173 cf — 78 cy G) C3 

— a (88 cf — 212 coe, + 44¢e,.¢, + T1c2)c, =0, (18d) 
30 [a (2 cf — 2c eg + 6) + lo eg + 5 a (40 c, c, — 20 c, cy — 44 cc? + 34 €, cg — 7 €2) cs — 35 cg 
+ a (240 3 ¢, — 120 cy eg + 160 c3 — 668 cy cy + 236 cy cy + 273 cf — 118 c Cy) cy 
+ @ (180 cs ¢. — 372 €3 ¢, + 129 cs cg — 476 c5 + 734 cy c, — 142 cg cy — 116 €?) cg 
+15a(10c,—9¢)2=0, : (18e) 


42 [a (2 cf — 2 cy eg + €G) + 1] ec, + 6a (48 c, c, — 24 ¢, cg — 52 c? + 40 c, cy — 8 €2) cg — 48 C, | 


+ a (360 cg cy — 180 €3 ey + 240 c5 — 984 c, c, + 344 cg cy + 395 c? — 166 c, G) 5 

+ @ (192 cy cy — 96 ¢4 eg + 576 3 Cp — 1164 cg c, + 396 cs co — 744 c2 + 1120 c, c, 

— 208 ey ¢g — 171 cf) cg + @ (108 3 — 828 cy c, + 615 cg cy — 111] Cy eg + 754 c2 

— 430 c. cy) c; — 84ac8 = 0, (18f) 
56 [a (2 cf — 2 cy ey + cg) + 1] eg + 7 a (56 c, cy — 28 cy cg — 60 c? + 46 c, cy —9 €2) c, — 63 c, 

+ a (504 ¢3 ¢, — 252 cs cg + 336 cy — 1360 cy cy + 472 cy cy + 539 c? — 222 c, Co) es 

+ @ (560 cy cy — 280 cy cg + 840 cg cp — 1672 cy cy + 562 cg Cy — 1068 c} + 1582 c, ¢, 

— 286 c, ¢) — 236 ci) c; + a (448 cy cy — 888 cy cy + 296 cy cy + 504 c2 — 2524 c, c 

+ 1834 ¢; ¢, — 318 cs cg + 1124 cf — 622 cy ¢) cy — a (468 c2 — 1235 cy c,+ 336 Cy Cy 

+ 395 c3)c, = 0, (18g) 
usw. Die Auflésung dieser Gleichung ergibt 


Co = sh [a (4c} — 5 ec + 303) +3), (19a) 
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C. 
a =e [a (44 cf — 148 cicy + 252 cf c§ — 254 c3 c8 + 163 c? ch — 63 c, c5 + 12 c8) 
+ a (84 ci — 188 cP cg + 214 c? c2-— 126 c, c38 + 36 cf) + a (51 ci — 63 cy cy + 36 c2) 


+12], (19b) 
Cc 
ee ae [a> (640 c{° — 3432 cf cy + 9352 cf c§ — 16412 cf c8 + 20 352 c8 ch — 18562 c c3 
+ 12622 cf cf — 6339 c} c% + 2268 c} c8 — 525 c, c§ + 60 ci) + at (1920 c8 
— 8244 cj cy + 17976 cf cf — 24800 c& c8 + 23528 cf cd — 15653 c3 c5 + 7200 c? cf 
— 2100 c, 3 +300 c&) + a3 (2224 c§ — 7322 8 cy + 12218 cfc? — 12.289 c3 c2 
+ 1992 cf c§ — 3150 c, c3 + 600 c8) + a? (1312 cf — 2975 c3 cy + 3456 c? ec? 
— 2100 c, cf + 600 cd) + a (396 c? — 525 c, cy + 300 c2)+ 60], (19c) 
c 
C5 = zyatna a” (11872 cl! — 86560 cl? cy + 320336 cl? c2 — 781.024 cl! 8 + 1388144 cl? cf 


— 1892656 c? ce} + 2036520 c8 c§ — 1755696 c? cf} + 1219874 cf c8 — 618474 c§ cf 
+ 302489 ct cf? — 103 980 c} cit + 26295 c? cl? — 4380 c, cf3 + 360 cf?) 

+ a (48.016 cl? — 300960 cl cy + 955392 cl c2 — 1984272 c? ck + 2975928 c8 cf 
— 3378624 cf c§ + 2974112 c8 cS — 2047740 cic? + 1098321 c# 3 — 449970 c3 c8 
+ 134430 c? cf? — 26280 c, cit + 2520 ci”) + a5 (81120 c}° — 427664 c? cy 

+ 1141272 c8 2 — 1968888 c? c3 + 2414464 cb ch — 2191980 c3 c8 + 1494906 ct cf 
— 760080 c? ci + 277725 c? c8 — 65700 ¢, c8 + 7560 ci) + at (76072 cf 

— 324.656 c7 cy + 703096 c$ cz — 966636 c3 c3 + 919682 ch cd — 620220 c? c3 

+ 292500 c? c6 — 87600 c, cZ + 12600 c8) + a3 (42870 c§ — 140922 cf cy 

+ 235.485 cf c2 — 240180 c3 c2 + 161025 2 ch — 65700 c, c§ + 12600 c) 

+ a® (14877 cf — 34.050 c3 c, + 41070 €2 c2 — 26280 c, 3 + 7560 c?) 

+ a (2955 c} — 4380 c, cy + 2520 c2) + 360], (19d) 


usw. mit der Abkiirzung 
N=a(2e7 —2e¢4 +2) +1. (19e) 
Wir miissen nun fiir die Lésungsfunktion Z (s) die Konstanten cy und c, gema8 den Rand- 


bedingungen (11) bestimmen. Mit der dimensionslosen unabhangigen Veranderlichen s folgt 
fiir die radiale Spannung nach (15) 


o,=5(3K _ G)z ~ (3K + 4G) (1 = 8) % +4 (3K +460 {fe—(1—s)2'P 
as (1 aes 8)* og (20) 


wobei Striche Ableitungen nach s bedeuten. Fir r = r, wird 
B10) Se, 2: (0) = cy. (21a) 


Weiter setzen wir zur Abkiirzung 


S=2(1—T)—Se,or, = RF (1—B\= 0-9 FW + Nero" (21d) 


Ta w=0 Ta n=0 
fe r 

mit o=1-—-. 
Ta 


Unter Benutzung der Beziehungen (20) bis (21b) entstehen dann aus (11) fiir cy und c, folgende 
Bestimmungsgleichungen: 


ee SKE AC) op (SK E46) a [(og 22)? = YS es, 
3 3 9 (22) 


2K +6)S 10K +46)5, 443K 446) a[(5,—S)— Sf] =—p- | 


Da die Auflésung des Gleichungssystems (22) zu umstandlich wird, behandelt man die 
Randwertaufgabe wie eine Anfangswertaufgabe. Von einem angenommenen ¢ ausgehend, 


126 Jindra: Der ebene Verzerrungszustand des dickwandigen Rohres usw. Ingenieur-Archiv | 


berechnet man zunichst ¢) aus der ersten Gleichung (22), die man in der Form schreibt: 


(co — 4)? + 3a (c —%) +28 = 0 (23) 


mit 


3K 2G 
a. Bee aaGe 


ads! 9 Da 5 R206 b 
eerrite eee r3oKa ac & @ ci) « | 


Oh = 
(23a) 


Sobald die reduzierte kubische Gleichung (23) aufgeldst ist, lassen sich weitere Koeffizienten | 
der Reihenentwicklung nach (19a) bis (19e) und somit die Ausdriicke (21b) berechnen. Aus der | 
zweiten Gleichung (22) folgt dann der zu dem angenommenen c¢, gehérige Innendruck py. Die | 
Rechnung wird fiir verschiedene Werte von c, wiederholt und durch Interpolation ein solches c, | 


ermittelt, fiir welches c, und c, die vorgegebenen Randbedingungen erfillen. 


Die Tangentialspannung o, findet man endlich gemafB (15) mit (16) zu 


op =~ BK 46)3—2 (8K —26)(1—s)¥ + | (8K +46) a [ze + (1 — 5) 79). (24) 


3 


Die Randwerte der Tangentialspannung o, an den Stellen r = r, undr = ry ergeben sich dann | 


mit (21a) und (21b) zu 


2 s 1 il : 
(Op)rare = J BK + 6) —3 (3K — 26) 4 + 9 (3K +46) a(c§ + ¢}), | 


(24a) 


bo 


Cae ON eA OG Koa), +e GK +46)a(55 +5). | 


Ebenso errechnet man nach (15) mit (16) die an den Stirnflachen des Zylinders anzubringende 
Langsspannung o, zu 
2 


o, =~ (3K —26)% — = (8K —26)(1—s)2’ —4 (8K +46)a {#8 + [2 —(1—s)7'}. (25) 


Wir zeigen noch, daB die Lésungen (20), (24) und (25) fir a — 0 in die bekannte Lésung 
der linearen Theorie itbergehen. Nach (19a) bis (19e) erhalten wir fiir die Koeffizienten der 
linearen Theorie 


1 
Giana (eae) er fiir uw > ‘Ls 
Die unendliche Reihe (17) wird dann 


= isco) 1 1 os) 
a(S) == hese =ant, Fas Oh ae tin (i ed) sl, 
t 


0 1=0 


Nun gewinnt man aber durch gliedweises Differentieren der geometrischen Reihe die Beziehung 


d Lf . Cl Se 1 
dx aie ee Pas (l— x} 
Damit folgt 
~ 1% LAlts il cx 
(Cesta ce Neae (a nace 
oder 
1 ra 
2*(r) = 5 — 2 ct (1 = S ® (26) 
Die Randbedingungen (11) verlangen bekanntlich 
be ol Poe Paes Li Pes Pate LP eRe 
0 P36 RU eee Ke eR cue aq) 
mit R=, (27a) 


9 
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he Spannungskomponenten der linearen Theorie nehmen damit in der Tat die bekannten? 
erte an 


1 va ; 1 te 
Toma |e Pate Pe 7 Ps) slo Se = fay Py — Pa R* + (Py — pe) Za], | 
(28) 


eo R21 G pe Pal 
O; 


o* = 


31-26 (hk? = 15: 
und insbesondere die Randwerte 


1 
(89) a1 = PL] [2 Po as Ce Se a Ss [Petit yi 20p, R*) 2) (29) 


; Als Beis piel eines Werkstoffes, das ein nichtlineares Verhalten zeigt, betrachten wir 
hier gezogenes Kupfer, fiir welches 


K =1,37-10%kg/em?, G = 0,46-10%kg/em?, y, = — 0,18- 10° 


errechnet wurde®. Es handelt sich also dabei um sehr kleine Abweichungen vom Hookeschen 
Gesetz. Abb. 2 zeigt den nach den Formeln (20), (24) und (25) berechneten Verlauf der Radial- 
spannung o,, der Tangentialspannungen o, und der Langsspannung o,, und den Verlauf von oF 


6; 60; Oz 


und oj nach (28) lings eines Halbmessers (und zwar fiir 
R=1,25 und py) = 200 kg/cm?, p, =0). Man sieht, daB die 
Spannung infolge der Nichtlinearitat gleichmaBiger verlauft 


(5p)r-r, 


600 


400 


200 


[0 
0 50 700 160 200g [om 


Abb. 2. Spannungen im dickwandigen Rohr. Abb. 3. Die Abhangigkeit der Tangentialspannung vom Innendruck. 


und so das Material besser ausgeniitzt wird. Um eine Anschauung von der Art der Abhangig- 
keit des Ausdruckes fiir (o,),-,, vom Innendruck py zu erhalten, wobei wieder pz =0 gesetzt 
ist, ist in Abb. 3 (c,),—,, abhingig von p, fiir verschiedene Werte von R graphisch dargestellt. 
Mit steigendem R weichen die Kurven (a,),—,, immer mehr von den Geraden der linearen Theorie 
ab, die gestrichelt eingetragen sind. Der Einflu8 der Nichtlinearitat wird also bei gréBeren 
Halbmesserverhaltnissen starker. 


4, Lésung durch sukzessive Naherung. Fiir gréBere Werte von R konvergiert die Potenz- 
reihenentwicklung von Ziff. 3 recht langsam. Wir geben deshalb noch die Lésung der nicht- 
linearen Differentialgleichung (13) durch sukzessive Approximation an. 


Zu diesem Zweck bringen wir (13) auf die Form 
on rie a(2rz’ + z)rz” 
ee ya ~ Qar?z? + Qarzz’ +az?4+ 1° et) 


Striche bedeuten jetzt Ableitungen nach r, 


1 Vgl. z. B. S. Timoshenko, Theory of Elasticity, 1. Aufl., S$. 56, New York 1934. 
2 Nach Versuchsergebnissen von K. Claus, Bestimmung der Poissonschen Zahl reiner Metalle und ihr 
Gang im Periodischen System, S. 39. Diss. Stuttgart 1953. 
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Fiir die Naherung nullter Ordnung behandelt man die Randwertaufgabe der homogenen 
Gleichung 


rz)" + 3 2(0)’ — 9, (31) 
deren Lésung sofort geschlossen angebbar ist: 


Oe a (32) 


Sie stellt auch die Lésung der linearen Theorie dar, die sich fiir a = 0 ergibt. Die Integrations- 
konstanten A und B sind aus den Randbedingungen zu bestimmen. Unter Beriicksichtigung 
von (11) und (15) ergeben sich zwei nichtlineare Gleichungen fiir A und B. Sie lauten 


*(3K +6) A —2GB’ +4 (3K +46)a[(A — B’)) — 8B] = — p,, | 
(33) 
+ (3K +6) 4 —2GR°B’ += (3K +4G)a[(4 — RB’)? —8 RB's] = — py | 
mit den Abkiirzungen 
/ B a 
B=5, erst (33a) 


Den Lésungswert 2() setzen wir in die rechte Seite der Differentialgleichung (30) und er- 
halten somit fiir die erste Naherungsfunktion z(!) die Gleichung 


4 a B? ( Ar? — 3 B) 1 


(Cy @)? = é 
i) ae (a A?+ 1)r4*—2aA Br? + 5a Br? - 


(34) 


Diese Differentialgleichung integrieren wir nach der Methode der Variation der Konstanten, wo- 
bei die auftretenden Integrale geschlossen auf elementarem Weg ermittelt werden kénnen. 
Mit den Integrationskonstanten C und D ]a8t sich die Lésung von (34) in der Form schreiben 
D 
Baer eT) 2 (35) 
wobei zu setzen ist 
: 1 B 
f= 50{(4 ys S In 
u 
V—a(5 + 4a A?) 


(aA* +1) 20454507) —(15 + Mad? — 1004 2) 
r iB 


(a A? + 1)r2?—aAB+ BV—a(5 + 4a A?) 
(a A? + 1)r2?—-aAB—BY—a(5 + 4a A?) 


, 66 
Die Ableitung von f(r) ist 


PHO 546 +8In|(a.4? 41) 20454505 


gaAy hn, eA ee 
AB (a A? + 1)r2?—aAB— BV—a(5 + 4a A |} 


Ves - 4e A) 
Nach (15) erhalten wir jetzt fiir die Naherungen erster Ordnung der Spannungskomponenten 


2 F 
6) = 38 K + G) z@) + = (3 K+ 4G) r2Q)’ + 5 (3K + 46) a[(2 4 r2()')3 + 73 2(0)’3], } 


eens 1 et I 
of) = 4 (3K + G) 2(0) tay le Roe 26) rz) cae, (3K + 46) a [23 — 73 20)’3] , | (37) 
6) = Sik — 2G) (220) + rz’) — $ (3 K +46) a [23 + (20) + r2(0)’)3] . 
J 


Die Randbedingungen (11) liefern zwei nichtlineare Gleichungen zur Bestimmung der Inte- 
grationskonstanten C und D. Sie lauten 


5 (3K +6)(C + D' + f(r] —4(8K +46) [2D! — ref'(ra)] 


1 ; (38) 
++ 9 3K +46)a {(C Dy (52) +rof’(ra)}? — [2D’ — Taf '(Ta)]*} = — Pas 
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+ (83K +6)[C+ RD’ + f(r)] —+ (3K +46) [2 RD! — ny f'(r)] | 


1 
3 


1 (38) 
Ss 73 K+G)a {[C —R? D'+ f(t) +70 f ‘(79)]?-— [2 R? D’ — tof (To) }* ae eal Ke | 
D=5 (38a) 


gesetzt ist. 


Somit hat man die Naherungslésung erster Ordnung fiir die Differentialgleichung (30) unter 
den Randbedingungen (11) gefunden. Diese Rechenvorschrift der sukzessiven Naherung 
kénnte weiter fortgesetzt werden, so dab die Verschiebung sowie die gesuchte Spannungsver- 
teilung beliebig genau angendhert werden kénnte. 

Zum Schlu8 danke ich Herrn Prof. R. Grammel sowie der Deutschen Forschungsgemein- 
schaft fiir die Méglichkeit, diese Untersuchung durchfiihren zu kénnen. 


(Eingegangen am 21. Mai 1954.) 
Anschrift des Verfassers: Dr.-Ing. F. Jindra, Stuttgart-O, Gerokstr. 70. 


130 Okubo: Die Formzablen tordierter Wellen mit mehreren Nuten. Ingenieur-Archiv 


Die Formzahlen tordierter Wellen mit mehreren Nuten. 
Von H. Okubo. 


1. Einleitung. In der vorliegenden Arbeit wird ein einfaches Verfahren zur genaherten Be- 
rechnung der Formzahlen fiir tordierte Wellen mit mehreren Nuten mitgeteilt. Ausgehend von 
der exakten Lésung fir halbkreisformige Mehrfachnuten, ermitteln wir zuerst die Formzahl 
fiir unendlich kleine Mehrfachnuten (Ziff. 2) und leiten daraus ein praktisch brauchbares Nihe- 
rungsverfahren fiir endlich ausgedehnte, aber hinreichend kleine Mehrfachnuten her (Ziff. 3). 
Seine Giite wird durch Vergleich mit den exakten Ergebnissen aus numerischer Rechnung 
gepruft. sete : 

Die von uns gefundene Formzahl unterscheidet sich schon bei unendlich kleinen Mehrfach- 
nuten wesentlich von einem von H, Neuber 1 angegebenen Ausdruck. Daher sei zunachst dieser 
kurz entwickelt. Fiir eine Einzelnut mit starker Kerbkriimmung 1/9 und der Tiefe t betragt 
die Formzahl Vi/o. Reihen sich unendlich viele solcher Nuten im gleichen Abstand b aneinan- 
der, so wird die Formzahl 

= : b mt 
a= Vrtie mit y=—39-; (1) 
der Faktor y < 1 ist die Entlastungszahl. Man erhalt also a aus der Formzahl fiir die ent- 
sprechende Einzelnut, indem man diese im Verhaltnis 


R, = yy (2) 
vergréBert oder auch, indem man bei der Mehrfachnut die wirkliche Tiefe t durch eine wirksame 
Tiefe y t ersetzt, im tbrigen aber den Ausdruck fiir die Formzahl der Einzelnut beniitzt. Auf 
Grund von rechnerischen Untersuchungen schlieBt H. Neuber, daB sich dieses Vorgehen nicht 
nur bei starker Kerbkriimmung, sondern ganz allgemein beim Ubergang von der Einzelnut | 
zum mehrfach genuteten Rand bewiahre. Er setzt so z. B. fiir flache Mehrfachnuten endlicher | 
Krimmung als Formzahl 


a =1+ Yyto (3) 

an, d. h. er vergréBert die Formzahl der entsprechenden Einzelnut im Verhaltnis 
Raitt 4e 
1+ tle 


Die in (3) ausgedriickte Neubersche Regel ist heute in der Praxis in Gebrauch. 


(4) 


2. Unendlich kleine Mehrfachnuten endlicher Kerbkriimmung. Wir betrachten eine Folge 
von sehr kleinen Nuten, deren Tiefe t und Kerbkriimmungshalbmesser @ sich nur wenig von- 
einander unterscheiden (Abb. 1). Besitzt die auf Torsion be- 
anspruchte Welle vom Halbmesser Eins m solche Langsnuten, _ 
so ist die gré8te Schubspannung im Kerbgrund der Nuten ? 


m—1 2 
mc le 
Abb, 1. Unendlich kleine, fast halbkreis- {hanes EW) le a ( ) ; 
formige Mehrfachnut. 1 —c%m 


(5) 


hier ist G der Schubmodul, w die spezifische Verdrehung der Welle und c = 1 — t. Dies ent- 
spricht bei m sehr kleinen Langsnuten einer Formzahl 
Tmax me 1 (1 = 6) 
ica cape ee VEC. (6) 
Daraus erhalten wir durch Grenziibergang m —> oo bei festgehaltenem Wert m t leicht die Form- 
zahl fiir unendlich viele, unendlich kleine Nuten zu 


0 Sane Qt) (7) 


‘ Fiir das Folgende vgl. H. Neuber, Kerbspannungslehre, S. 139 u. 141, Berlin 1937. 
2 H. Okubo, J. Appl. Mech,, Trans, ASME, 17 (1950), S, 359. 
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In Polarkoordinaten r, } mit demin Abb. 2 angegebenen Polin der Nutenmitte lautet dabei 
die Gleichung fiir das Nutenprofil 


add of lbh ages 
set Bal, : ) 


2ut 4nr 2a0r f 
Sees — cos @ — cos } ‘ 
ae ee “yh AR eeet Lie ore eel ued (8) 
und der Kriimmungshalbmesser im Kerbgrund der Nut betragt 
pee St, CANAL ut 


Wenn der Nutenabstand 6 grof gegen r oder die Nutentiefe ¢ ist, vereinfacht sich (8) zu r = t, 
und es liegen dann halbkreisférmige Nuten vom Halbmesser t vor. Fiir nicht hinreichend groBe 
Werte b ist dagegen das Verhiltnis 9/t +: 1 und ist aus (9) zu be- 
rechnen; doch zeigt Abb. 3, daB der Wert o/t nur fiir sehr nahe 
nebeneinander liegende Nuten (b/t klein) merklich von Eins ver- 
‘schieden ist. In allen anderen Fallen liegen praktisch halbkreis- % 
férmige Nuten vor. io 

Fur das Formzahlverhaltnis bei unendlich vielen Nuten folgt 
nunmehr aus (7) und (3) (mit y = 1) 


R, = ( E 2 - n=) ! Vt/o) b (10) Abb, 2. Polarkoordinaten, 


b Sin —— 
Sint b 


wo 0/t tiber (9) in Abhangigkeit von b/t gegeben ist. Abb. 3 zeigt den Verlauf von R, und zum 
Vergleich die Neuberschen Werte R (4) als Funktionen von b/t. Dazu ist noch die Kurve R, (2) 
gestrichelt eingezeichnet, die fiir stark gekriimmte Mehrfachnuten gilt. Es mu8 auffallen, daB 
sich die Neuberschen Werte R (4) durchweg stark von den Werten R, unterscheiden und auBer- 
dem, entgegen der Erwartung, nicht in das zwischen den Kurven Ry und R, gelegene, schraf- 
fiert gezeichnete Gebiet fallen. Die Neuber- 
sche Regel zur Gewinnung der Formzahlen fiir 
Mehrfachnuten fiihrt also in unserem Beispiel 
auf Formzahlverhaltnisse R, die sich nach der 
verkehrten Seite auswirken. Daher lat sich 
diese aus den Verhaltnissen bei stark gekriimm- 
ten Nuten gewonnene Regel nicht allgemeinauf- — 
rechterhalten, sondern es mu der Sachverhalt 
bei Mehrfachnuten erneut untersucht werden. 

Es sei noch bemerkt, daf unser Ergebnis 
(10) sowohl fir Langsnuten als auch fir gleich 
gestaltete Querrillen (Ringnuten) gilt, daes sich 
um unendlich kleine Nuten handelt. Aus ihm 
wollen wir nun eine brauchbare Naherungs- 
lésung fiir endlich kleine Querrillen herleiten. 


3. Endlich kleine Querrillen. S. Momma’ 4, 
hat an Hand numerischer Rechnungen gezeigt, : Z i $ : ei Rint mie 
daB fir verhiltnismaBig kleine Querrillen der — “”” Nuiburcisfirmige, onendlich kleine Mehrfachnuten, 
Faktor R mit abnehmendem Krimmungshalb- 
messer ¢ nur sehr langsam zunimmt, sofern man hierbei die Werte von Rillentiefe und -abstand 
fest 14Bt. Man erhalt daher genauere Werte des Formzahlverhaltnisses R, wenn man zwischen 
den beiden Grenzkurven R, und R, in Abb. 3 interpoliert; denn nach den Mommaschen Rech- 
nungen fallt der Wert von R fir beliebige Werte 0/t < (@/t)) in das schraffierte Gebiet der 
Abb. 3, wobei (0/t), in Abhangigkeit von b/t durch die 9/t-Kurve in Abb. 3 gegeben wird. Wir 
erlautern die Verhaltnisse an Hand von Abb. 4, die nach Rechnungen von Momma fir b/t = 3, 
t/d = 0,05 mit d als kleinstem Wellendurchmesser entworfen ist. Die Grenzwerte R, und R, 
fir unendlich kleine Nuten mit b/t = 3 sind der Abb. 3 entnommen. Wie man sieht, liegen sie 


1 §, Momma, J. Appl. Phys. 24 (1953), 5. 959. 
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ziemlich nahe bei den entsprechenden Werten fiir endlich kleine Nuten. Die vergleichsweise 
kleine Unstimmigkeit zwischen beiden Werten R an der Stelle g = 0 ist hauptsachlich den von 
Momma gewahlten Nutenprofilen zuzuschreiben. Ihre Form bleibt nicht konstant, sondern 
aindert sich so, da® die Spannung im Kerbgrund der Nut bei abnehmendem Kriimmungshalb- 
messer etwas griBer als bei Nuten ausfallt, deren Form sich mit g nicht andert 1, Daher ist der 
richtige Wert fiir das Formzahlverhaltnis R an der Stelle 9g = 0 etwas kleiner als der in Abb. 4 
eingezeichnete Mommasche Wert anzusetzen und wird so naher bei Ry liegen. Bei geniigend 
kleinen Nuten werden daher die Werte von R, und R, praktisch nicht von dem GréSenverhalt- 
nis der Nuten zur Welle abhangen. Nehmen wir weiter an, da sich R linear mit g/t andert, so 
kommt 
R= Ry + (Ri — Bi) (7) (11) 
0 
Abb. 4 zeigt fiir b/t = 3, daB sich die Interpolationsgerade (11) der genauen Mommaschen 
Kurve sehr gut anpaBt. Somit laBt sich das Formzahlverhaltnis R mit Hilfe von (11) fir 
endlich kleine Nuten naherungsweise be- 
stimmen, womit sich die Formzahl fir mehrere 
Querrillen aus der Formzahl o;, fir eine 
Einzelrille gemaB 
OE = Roako (12) 


ergibt. Die Berechnung von dj bereitet 


o_ 
g auch bei maBig groBen Nuten keine be- 
; zt sonderen Schwierigkeiten, da hierfiir ein- 


OF 1 i =e ell A 
0 7 O4 06 08 10 fache Formeln vorliegen ?. 
Abb, 4, Zur Interpolation des Formzahlverhaltnisses R bei Die Formzahl wird gewohnlich als Ver- 


mehrfachen, endlich kleinen Querrillen. 


haltnis der SpannungsgréBtwerte bei gleicher 
Torsionsbeanspruchung der Welle mit und ohne Nuten definiert. Es ist aber zu beachten, dab 
dann die Beziehung (10) bei mehreren Langsnuten endlicher Abmessungen auch nicht nahe- 
rungsweise gilt. Vielmehr ist der auf das Torsionsmoment bezogene Wert des Formzahlver- 
haltnisses R oft viel gréBer als Eins, weil namlich die Torsionssteifigkeit von Wellen mit vielen 
solchen Nuten stark herabgesetzt wird. 

Dagegen bleibt die Beziehung (10) auch fiir maBig groBe Nuten bestehen, wenn man die 
Formzahl nicht auf das Torsionsmoment, sondern auf den Torsionswinkel bezieht. Beispiels- 
weise erhalt man fiir mehrere halbkreisférmige Langsnuten mit c = 0,63 oder alsot = 9 = 0,37 
naherungsweise folgende SpannungsgréBtwerte 3: 


Zahl der Nuten 1 2 3 | 4 
m | 
EminlO'as | 1:63:0) 8G Tb oe? bee 0 SB oF als Sao 


Messen wir den Nutenabstand b langs des Kreises vom Halbmesser c, der die Nuten im Kerb- 
grund berithrt, so wird b = 2 wc/m oder b/t = 2c/mt. Dies fihrt auf die Wertetabelle 


Nee eeeeaser slum dben te te 
Ro) .0,93' 1 0s858 ah aoe 


Die so fiir Nuten maBiger Ausdehnung erhaltenen Werte R stimmen fast mit den fiir unendlich 
kleine Nuten giiltigen Werten R, aus Abb. 3 iiberein. Die Beziehung (10) gilt also naherungs- 
weise ebenso fiir Langsnuten wie fiir Querrillen selbst bei maBig groBer Nutenausdehnung, so- 
fern man die Formzahl jetzt auf den Torsionswinkel, im friitheren Falle endlich kleiner Quer- 
rillen dagegen auf das Torsionsmoment bezieht. 


(Eingegangen am 11. Juni 1954.) 


Anschrift des Verfassers: Professor Dr. H. Okubo, Faculty of Engineering, Nagoya University, Chigusaku, 
Nagoya, Japan. 
1 Vgl FuBnote 1 von S$. 131. 
°F. A, Willers, Z. Math, Physik 55 (1907), S, 251; R. Sonntag, Z. angew. Math. Mech. 9 (1929), S. 33 
H, Neuber, a.a.O., S.90; H. Okubo, J. Appl. Mech., Trans. ASME, 19 (1952), S$. 16. Einen weiteren 
Beitrag gibt die Arbeit des Verf. ,,Stress concentration factors for a circumferential notch in a cylindri- 


cal shaft‘, die in den Memoirs of the Faculty of Engineering, Nagoya University 6 (1954) Nr, 1 erscheint. 
3 Vel. FuBnote 2 von S, 130. 
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Die strenge Theorie anisotroper prismatischer Faltwerke. 
Von D. Riidiger. 


1. Einleitung. Mit der vorliegenden Arbeit! wird die Theorie der beliebig anisotropen pris- 
matischen Faltwerke, also die Theorie der aus einzelnen, anisotropen rechteckigen Platten, 
Zylinderschalen und Staiben zusammengesetzten Tragwerke, deren in Abstainden angeordnete 
Binderscheiben einen raumlichen Spannungszustand erzwingen, aufgestellt (Abb. 1). 

Die ersten Untersuchungen iiber raéumliche Tragwerke, die aus ebenen Tragflachen gebildet 
werden, sind wohl zuerst im Jahre 1920 von H. Schwyzer ? fiir U-, I- und Kastenprofile des 
Stahlbaues durchgefiithrt worden. Als Vertraglichkeitshedingungen der an ausgezeichneten 
Punkten zusammenstofenden Scheiben wird die Gleichheit der Spannungen verwendet. Der 
Gedanke, prismatische Faltwerke im Stahlbetonbau als tragende und zugleich raumab- 
schheSende Konstruktionen zu verwenden, geht auf ein Zusatzpatent der Dyckerhoff & Wid- 
mann AG. aus dem Jahre 1927 zuriick. Diese Tragwerke werden als Abwandlung des Zeib- 
Dywidag-Schalendaches dargestellt, indem die stetige Kurve durch ein um- oder einbeschrie- 
benes Vieleck ersetzt wird *. Zu der gleichen Zeit hat G. Ehlers * auf einem Vortrag vor Inge- 
nieuren der Wayfi & Freytag AG. auf ein neues Konstruktionsprinzip aufmerksam gemacht, 
bei dem Scheiben und Platten zu raumlichen Tragwerken zusammengesetzt werden, ohne 
speziell an prismatische Tragwerke mit Uberhéhung gegeniiber der Mittelkraftlinie zu denken. 
Der Name Faltwerk ist zuerst von H. Craemer ® in einer Arbeit iiber Scheiben und Faltwerke 
als neue Konstruktionselemente des 
Eisenbetonbaues gepragt worden. Die 
Theorie der prismatischen Faltwerke 
ist grundlegend durch die Abhandlun- 
gen von G. Ehlers*, H. Craemer °, 
E. Gruber ®, G. Griining’, J. Goldenblatt 
und E. Ratz®, R. Ohlig® geférdert 
worden?®, 

In allen diesen Arbeiten wird der Berechnung ein Gelenkwerk zugrunde gelegt, bei dem der 
biegungssteife ZusammenschluB der zwischen den aussteifenden Binderscheiben und den Falten 
liegenden Platten durch eine scharnierartige Verbindung ersetzt wird. Eine Vertraglichkeit der 
Platten an den Kanten wird durch die Bedingung gleicher Spannungen gefordert. Dazu wird 
die Flachenbelastung statisch durch eine in den Kanten wirkende Linienbelastung ersetzt und 
die Notwendigkeit gleicher Spannungen durch Schubkrafte erfillt. Die durch den Verschie- 
bungszustand an den Kanten sekundar eintretenden Biegemomente werden auf Grund der 
winkeltreuen Verformung unter Vernachlassigung der Biegemomente und Querkrafte in Rich- 
tung der Erzeugenden sowie der Drillmomente beriicksichtigt. Diese von der Schalentheorie 
entlehnte Trennung der ohne Zweifel zusammenhingenden statisch unbestimmten Aufgabe 
tragt in die Berechnung eine gewisse Willkiirlichkeit hinein. Man kann wohl von den ,,Mem- 
branlésungen“ der Schalen sagen, daf diese in den meisten technisch interessierenden Fallen 
praktisch mit den Partikularlésungen der vollstandigen Differentialgleichungen zusammen- 


Abb. 1. Querschnitt einer Faltwerkkonstruktion. 


1 Dresdner Habilitationsschrift. 

2 H. Schwyzer, Diss. Ziirich 1920. 

3 Fr, Dischinger u. U. Finsterwalder, Bauing. 9 (1928) S. 811. 

4 G. Ehlers, Bauing. 11 (1930) S. 125, 319, 369; Beton und Eisen 29 (1930) 5S. 281, 291. 

5 H. Craemer, Beton u. Eisen 28 (1929) S. 254, 276, 269, 337, 339; 29 (1930) S. 276; 36 (1937) 5S. 297; 
Bauing. 21 (1940) S. 268; Oster. Ing.-Arch., 4 (1950) S. 350. 

6 E. Gruber, Abb. Intern. Ver. f. Briicken- u. Hochbau 1 (1932) S. 225; 2 (1932) S. 276; 3 (1935) 
S. 134; 7 (1944) S. 139; 12 (1952) S. 167; Bautechnik-Archiv 7 (1953) 5. 62; Bauing. 23 (1942) S. 323. 

7 G. Griining, Ing.-Arch. 3 (1932) 5. 319. 

8 J. Goldenblattu. E. Ratz, Beton- u. Eisen 33 (1934) S. 369. 

9 R. Ohlig, Diss. Darmstadt 1934; Ing.-Arch. 6 (1935) S. 346; 12 (1941) S. 254. 

10 Hierbei sei ein Buch von W.S, Wlafow, Moskau, 1949 in russischer Sprache erwibnt, das u.a. auch 
die Theorie prismatischer Faltwerke behandelt. Die Berechnung wird etwa in der von Hartenbach" 
vorgeschlagenen Form durchgefihrt. 

11 Hartenbach, EMPA-Ber. 99 (1945) S. 162. 
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fallen; fiir die Lésung des Gelenksystems des Faltwerkes gilt das aber nicht. Eine strenge Er- 
fillung der Vertraglichkeit der Verschiebungen lings der Kanten ist nur durch die gleichzeitige 
Wirkung simtlicher in den Schnitten vorhandenen Schnittkrafte gegeben. Damit ist die 
Theorie aufgezeigt. 

2. Geometrie. Abb. 2 zeigt den Querschnitt einer Faltwerkreihe. Die Kanten werden mit 


groBen Buchstaben J, K bezeichnet. Die zwischen den Kanten in verschiedenen Ebenen 
liegenden Platten der Dicke h; erhalten arabische Ziffern 1, 2,3 ....1. Der Abstand zweier 


Abb. 3. Anordnung der Koordinaten- 
Abb, 2. Bezeichnungen zur Faltwerktheorie. eysteme (xp Xp “)- 


Kanten JK ist die Breite a; der Platte. Ihre Lange ist durch den Abstand L der aussteifenden 
Binder bestimmt. Die Winkel, welche die Plattennormalen mit der Vertikalen einschlieBen, 
heiBen w;. Die Winkelhalbierenden der zusammenstoBenden Platten bilden mit der Verti- 
kalen die Winkel yj, yx. An der 
Kante J bestimmen die beiden 
anliegenden Platten mit den 
Schnittufern i und k den Winkel 
2yy. Jedeim Raum liegende Platte 
I wird ebenso wie ein Zylinder- 
schalenteil von jedem Ufer i und 
1 k durch ein rechtsgewundenes 
Koordinatensystem %;, y;, 2; nach 


Abb. 3 beschrieben. 


3. Allgemeine Liésung. Zur Be- 
rechnung dieser Tragwerke wer- 
den die durch eine beliebige Bela- 
stung in den Platten (Scheiben) 
und Zylinderschalen auftretenden 
Spannungen in iblicher Weise 
zu SchnittgréBen (Langskrafte, 
Schubkrafte, Biegemomente, Drill- 
momente) zusammengefaBt. Wird 

A das in Abb. 1 dargestellte Trag- 

; werk an den Verbindungsstellen J 

aufgeschnitten, so treten an jedem 

Jes. TAs be ERE, 20S oN GS SRR ran et om Ufer vier SchnittgréBen als AuBere 

Pe ry elimition der erbin pe anor un er erschie ungen In den Krifte heraus. Diese sind: Quer- 

kraft + Drillmoment, Langskraft, 

Schubkraft und Biegemoment. Wenn es nun gelingt, diese 4uBeren Krafte so zu bestimmen, 

daS fiir die Ufer des gefithrten Schnittes die geometrische Vertraglichkeit der Verschiebungen 

besteht, ist die Aufgabe gelést. Zu diesem Zwecke werden an den Schnittufern i und k wir- 

kende Krafte zu vier statisch unbestimmten GréBen Xj,, Xj,, Xj, Xz, zusammengefaBt. 
Gleichzeitig wird der Verschiebungszustand 4y,, dj,, Oy, und Oj, definiert (Abb. 4). 

Die in Abb. 4 definierten Uberzahligen sind Funktionen von x; Das gilt auch fiir den Ver- 

schiebungszustand der Schnittufer. Da es nun méglich ist, die Verschiebungen der Ufer i und k 
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A 


I. Gleichungssystem 


B 


Cc 


KA XA XA XA, 


Xp.Xp,Xp.Xp, 


X¢,Xe,%e,Xc, 


XmiXmu,Xu,Xm, 


il 


I 
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langs der Kanten unter Beriicksichtigung der an den Binderscheiben vorgeschriebenen Rand- 
bedingungen in Ubereinstimmung zu bringen, ist die Lésung gegeben. Sie‘wird auf die Berech- 
nung eines linearen Gleichungssystems zuriickgefihrt. Die vier iiberzdhligen Schnittkrafte 
eines Punktes J rufen an den anliegenden Schnittufern k und i der benachbarten Kanten H und 
K Verschiebungen 0x,, Oxy; On On, bzw. OK, Ox, Oxy Ox, hervor, so daf damit zur Be- 
rechnung ein zwilfgliedriges System von linearen Gleichungen aufzulésen ist, was bei den 
heutigen Hilfsmitteln keine Schwierigkeiten bereitet. Im Falle des Tragwerkes nach Abb. 1 
ist das Gleichungssystem nachstehend angegeben. In der Matrix sind die Doppelstriche gleich- 
bedeutend mit der Summe der Verschiebungen von zwei Schnittufern. Die einfachen Striche 
besagen, daB nur der Anteil eines Ufers 1 oder k zu beriicksichtigen ist. 


4, Vor- und Belastungszahlen. Zur Berechnung der Vor- und Belastungszahlen werden in 
Abb. 5 die Schnittkrafte und die Komponenten des Verschiebungszustandes definiert. 

An den Schnittufern i und k wirken langs der Kanten in Richtung x die in Abb. 5 definierten 

Schnittkrafte M,,, Nyy, Qyz 

ee = Q,, + M,2z2, Ny, sowie die 


Neer 
z — : 
¥ ee > Verschiebungen @,, x, Vy ys Uggs Ux x 
PDC EG ; 
\yg = a) Vorzahlen. Eine Be- 


Sz ‘ lastung 


Ny Uz2 
a Pie ES ee he 8 aad NA 
Qos) Mey ie ae 
Z y A Qy: ae ° Noe a 
My an einem Schnittufer einer Plat-— 
eae te oder einer Zylinderschale be- 
fs wirkt am gleichen Rande 0 Ver- 


schiebungen R.. und am gegen- 


Abb. 5. Positive Definition der Schnittkrafte und Verschiebungen am Platten- ae : a 
Bea ye itesdaciie tales 8 uberliegenden Rande 1 Verschie- 


if 
bungen R.., die nachfolgend 
definiert sind. Dabei sind alle ibrigen an den Randern 0 und | heraustretenden Schnittkrafte 
Null. 


Die an einem Ufer angreifende Belastung bewirkt Verschiebungen 


am gleichen Ufer: am gegeniiberl. Ufer: 
0 ul 
M ys = 1: Ry = Sethe ae Ry, = Sal Bee 1 Bape 
0 1 
Roy = + Ew, Ry = + E v,, 5 
0 1 
IN ope es AN Bees OE Ue Ry = — Evyy, 
a 
Ryo = + Bye,, » Ryo = + E Vens 
a 0 1 (1) 
On =i Ky, = 4 ba. Ris = — BOya, 
) 1 
Reg ee Paes Ry = + E vz, 5 
0 a: 
IN peak oe eee Ro Ev yy 
1 
Rig = +B vex Ry, = + Eve: 
Zwischen den Verschiebungen nach Abb. 4 und Abb. 5 bestehen die Beziehungen 
am Ufer i: am Ufer k: 
6,.=+EF,., 0) = ED,2, 
0, = — Ev,,cosg +E v,, sin ©, 0g = — Evy, cosg + E vz, 8in | 
03 = + BE v,, sing +> Kv, €os @ , Ne — Ev,,sing — E v,, cos , (2) 
f= Eas a pate —— en ere | 
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Die Einheiten der iiberzadhligen GréBen X,; (nach Abb. 4) nacheinander wirkend liefern nach 
Abb. 5 an den Schnittufern die Randbelastungen 


Xz 


1 


a5 


Aus den Gleichungen (2) und (3) ergeben 
laufig die Anteile der Vorzahlen an den Schnittufern i und k. 
sammengestellt. 


Os a 


OJo, 
Oe 


0 
= + Ry sing , 
0 


Ry ? 


+ Ry, cos @ , 


0 


am Uferi: 


= 1 lhiefert Miyg= +1, 
Nyy = 0, 
Qyz = 0, 
IV ix — 0 ° 
, = 1. liefert a 
Nyy = + cos 
Qyz =—sng, 
oe OR 
= 1  liefert M = (Oe 
Nyy = — £10 0, 
Qy, ='— cosy, 
Oey = 0, 
=1 liefert M,,.= 0, 
Ae On 
Qyz a 0, 
Nae ks 


Il. Vorzahlen. 


Die Uberzahligen Xj; am Uferi 
wirkend liefern die Verschiebungen 
Oynj am Ufer i: 


= — Ry,sing, 
0 0 
= — (xe R,. cos? g + Reg sin? 9) > 
/ 0) 0 


— sin y cos p (— Roy + Rss) ’ 
= Re cos@ , 
0 
— Ry cos@ , 
— sin Y cos Pp (IR Rag an ie 


—(+ ze sin? g + ibe cos? 9) ; 
4. ie sing, 
0, 
0 
— R,, cosp , 
0 
+ Ry sing, 
0 
— Ry» 


am Ufer k: 
MM, +1, 
Nyy = 0, 
Oy, =o Os 
IN 4 a8 ie 
My x ee 0, 
1 + cos, 
Oy, =— sing, 
ih Poe = OF i 
3 
eats. (3) 
Nyy = -b sin @, 
Oy: “=a + cosg, 
Die = 0, 
MS er 0, 
Ney = 0, 
On — 0; 
Nye =a 


sich unter Beachtung der Randwerte (1) zwangs- 


Sie sind nachfolgend zu- 


Die Uberzahligen Xj; am Ufer k 
wirkend liefern die Verschiebungen 


OSs = 


OF, 
OF; 
Oh 
Oe 
O Fes 
DJs 
Oy. 


Osn; am Ufer k: 


0 


Ss Ry ? 


0 
+ Rg sing , 
0 


— R,, cosg , 


9 


= Ris sing , 


O13 = 


DJo5 
OF53 
OF. 
Osu 
OJes 
OTe 
OF 


— (R,» cos? @ + Ry sin? 9), 


+ sin @ cos p NE Ry + Rel 
0 

+ Ry cose, 

a= fie COs Pp , 

+ sin p cos p (258 Ry. + Rae 


— (Rg sin? @ + ihe cos? @} , 
a ne sin , 

0, 

0 
+ R,, cos, 

0 
+ Ry sing , 

0 


== Kay: 


10 
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III. Vorzahlen. 


Die Uberzahligen Xj; am Ufer k Die Uberzahligen Xj; am Ufer 1 
wirkend liefern die Verschiebungen wirkend liefern die Verschiebungen 
Ox,j am Ufer i: Ounj am Ufer k: 
7 1 1 
| OKy = + Ry ? On, == -+ ah ° 
: . 
Xy, = le * OK a a Ro sin PY ? Ont =F ao Ro sin 7 ) 
| OK, = + R, COS (p , One: = — Rs; cos@ , 
OK ae 0, Ong == 0, 
1 1 : 
Ox, = — Ry sing , Ong =o Ra a : 
uf 1 . 
xX 1 | OK, == ee Ry. cos? p + Rgg sin? 9) ; On, = = (4 R,, cos? os Ry ae ¢) - 
= 1s 1 1 
ie | OKy = sin @c0s 9 (— Re eRs OH. = + sin g cos y (— Ryo + Rss) ; 
1 1 
OKy Say Ry cos P 2 OH. —— ales Ryo cos 7 > 
1 1 
OK4s a a5 Ry cos Pp ° i ; OHys SS Ry cos YP > : ‘ 
x 1 OK, = + sing cos y =: Ro + Rs) : On, = — ae cos Pp = Re + Ree 
= 1 
: | OK = + (+ R, sin’ y -- Rgs cos® ¢) ’ Og. aay & Roo sin? p + Rss cos? 9) , 
s > . 
OK 4g Sp Ryo sin P > O43 = — Rp sing, 
OKs = o Ons = Ys 
OKo4 = + Ry, cos@p , Ono, = — R,, cos, 
Ay, == II ¢ 1 i 1 ; 
Oky = — Rysing , Ou = — Rog sing , 
1 
OKs = = Ry ? Ong a ae Ry, . 


b) Belastungszahlen. Der Verschiebungszustand infolge einer Flachenbelastung 
wird in dem Koordinatensystem am Ufer i beschrieben. Entsprechend den Gleichungen (1) 
0 al 
werden die Randwerte R.. und R.. fir die Einheitsbelastungen der Komponenten Poo und 
pzz definiert. 
Die Flachenbelastung bewirkt an den Ufern i und k Verschiebungen: 


Belastung Wirkung am Ufert: Wirkung am Ufer k: 

0 1 
Pyy =1: Roo =+Evy,, Roo =—Evyy, 

Rig = Eves 5 Ra = se One 

0 1 (4) 
Dee Ry = +E Oyx ; Ry = — E,., 

0 1 

Reg o£ oe Rog = = Bw yy: 


Damit kénnen die Belastungszahlen fiir die Schnittufer 1 und k des Gleichungssystems be- 
rechnet werden. Unter Beachtung von (2) ergibt sich 


am Uferi: am Uferk: 
i 0 ; ’ 1 
Ofer Fi98 On, = + Ry» 
/ 0 0 " 1 ue 5 
9 Jog sis (= $20 cos p + Reso sin g) ’ Oty = + i= Ryo cos y + Ry sin ¢) ’ (5) 
\ . - 7 1 
9 Tap = le Ry sin p + Ry cos @) , Cine a a Ry sing + Roy cos ) A 
1 
One a Ras Py OH —— Ryo 


Mit der Ermittlung der Vor- und Belastungszahlen des Gleichungssystems ist die Lésung des 
zusammengesetzten Tragwerks in allgemeiner Form gegeben. 


5. Symmetriebetrachtungen. Bei der Aufstellung der Matrix und zur Abkiirzung der Berech- 
nung sind Symmetriebetrachtungen von Nutzen. Ist der Spannungszustand aus einer Flichen- 
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belastung zur Achse y = + a/2 symmetrisch, dann bestehen, wie man aus der Abb. 4 ablesen 


3 0 nt 
kann zwischen den Randwerten R.. und R.. die Beziehungen 


0 1 
+ Rp = + Ry, | 
1 
R 


=i. ra R, o) 
aera (6) 
sig Rs seals Roo ’ 
+ Ry = + Ry - 
Im Falle eines zur Mittelfliche antimetrischen Spannungszustandes gilt entsprechend 
0 1 
iis Ryo oe Ryo ’ 
Rey == == Rog 
ti *20 $20 (7) 
2% Roo ae Reo ’ | 
sfoathag age 


Fir eine im Raume liegende Platte mit beliebig geneigter Achse y kann der Dehnungsspan- 
nungszustand antimetrisch, der Biegespannungszustand symmetrisch eintreten. Dann gelten 
fiir die Beziechungen am Ufer k in Gleichung (5) unter Beachtung der Gleichungen (6) und (7) 
auch die Ausdriicke 


g } 
On, = + Rio» F 
07.6 = (i Ry cos p + Roo sin 9) : (8) 
0a i Re sin p EE R., dos)'s | 
On, — Ry. J 


Betrachtet man ein im Querschnitt symmetrisches Tragwerk, dessen Symmetrieachse in der 
Mitte einer Platte oder einer Zylinderschale liegt, so gelten fiir die an den Schnittufern k und i 
der Kanten J und K angreifenden Uberzahligen nach Abb. 4 die folgenden Gleichungen: 


bei symmetrischer Belastung 


KK = a AS,» 
Xk, — =e EN (9) 
Ane = XG,, 
Ag, = —Xj,> 
bei antimetrischer Belastung: 
Xx, nae Xj; | 
Staab (10) 
AK, oe =e Mas | 
Nie ey, 


Fallt die Symmetrieachse des Tragwerkes mit einer Kante J zusammen, dann werden die vier 
Uberzahligen Xj am Ufer so bestimmt, daf fiir symmetrische Belastung die Bedingungen 


Os, ee 0 2 | 
Digan? (11) 
er ay | 
oy Al) 
und fiir antimetrische Belastung die Gleichungen 
Ge => 0 9 | 
nO ges (12) 
6;,°= 0; 
Oy = 0 


erfiillt sind. 
10* 
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% 


Die strenge Theorie des in Abb. 1 dargestellten Tragwerkes ist somit auf die Berechnung 


0 
der Randverschiebungen R.. und R .. fir Einheitslasten 
My gale, IN eee bes Opa sh Hh eral 
und fiir die Flachenbelastung 
Pyy = 9 Pzz = 1 


zuriickgefiihrt worden. Zur Ermittlung dieser Verschiebungen stehen die bekannten Theorien 
der anisotropen Scheiben, Platten und Zylinderschalen sowie die Theorie der Stabe zur Ver- 


fiigung. 
Fir die hier am meisten interessierenden isotropen Scheiben und Platten soll die Berech- 
nung der Randverschiebungen durchgefiihrt werden, da sich in diesem Falle die Rechnung 


) 1 
allgemein bis zu expliziten Werten R.. und R.. durchfihren abt. 


6. Das isotrope Faltwerk zwischen Endbindern, Zur Berechnung der isotropen Faltwerke 
stehen fir die Dehnungs- und fiir die Biegespannungen die Differentialgleichungen 


— y2 
AAF=0,  Adv,, =F”? p(x y) 
| ee ue (13) 
mit ~ Ox 1 Gy? ps ey | 


zur Verfiigung. Die in Abb. 5 definierten Schnittgré8en und Verschiebungen ergeben sich aus 
der Spannungsfunktion F und aus der Normalverschiebung v,, in der bekannten Weise: 


i ae M Eves 
xx = + Fyy, “y= "Typ pa (TF Petes $Y Mariya) » 
E hs 
Na = + Tatas Me eS 12 (1 — 2?) ( V Vz x0 x a Cee 
EhS 
Nay = Nyx = — Frey — | pyy(x) dx , Moe = Mop 12 (1 — ») (1 =?) trax ys 
Eh? 
E S45, SU ain =F ey =p ep 2 OF: = ee 12 (1 — 2?) Av zai 2 
(14) 
E h8 
Eheéyy = yyy = + Fixe — Fryy, Qy; 12 (1 — v2) Avezy » 
Eh | ate Q ae ( rv 
Vay — ar Vealy a Vy y\x dh NOT = y?) peice 
= —2(l a I See + (2 =») tesizyy) + 
a ihe 
we Se 12 (1 — ¥?) ( peta © —yv) 


E0,,.=+Eu, 


Der Verlauf der Belastungskomponenten in der x-Richtung wird durch unendliche Reihen 


dargestellt: 
— ~ Ne x “ = Tt 3 
Oy Digs > Pyn COS] — he = > Pen cos 7 oe (15) 


Als Lésungen der Differentialgleichungen (13) eignen sich im Sinne von (15) die Funktionen 


oo 


F =f) cos = ate ee = > w,(y) cos n= : (16) 


w= 0 n=0 


Fir eine zweckmaBige Durchfiihrung der Lésung werden die folgenden Abkiirzungen und 
dimensionslosen Koordinaten eingefihrt. 

_ nna A x A Rife A 2 

A= re =) ee p = 1h, CD == Ie % B 


cs @ a a ee 


(17) 


Mit den Gleichungen (16) u. (17) ergeben sich fiir die Schnittkrafte (14) die nachfolgenden Be- 
ziehungen (Striche bedeuten Abteitungen nach 7, die Summenzeichen sind nicht geschrieben); 
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Eh’ A fy ) 
Ma is (ea yw-+o'’)cos€é , 
Vea Ty, Eh? A my) 
1s oie maces cos é, M,, = 4 Ra) a w—vw'')cosé, 
- A Eh A 
Nyy =—— pcos, Dig 12) amie sin € , 
Aen a : Eh’ 7 he 
1 Se = ee, i) ot 
Bye | (4 ~# 7 Pyn) sin €, Or; 2 >) = 1G @ — w’’) sin &, (18) 
A wr Eh’ Me , mr 
Eheg=+—(+ +) cos €, Qys =e ge + w )cos&, 
A W 5 Eh A 
Ehéyy=— —(+ +9 ) cos €, Que = Fatih at w — a’ (2—»)) 
sin sae 
$ Eh? ie ; 
ee ay 
+"""| cos &, 
E?=-+ Eo’. J 


Fur die Funktionen gy und w@ ergeben sich wegen der Doppelwurzeln der charakteristischen 
Gleichungen die Lésungen 


p= +2 fb em [At Aa (n+ 0) + em" [4o + Aa(n + 0] | 
a 12 (1 — 12?) 5 (19) 
om i eB + Bala t+ 0) +e" [Bs + Bulg + OD) «| 


Setzt man die Funktionen (19) und ihre Ableitungen in die Gleichungen (18) ein, so sind die 
Schnittkrafte und Verschiebungen als Funktion der Koordinaten &, 7 und der Integrations- 
konstanten A,, B; (i = 1... 4) gegeben: 


De a 1 et" (4 A, + 4, (9 + 2)) + e-7[+ A, + A, @— 2)]} cos €, } 
BN hee {+ etn[+ A, + A,-n] + e-"[+4+ A; + A, 7} cos é, 

eee ee fa A a (1) ee t= . Ay (= 4 4)]} sine 

M,, = — {+e+7[4+ B, (1 —») + B,(+7 (1 Aiea ard [lero 2 le) 


oe (+7 (1 —v) — 2)]} cos, 
= + {+et7[+ B, (1 —») + B(+7(1 vy) —2y)| + e—-"[+ B; (I —») 
: 3 + By (+7 (L =») + 29)\]} cos €, 
M,. = +(1 —»){ + e+"[+ B, + B(+9 4+ 1)] 
te" [— B, + B,(—9 + ]} sing, 


s 
| 


Ore = +24 {4 e+" B, —e-” By} sing, 


2A = : 
0. = -—— te eum. Bel tea By cose, ’ 


Q.. = +4 {+ et [+ Bl ») + B, (4 
Weedon (1 —y) + 
Ose =4+A{+erm[+ B (1 v) + B, (4 | 
—e"[+ B, (1 —») + By(+7 (1 —») + (1 +2))]} cosé, 
E ve= +p5tt em (t+ AL +») + 4a{t 0 +» + 3] 
4 e-"[-+ As (1+ ») + Ay f+ 4 (1 + ») — If sine, 
{4 etn[+ A, (1 +7) bk {+n(1 +7) —(1 —»)}] 
+ e-"[— Ay (1 +) + 4a {— (1 +7) — (1—»)}]} cose, 
Ma eis foam) wwe 


(+ n(l—yv) —2 (2 v))]t sin & , 


a 
Evyy = — =: 


Ey. = + Go {tet By + Baly +] +e "[— By + Bi(—n + I} cose. | 
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Die Bedingungen an den Randern 0 und 1 der Scheiben bzw. Platten sind im Sinne von (1) 


festgelegt: 


1 = 9 Nyy =1, Nyx =0, 
Pash Ney = 0, Nyx =9, 
2.7, == 0 Nyy =9, Nyz=1, 
n=A Nene IN == W% 
1) Sy oo (21) 
3.4 = 0 My, =1, Qy, = 0, 
i =A My, =, OF = 
4.7 =0 Myx =9, Qy,=1, 
fs M,, =0. Q,,=0. J 


Zur Erfillung der Randbedingungen 1. und 2. in (21) ergibt sich fiir die Integrationskonstanten 
A; (i =1..4) das Gleichungssystem: 


A, A, A, A, 

+ 1 0 | +1 0 == Nyy 
si cat es | +1 + Nyx 
ets +hett +e hens 0 
af AS ba A) Oe, Coke Realtor : 


Die Forderungen 3. und 4. in (21) bedingen das nachstehende Gleichungssystem zur Berechnung 
der Integrationskonstanten B; (it =1..4): 


B, B, B; B, 
+ (1 —») Oy: +(1—y) =P 3) Bea 
Tee?) = (ler 9) =e) Saga) +7 Oy: 
+ (1 —y)et4 | + [A (1 — v) + 2] et4 | + (1 — v) e*| + [A (1 —v) — 2] ee 0 
(lop) oF * | FTA 9) (1 foetal = pet a[Aitony) (2) cet ae 


Die Auflésung der Gleichungssysteme liefert fiir die Bedingungen (21) die folgenden Lésungen 
(22) fiir die unbekannten Konstanten: 


lL. N,, = 1: 2. Nye = 1: 
A =4 ca 2/2 + a4), A, = Fae 
A, = ene 2=+77h. 
A, = 4 roe 2 A? OX») Ay= +7528, 
A= —7am, Ae cats 
3. My, =1: 7D Bi jes r (22) 
Ese Tay ta); Bal a ear Ga; =P + Ba) 
ie CSIy oN Oite ENNELY p ay 
isa Nene (- 3 1344) ? B, +a7=wal tbs eae 
B, = Sa B, +a na, 
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Die Abkiirzungen in (22) haben die folgende Bedeutung: 


oj = +Aet4 + Sind, Bi = tAet+4 — Sind, 
=+Ae7* + Gind, Bo = tAe+ — Gind, 
a= +s" Sind —Act4, fy = +2 =" Sina tA cr, 
(23) 
a= +722 SinA—de*, = fy +E SindA te, 
A, = + Gint'a —2. da= + (75) Sine a — 78. 


Mit den festgelegten Konstanten kénnen nun aus (20) die Gleichungen der Schnittkrafte und 
Verschiebungen fiir die Randbelastungszustande 


Ney = 1, Noga lL, My,=1, Qy, = 1 


berechnet werden: 


Noy == 1: 
Nyy = = tee Perth ay (1 = 9) | | 
+ e+@— [+ A? e—* — a (1 + 7)]} cos &, 
Nex = Tey {+ e@—) [4 A? e+4 — a, (1 +-)] 
+ et+@—n) [+ 22 e—4 + ay (1 — 9)] ]} cos €, 
Nyz=+ +1 {+ e—@—) [+ 22 e+4 — a] + e+?) [22 e—* + ay n]} cosé, 
ui 
1 2 
poy [eoeeLtre tn (eta ° 


ee? |p Ae + a, (T35+7)]} eos 


) 


+ e+@—n) |+ Re — Oy ( rH : Slik sin ¢ ; 


Eee ag te te —% (+9 +75 


Nowe = 55 {+ e441) [— J? et4 + Bin] + et 0) [Ae — By n]} cos & , 


Neo = + yiy (het) [Be 48.2 +0] 
1 + e+@—n [+ 22 e 4+ B, (2 — )]} cos é, 


Nyx = att a Anette By (lt) 
+ et @—1) [— A? ee — B, (1 —7)]} cosé, 
a (1 =») ia} de Ss) 1—y (25) 
LE ae SE {4 nay |-4 : a Pale seg) 
+ e+@—1) ie Ae? + Bo ( a + n)|| cosé, 
E vs. = 4 aT ieee tant ooh Gare Te. n) 


——— 


+ e+@—1) E: Ve pa(5 a —n)| sues 
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Moot 
Myx = cat me berets ea (A — 7) “ape 
a ue bt ag +9 {2 (2 — 1) esi cos, 
May = — ga ft er@ +7 rag + a4 f24 — 9) — FF 
+ mes Fag + 9542 (A — 1) +- rT} hoo e 


Pi 
ES Sapper 
Mie ae ( eer 


canine 


“3 + a4 {2 (A — ») +1}] 
dg dig {2A 2) —Uljsing, 


=| 


qh 


a 


+ med pais em 
OF te het oe joes - Os +4 {2 (A —n) — 
beta [fet 
Caer 2 {+ cee aes oa = Og + a, {2 Aim) 3 
+ ectan[ — Tag + 05 {2 (A =) — 
== ii {+ et+G@—) a, + e—@— ag} cos &, 
HA) oy — e—U—") ag} sin E , 


Cree rie I —») A, {+e 


on ~ My {2 (A — 7) +p F}| coe & 


(26) 


3 


eae 


a 


E U2 — ahBAAT oe ae + O4 {2 ( (A 1) ) 1}| 
ae rete ie + a {2 (A —») + 1} eos en 
| Lie | —1 E 
E By, = 4 AUER. ent M+ Tra +34 {2 (4 — ») +y 
4. e—(4—n) |+ : Late — Os {2 (A — ”") — 1} cos & ; | 
Q,.,=1: 
sematalteHEtEa nba 0 2g | 
ae 17 Bs + Bs {2 (A — 0) — ya 5}|feosg. 
a 3 
Mey = —qrq{tere—)[4 ita — B20 — 9 Ey 
$00] TSBs + AaB 0) + Ea cose, 
eocame Hef een = pat (2 (A — 9) +3] 
The eet ee 2G —1) — Hsing, 
Q, ewes oo — Bf? (de) = 


see BO acres 


Ba + Ba {2 (A = 0) + 


a tH cosé, 
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cua eee Be sees) —Bf2 (A ee 


b= | (27) 
ee + 7—3>y])_. 
jer =1) ee I = Hi se Pa {2 ( (A — 7) “TT heme . 
Oy. = ea) Tage pCa, oe —Eaa) Bet eos, 
0. = — cere {+ e+@—n) B, + e—@—”) B5} sin E , 
(1+ y)a pee 3+ 4 : 
Ce aly A: f Lert »| peda - 1] 
+ ¢—“— JF fa By + Bs {2 ( (A —n) + I} eos & 
- (1 + ») LN Bay \ 
oh a ieaa, {4 ett D ehaceaet Bee Hy Aa i 
=e-U-) fre: pci” Ba + Bs {2 (A — 7) — HI cosé, 
4 J 
Fir die Flachenbelastung mit den Komponenten 
Pyy = Pyn cos é, Pzz = Pain cos & 
gelten die Randbedingungen 
"= 0: : 2 
i ee Nyy = Nyx = Myz = Qys = 0. (28) 
Fir die Differentialgleichung (13b) wird noch die Partikularlésung 
ees a? (1 — >?) _ 
Eh t.s=-+- psn ce cos é mit He (29) 


 benétigt. 


Die Randbedingungen (28) werden erfillt, wenn die Konstanten A; und B; (i = 1..4) den 
folgenden Gleichungssystemen (30) geniigen: 


+(1—y)et4; + [Al —v) + 2] e tA | + (1 — pr) tfA(l —v) —2]e7 | — vy pin _ 
+(1—v)e+4] + [A(A—v) —(1+9)]et+4) — (1—v)e | — [ACL —») + (1-+-r)]e 0 


eae a ees mr } 
= 0 fad 0 0 | 
Se a —1 | sal +4Pyp 
+ et? +het +e +he- 0 
+ et +(1+A) —e ce (esd) +7 Pyn 
(30) 
B, B, B, B, | 
(1 — ») +2 Ail — 9) at +9 Pen Fp 
+ (1 —») (le) a=) (ey) 0 
J 


Die Auflésung der Gleichungssysteme ergibt nach langerer Rechnung die Werte der Integra- 


tionskonstanten 
a 
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BS SHS Tae ae 
y a 
Ay =) Pong 7 yet? 
a 
A, SS Pyn 24,4"? 9 
B es Da i a? 1 =" Vv yp) 3] | 
t oa Pie 992 (229) d,\1 = 9? : (oa | 
a? 
B, = Pe Pen Cy tear 
mM: a? Ue cae 
a2 
B, = -|- V Pin 2rd —») dA, V2° 
Dabei sind die folgenden Abkiirzungen definiert worden: 
41 =+ Les 
pea ae l= Cares, 
A, = + Gind —A, (32) 
3 Dems 
Ay = +7=" Gina —A. | 


Integral der Plattengleichung (29) die Schnittkrafte fiir die Flachenbelastung berechnet werden: 


Aus den Gleichungen (20) kénnen mit den Integrationskonstanten (31) und dem partikularen : 
s } 
Nyy = Pya pq, em (—A+ ny] te"[+4 + yon ]} cosé, 
NE + Pyngzq {+ e[—A + 1 (2 +7] +e" [+4 +72 (2 — mf cond, 
1 ; : 
Nye = + Pong | tga. lt ett {—4 + (1 +0} 
+e {—A— 79 (1 — 9) — Hsin g, 
7 = (Lise aa? 
| SOPs = Toe re 3 A —- Mil (n = I ate al 
eas f AG I ce )}peoe 
Weg Sm yo aay at | aes n(j aL +n) | 
2 z 
et Ae alters nh} sin € , 
2 1 / 
Myx = +9 Pen es It te etl 4 n (4 n) a ) 
3 v 
= oi {a Yo Gi as : — n) == a cos &, 
Mey = + Pont l tyy (tem {+n +2) — 3 
—e—{+ y, (1 —n) +3}} cos € , 
=, J a’ [ TY, ) 2 7 
Min Pen we Tn 1 @ | Ly, (+ Oe ss — "| = i (33) 
bem| + 79 bra) + allsing, 
a 
Oye SS Ue Fie ae ert V1 + e—1 Yo} cosé ; 
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0... = 4 Y Pen 5 ie) a {+ ety, — ey} sin E , 


Qy: = +9 Pin grag {+ ern (+ y,9 —A) —e—" (+ yon —A)\cosé, 


4) 

—e-n 5 es = n) 4 A)| sin § ; 
a ar r TOeER pale fear ez es +n} 
2, tt fiers ate nhl eos 


0 = eee +4 5 
On +9 Pen syalte slaarels 


E v,, 


I 
+ 
fe 
3 


1 ) 
E852 = +9 Pony pqeed, {+e “is A+n( 


okay ee 


Damit lassen sich nun die Randwerte R..und R.. aus den Gleichungen (24), (25), (26), (27) 
und (33) unter Beachtung ihrer Definitionen (1) und (4) berechnen: 


0 1 

Rog = — pz q (+ Sin 24 +24), Ry = + gaz (+ AGO} A + Gin), 

0 — . ) 1 : 

Pe ==) = ae ala Sinta +), Ry = — pa Sind, 

> Le 9) 1+ 1 2 : 

Ry = oz (+ Sin?A +25 Tse Ry = — py Sind, 

0 1 a a . 

n= — gz, + (Gin 2a — 2A), lige ma ACoj A — Sind), 

oa 1 (1 + »)/? eee; aap 

Ray hi?x A, (1 —») Ba Jaa, re Sind, 

3+ A. 1l—vpy 
ee, OU ef) 

0 1 l+y» 1 f= ieee 

arrears ua hiax4,1—y 

x (+ p= iein 22 —24), «(+ 2 Sind — 26044), 

w a 1l+yp ae 2a 1+ y» 
Ry hax A, 1—y UE hi®x A, 1 —v , (34) 

x(+ $e ein 2a +22), x (4 SEteina 4.266 a). 
9 1 (1+) 1 Beli iia 
pis = F iq» A, (1 — ») Rae Ahz A, Teen 

3+ %% l1—yp 

x (+222 ina 1289), 

0 na 1 pyn @ 

Boy = +5 9¢4,(1 —”) , Ro =— Fag (l—?) Ky 
ae x (+ Gina ere 
l—vpyp 

A n2a 

ss — Pre (+ Ooh masts Ryy = + B25" (+ CofA — 1), 

0 2(1 — v)? 1 ss a? (1 — y?) 

Bt es Rg = + Pex hitx 

: EE NCHA (a a ising — A), 
vs f cereus ai Sells a) d “ Bike AR 


0 9 1 1+ y 2a J. 
Ry = — pes? eer, OA —- 1), mnt es Ae —vhiu Ay (CofA Le 
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7. Sonderfall verschwindender Kantenwinkel. Die vorliegende strenge Theorie der pris- 
matischen Faltwerke gilt auch fiir den Sonderfall mit verschwindenden Kantenwinkeln Py 
Das Tragsystem zerfallt dann je nach Art der auBeren Belastung in ein Plattenwerk bzw. in 
ein Scheibenwerk. 

Im Falle einer Belastung senkrecht zur Ebene ergibt sich ein System von durchlaufenden 
Platten, die in Abstanden durch elastische Trager unterstiitzt sind (Abb. 6). Damit ist die 


Z 
Abb. 6. Querschnitt einer iiber elastischen Tragern durchlaufenden Platte. 


Frage der mitwirkenden Plattenbreite von Tragern, die mit kreuzweise bewehrten Platten 
monolithisch verbunden sind, bei Beriicksichtigung ihrer Verdrehungssteifigkeit eindeutig 
geklart. Die mitwirkende Breite wird im wesentlichen von dem Verschiebungszustand der 
Platten zwischen den Tragern ab- 
hangen. Die Lésung bei beliebiger 
senkrechter Belastung ist durch ein 
lineares Gleichungssystem fiir die 
unbekannten AnschluBkrafte  ge- 
geben. 


Z Abb. 8. Querschnitt einer Faltwerkreihe. 


, 


ys 
SS 
S 


Mery [mt] 


7 
Yea gene Mey (rn) 


Abb. 7. Querschnitt einer anisotropen Scheibe. 


Im Falle einer Belastung in 
Richtung der Ebene entsteht das 
Tragsystem einer anisotropen Schei- 
be (Abb.7). Die verschiedenen Quer- 
schnittswerte der Einzelteile finden 
exakt Beriicksichtigung. Zur Be- 
rechnung werden die Tragglieder 
durch Langsschnitte voneinander 
getrennt und die geometrischen 
Vertraglichkeitsbedingungen der 
Schnittufer durch die Wirkung iiber- 
zahliger GréBen erfiillt. 


r yo Uae Ny y des Faltwerkes im 
8. Anwendungen. Die dargelegte Miteclquevsebenee 0 
s ‘ ———— Schnittkrafte im Faltwerk; 

Theorie der amsoriope™ Faltwerke eae Schnittkrifte der dem Faltwerk einbeschriebenen Zylinderschale. 
hat gegeniiber den bisherigen Be- 
rechnungsverfahren nicht nur hinsichtlich ihrer allgemeineren Giltigkeit, sondern auch wegen 
der Ubersichtlichkeit des Rechnungsganges manchen Vorteil. 

In Abb. 8 ist der Mittelquerschnitt einer von den iblichen Abmessungen abweichenden 


Faltwerkreihe zwischen zwei Binderscheiben dargestellt. Alle Platten werden durch die Form- 


Abb. 9. Schnittkraftverlauf My x M,;, 
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gréBe 1 = nz a/L = 2 beschrieben, so da fiir das erste Lastglied der Reihen (415). ant n= 1 
und der Lange L = z ein ,,kurzes Faltwerk“ (analog der kurzen Zylinderschale) oder fiir das 
n-te Lastglied ein ,,quasikurzes*¢ Faltwerk der Lange L = na vorgegeben ist. Die Abb. 9 
zeigt fiir die konstante Flachenbelastung von 0,310 t/m2 jeder Platte bei einer Randtragerlast 
von 0,130 t/m den Verlauf der beiden 
Biegemomente, der Normalverschie- 
bung und der Langskraft N, , im Quer- 
schnitt. Zum Vergleich der berech- 
neten GréSen sind in der Scheitelzone 


O17 
~0157 


die entsprechenden Werte der Parti- Mey (nt] 
kularlésung der dem Faltwerk ein- 
beschriebenen isotropen Kreiszylinder- 
schale mit dem Radius 14,3 m ein- 
getragen worden. Es zeigt sich er- S S 
wartungsgem4B eine gute Ubereinstim- 3 = $ 
mung der Ergebnisse. Fiir das Faltwerk S$ 
ist zu bemerken, da8B die Verschiebun- F 
gen v,, und damit auch die Biege- 
momente M,, zwischen den Kanten W,,[t] 
YY 
Nutelast 0500+ /m x 
~ 
2,00™ 200% 
Abb. 10. Abmessungen der durch Abb. 11. Schnittkraftverlauf M oo Mt. y? Nex Nyy im Mittelquer- 
Randtriager ausgesteiften Platte. schnitt der Platte des Tragsystems nach Abb. 10. 


gréBere Werte annehmen als bei der Zylinderschale, da die gesamte Flachenbelastung der 
Platten durch Biegung abgetragen werden muf. 

In der Abb. 10 wurde eine der das Faltwerk erzeugenden Platten mit zwei Randtragern 
biegesteif verbunden. Im Schnitt zwischen Randtrager und Platte greifen die vier Schnitt- 


eroBen M,,, N,,, Q,y: und Ny, an, die die geometrischen Vertraglichkeitsbedingungen der 
Schnittufer erfillen. AuBer dem Eigengewicht wirkt die 


Nutzbelastung 0,500 t/m? (Abb. 10). Der Verlauf der Schnitt- thn! 


krafte M,,, M.,, Nex und N,, ist fiir das erste Lastglied yee “eS Sasehnity. 
durch Abb. 11 dargestellt. Die Randtrager werden auf & 4 
schiefe Biegung mit Langskraft und Torsion beansprucht. 
Die Langsspannungsverteilung zeigt Abb. 12. 

9. Zusammenfassung. Zur Berechnung der Spannungen mon + Mot Jae 
und Verschiebungen in dem ausStaben, PlattenundSchalen- — 4}), 12. Langsspannungen ox zim Randtrager 
elementen zusammengesetzten Faltwerk werden die Kinzel- des Tragwerkes nach Abb.10. 


teile durch Langsschnitte in den Verbindungsstellen ge- 

trennt. Die in den Schnitten wirkenden Verbindungskrafte erfiillen die geometrischen Ver- 
traglichkeitsbedingungen der Formanderungen der Schnittufer, wenn diese einem linearen 
Gleichungssystem geniigen. Die Ableitung zeigt, daf’ in jedem Falle die Aufstellung der 
Gleichungen fiir beliebige Bedingungen an den Randern + L/2 des Faltwerkes méglich ist. 
Die Vorzahlen der Matrix sind neben den trigonometrischen Werten der Schnittwinkel Funk- 


tionen von Verschiebungsgréfen R.. fir Randbelastungszustinde der Faltwerkelemente. 
Diese kiénnen fiir beliebige orthogonale Anisotropie nach den bekannten Theorien der Stabe, 
Scheiben, Platten und Zylinderschalen berechnet werden. Die Auflésung der Bedingungs- 
gleichungen liefert die Verbindungskrafte in den Kanten, welche mit der Lésung fiir die 
Flachenbelastung den Spannungs- und Verschiebungszustand in jedem Punkt der Elemente 
bestimmen. 

Die vorliegende Theorie ist auch auf den Sonderfall verschwindender Kantenschnittwinkel 
anwendbar. Im Falle einer senkrecht zur Plattenebene wirkenden duBeren Belastung entsteht 
ein System mehrfach zusammenhangender Platten, die mit elastischen Tragern verbunden 
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sind. Damit kann auch die Frage der mitwirkenden Plattenbreite von Tragern unter kreuz- 
weise bewehrten Platten beantwortet werden. Wirkt die 4uRere Belastung in der Plattenebene, 
entsteht das Tragsystem einer orthogonal anisotropen Scheibe. 

Fiir die zunachst am meisten interessierenden isotropen Scheiben und Platten sind die Glei- 
chungen der Schnittkrifte und Verschiebungen fiir harmonische Randbelastungszustande und 
Flachenkrafte explizit angegeben. Der Verlauf innerhalb der Platten und Scheiben, beschrie- 
ben in dimensionslosen Koordinaten, wird im wesentlichen durch einen Parameter / bestimmt. 

Fiir ein isotropes Faltwerk zwischen Endbindern wurde gezeigt, daB der Spannungs- und 
Verschiebungszustand mit dem der einbeschriebenen Zylinderschale erwartungsgemaB gut 
iibereinstimmt. AuBerdem wurden die Schnittkrafte und Verschiebungen einer mit zwei ela- 
stischen Tragern monolythisch verbundenen Platte berechnet. Zu dem Biegespannungszu- 
stand der Platte addiert sich ein Dehnungs-Spannungszustand. 


ai (Eingegangen am 13. Juni 1954.) 
Anschrift des Verfassers: Dr.-Ing. D. Riidiger, Freiberg/Sa., Wernerplatz 9. 
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von Saint Venat-Mises und von Prandt]-ReuS bei ihrer Anwendung auf das Torsionsproblem. Kombinierte 
Torsions- und Zugbeanspruchung eines Kreiszylinders. — Ebener Verzerrungszustand: Probleme mit 
axialer Symmetrie. Allgemeine Bezichungen. Inkompressibles Material. Entlastung und wiederholte  - if 
Belastung. Uneingeschranktes plastisches FlieSen. — Ebener Verzerrungszustand: Allgemeine Theorie. 4 
Einleitung. Allgemeine Begriffe. Geometrische Eigenschaften der Gleitlinien. Randbedingungen; = |_ 
Naherungsweise Konstruktion von Gleitlinien. Geschwindigkeitsfelder. Scharen gerader Gleitlinien. . | 
Grenzlinien. Unstetigkeitslinien. — Ebener Verzerrungszustand: Besondere Probleme. Probleme des be- | 
ginnenden, stationaren und pseudostationaren plastischen FlieBens. — Ebener Verzerrungszustand: Ein- 
geschrinkte plastische Verformung. Traglastverfahren. Die Plattenanalogie. Analytische Behandlung oe | 
eines speziellen Falles. Das Prinzip der virtuellen Arbeiten bei ebener Verzerrung. Traglastverfahren bei =f 
ebener Verzerrung. ‘Lésungen nach der Saint Venant-Misesschen Theorie und das beginnende plastische 1 
FlieBen eines Prandtl-ReuBschen Materials. — Extremalprinzipe. Einfihrang. Cartesische Tensoren- Ie 
Extremalprinzipe der Misesschen und der Prandtl-ReuBschen Theorie. Traglastverfahren. Einige Be- i 
merkungen betreffend die Anwendung von Extremalprinzipien. — Namen- und Sachverzeichnis. Auf- : 
gaben und Literatur am Ende jedes Kapitels, Verzeichnis der im Text angefihrten eng- 
lischen Fachausdricke. Jey i 
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Aufgaben aus der Technischen Mechanik a? 
Graphische Statik, Festigkeitslehre, Dynamik fester Kérper. Von Dr.-Ing. Rudolf Sonntag, 0. Professor sf 
an der Technischen Hochschule Karlsruhe. Mit 324 Abbildungen. XI, 209 Seiten Gr.-8°. 1955. | 
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er Studierende der Ingenieur- Wissenschaften nt oft vor Schwierigkeiten, wenn es gilt, die verhiltnis- 
méaBig einfachen Lehren der Mechanik auf bestimmte praktische Fragestellungen anzuwenden. Diese 
auf langjahrige Unterrichtserfahrung aufgebaute Aufgabensammlung bringt deshalb die Lésungen der 
— grdBtenteils aus der technischen Praxis stammenden — Aufgaben in groBer Ausfiihrlichkeit, so daB 
auch der Gurchschnittlich Begabte oder der fiir Theorie weniger Veranlagte zu ae vermag, selbst dann, 
wenn die Aufgabe héhere Anspriiche an ihn stellt. 


Besonderes Gewicht wurde auf die Anwendung gewisser Satze und Methoden der Technischen Mechanik - 
gelegt, die, wie das Prinzip von d’Alembert die Satze von Gastigliano und Coriolis die Freimachungs- 
und Schnittmethoden, dem Anfanger gréBere Schwierigkeiten zm bereiten pflegen. Um schon den an- 
gehenden Ingenieur mit Problemen bekannt zu machen, die erst in jiingster Zeit in der techischen Praxis 
aufgetaucht sind, wurden in verschiedenen Aufgaben Fragen berihrt, deren Klarung durch die Technische 
Mechanik yu entscheidender Bedeutung fiir die Verwirklichung neuartiger Konstroktionsgedanken ge- 
wesen ist, 
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